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A new  concept  and  design  methods  for  the  radiation 
patterns  of  an  array  antenna  system  are  introduced  and 
developed.  It  is  well  known  that  the  analysis  and  design 
problems  of  the  array  antenna  system  are  analogous  to  those 
of  a discrete-time  signal  and  system.  Through  the  analogy, 
the  utilization  of  the  analytic  signal  concept  for  the  design 
of  an  array  antenna  system  is  fully  exploited  and 
investigated . 

First,  the  equivalent  one-dimensional  and  two-dimensional 
discrete-time  analytic  signals  and  the  associated  Hilbert 
transforms  are  defined  and  derived.  There  are  two  types  of 
two-dimensional  discrete-time  analytic  signals  that  are  found 
to  be  useful  for  the  planar  array  antenna  system. 
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Using  the  spectrum  property  as  well  as  the  spectrum 
relations  between  the  real  part  and  the  imaginary  part  of  an 
analytic  signal,  these  design  concepts,  are  theoretically 
discussed  and  physically  evaluated  in  order  to  improve  the 
beamwidth  resolution  of  a conventional  broadside  symmetric 
radiation  pattern  and  to  implement  the  sum  and  difference 
patterns  for  a monopulse  array  antenna  radar  system.  The 
achievements  depend  on  and  are  determined  by  the  realizable 
Hilbert  transformer. 

This  analytic  signal  concept,  along  with  the  modulation 
theory,  is  then  used  to  design  various  desirable  asymmetric 
sidelobe  radiation  patterns  which  include  linear  and  planar 
broadside  narrow  beamwidth  patterns  and  shaped  beam  patterns. 
By  comparing  the  characteristics  of  the  asymmetric  sidelobe 
patterns  to  the  conventional  symmetric  patterns,  we  found  that 
the  asymmetric  patterns  have  a large  advantage.  The  features 
of  those  design  methods  compared  to  the  others  are  conceptual 
simplicity,  computational  effectiveness  and  use  of  analytic 
design  procedures. 

In  summary,  the  analytic  signal  theory  provides  a new 
concept  and  methods  for  designing  various  desirable  radiation 
patterns  of  an  array  antenna  system. 
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CHAPTER  I 


INTRODUCTION 

1 . 1 Motivation  and  Objectives  of  the  Research 

The  fundamental  theory  and  principles  of  an  array  antenna 
system  are  so  widely  known  that  they  can  be  found  everywhere 
in  textbooks  and  published  literature  [1-11].  Since  the  far 
field  radiation  pattern  and  current  excitation  of  an  array 
antenna  system  can  be  formed  as  Fourier  transform  relations 
(discrete-time) , it  has  long  been  recognized  that  there  are 
analogies  in  the  design  problems  between  the  array  antenna 
system  and  the  discrete-time  signal  and  system.  Various  useful 
analysis  and  design  concepts  and  theories  have  been  exploited 
and  beneficially  shared  on  both  areas  [5,7,12-13]. 

The  analytic  signal,  first  introduced  by  Gabor  in  1946, 
has  been  playing  an  important  role  in  signal  processing  and 
communication  systems.  Specifically  speaking,  wherever  the 
Fourier  transform  is  applicable,  the  analytic  signal  concept 
always  provides  both  theoretical  interest  and  practical  useful 
results.  Surprisingly,  no  attention  has  been  paid  to  the  use 
of  the  analytic  signal  concept  in  the  design  of  the  radiation 
patterns  of  an  array  antenna  system.  Such  a design  idea  to  the 
author's  knowledge  is  still  ignored  by  the  experts  in  the 
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antenna  area  and  thus  not  found  in  the  published  literature. 
However,  there  are  several  properties  of  an  analytic  signal 
which  are  well  suited  for  the  design  of  the  radiation  patterns 
of  an  antenna  system.  It  is  the  purpose  of  this  dissertation 
to  exploit  and  exhibit  the  advantages  of  the  design  of  the 
radiation  patterns  of  an  array  antenna  system  using  analytic 
signal  theory. 

1.2  Areas  of  Research 

Basically,  there  are  two  main  trends  in  designing  the 
radiation  patterns  of  an  array  antenna  system:  either  to 
maximize  some  pattern  characteristics  [14-18]  or  to 
approximate  a prescribed  pattern  shape  [19-21].  This 
dissertation  concentrates  on  the  study  of  the  latter  trend  by 
applying  the  analytic  signal  concept.  The  function  of  the 
radiation  pattern  of  an  antenna  system  can  be  alternatively 
interpreted  as  a spatial  filter.  Consequently,  the  required 
shapes  of  the  radiation  patterns  may  very  much  depend  on  the 
applications.  Detection  and  discrimination  of  spatial  targets 
may  require  a narrow  beamwidth  radiation  pattern,  and  tracking 
a moving  target  may  require  an  even  symmetric  and  an  odd 
symmetric  narrow  beamwidth  radiation  pattern.  Some  special- 
purpose  applications  may  require  particular,  shaped  beam 
patterns,  such  as  the  cosecant-squared  shaped  power  pattern 
for  an  air  surveillance  or  ground  mapping  radar  system. 

The  analytic  signal  can  be  formed  by  a Hilbert  transform 
pair  [22],  and  there  are  several  interesting  properties 
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related  to  the  analytic  signal  which  are  well  understood 
[23-26] : 

1.  The  spectrum  of  the  analytic  signal  vanishes  for 
negative  freguency  while  its  dual  (its  complex 
conjugate)  vanishes  for  positive  frequency. 

2.  The  real  and  imaginary  parts  of  the  analytic  signal 
have  an  identical  power  spectrum. 

3.  If  the  real  part  of  the  analytic  signal  is  even 
symmetric,  then  its  imaginary  part  is  odd  symmetric. 

4.  The  analytic  signal  is  completely  described  by  its 
real  or  imaginary  part. 

By  carefully  associating  the  required  shapes  of  the 
radiation  patterns  of  an  antenna  system  with  the  spectral 
properties  of  an  analytic  signal,  we  find  that  Property  1 
indicates  the  beamwidth  resolution  can  doubly  be  increased, 
and  one  of  the  sidelobes  will  be  eliminated  compared  to  a 
typical  broadside  symmetric  antenna  radiation  pattern. 
Properties  2,3,4  provide  a new  idea  for  designing  the 
radiation  patterns  of  a monopulse  radar  system  using  the 
excitations  corresponding  to  the  real  and  the  imaginary  parts 
of  the  analytic  signal.  Those  design  concepts  will  be  fully 
exploited  and  analyzed  in  this  dissertation. 

The  design  of  a broadside  narrowest  beamwidth  radiation 
pattern  with  a specified  symmetric  sidelobe  level  for  an  equal 
interelement  distance  linear  array  antenna  system  has  been 
proposed  by  Dolph  [19]  and  extended  to  continuous  sources  by 
Taylor  [20].  Lopez  [27]  and  Elliott  [28]  have  pointed  out  that 
the  symmetric  sidelobes  of  an  antenna  radiation  pattern  are 
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neither  a requirement,  nor  even  desirable,  because  most 
operation  environments  have  a situation  where  the  effects  of 
the  interferences  (either  intentional  or  unintentional 
clutters)  are  unevenly  distributed  in  the  space.  For  instance, 
a ground  base  radar  system  is  desirable  to  have  a lower 
sidelobe  below  the  horizon  than  above  the  horizon  so  as  to 
minimize  the  effect  of  the  ground  clutter.  Thus,  by  using  an 
asymmetric  sidelobe  radiation  pattern,  there  would  be 
advantages  in  a narrowing  of  the  beamwidth  and  an  increase  in 
gain  relative  to  the  symmetric  case  with  both  sides  of 
sidelobes  having  the  same  sidelobe  level  as  that  below  the 
horizon  of  the  asymmetric  radiation  pattern.  However,  the 
design  of  an  asymmetric  sidelobe  radiation  pattern  is  not  as 
straightforward  as  symmetric  one.  In  [28-29],  Elliott  was  able 
to  generate  an  asymmetric  sidelobe  aperture  radiation  pattern 
by  independently  specifying  peak  sidelobes  on  both  sides  of 
the  Taylor  aperture  pattern.  However,  his  is  a special  design 
case.  Then,  he  has  published  several  papers  [30-34]  related 
to  the  asymmetric  sidelobe  narrow  beamwidth  and  shaped  beam 
radiation  patterns  by  adjusting  the  nulls  of  the  radiation 
pattern  using  the  numerical  iterative  method.  However,  skill 
and  luck  are  often  required  in  the  numerical  manipulations. 

In  this  dissertation,  using  the  analytic  signal  concept, 
we  introduce  a generalized  theory  to  design  the  desirable 
asymmetric  sidelobe  array  radiation  patterns.  The  meaning  of 
generalization  is  that  this  design  method  can  be  equally 
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applied  to  a broadside  narrow  beamwidth  radiation  pattern  and 
a shaped  beam  pattern.  It  will  be  shown  that  this  design 
methodology  appears  to  have  several  advantages  compared  to 
others.  It  is  conceptually  simple  and  straightforward.  It  is 
computationally  inexpensive.  It  provides,  for  the  first  time, 
a design  method,  related  to  an  asymmetric  radiation  pattern, 
which  is  supported  by  a generalized  theory,  implying  that  this 
design  method  can  be  analytically  carried  out. 

One  of  the  most  widely  used  antenna  radiation  patterns 
is  the  shaped  beam  pattern.  Various  techniques  have  been 
suggested  for  synthesis  of  the  shaped  beam  pattern 
[13,21,35-42].  However,  the  equiripple  (minmax)  design 
criterion  (Dolph-Chebyshev  synthesis)  has  never  been  applied 
to  the  synthesis  of  an  asymmetric  shaped  beam  pattern,  because 
of  the  complex  polynomial  representation  of  such  a pattern. 
In  this  research,  by  applying  the  analytic  signal  concept,  we 
will  show  the  feasibility  of  applying  the  equiripple  design 
criterion  to  the  design  of  asymmetric  shaped  beam  patterns, 
although  the  realized  asymmetric  shaped  beam  patterns  may  not 
have  the  equiripple  result  because  of  the  operation  of  the 
Hilbert  transformation. 

Various  practical  examples  are  chosen  and  simulated  to 
demonstrate  the  applications  and  to  show  the  feasibility  of 
the  proposed  design  methodology.  It  will  be  shown  that  the 
contributions  of  the  analytic  signal  concept  to  the  design  of 
array  radiation  patterns  are  practical  and  valuable. 
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1 . 3 Organization  of  the  Dissertation 

In  order  to  equally  apply  the  analytic  signal  concept  to 
the  design  of  the  linear  and  planar  array  antenna  systems, 
Chapter  II  is  devoted  to  defining  the  equivalent  one- 
dimensional and  two-dimensional  discrete-time  analytic  signals 
by  way  of  sampling  theory  and  variable  transformation,  and  to 
deriving  the  associated  required  mathematical  tool — two- 
dimensional  discrete-time  Hilbert  transforms.  In  Chapter  III, 
the  broadside  narrow  beamwidth  array  radiation  patterns  using 
the  excitations  corresponding  to  both  theoretical  and 
realizable  analytic  signals  are  derived  and  quantitatively 
analyzed.  Several  properties  of  the  broadside  narrow  beamwidth 
radiation  patterns  using  the  analytic  signal  concept  are 
discussed.  In  Chapter  IV,  the  general  theory  and  procedures 
for  designing  a broadside  narrow  beamwidth  asymmetric  sidelobe 
radiation  patterns  of  an  array  antenna  system  using  the 
analytic  signal  concept  are  derived.  Examples  are  given  to 
demonstrate  the  introduced  design  method  and  to  compare  the 
performance  between  the  asymmetric  sidelobe  radiation  pattern 
and  the  typical  broadside  symmetric  radiation  pattern. 
Chapter  V extends  the  design  concepts  and  modifies  the  design 
procedures  presented  in  Chapter  IV  to  design  an  asymmetric 
shaped  beam  pattern.  In  Chapter  VI,  the  association  of  the 
generation  of  radiation  patterns  for  a monopulse  array  antenna 
radar  system  with  the  analytic  signal  concept  is  addressed. 
Then,  performance  is  analyzed  and  design  procedure  is 
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formalized  for  the  generation  of  the  radiation  patterns  of  a 
monopulse  array  antenna  radar  system  using  the  analytic  signal 
concept.  Finally,  conclusions  are  presented  in  Chapter  VII. 


CHAPTER  II 


THE  DISCRETE-TIME  ANALYTIC  SIGNAL 
AND  HILBERT  TRANSFORMATION 


2 . 1 Introduction 

The  fundamental  properties  and  applications  of  the  analog 
analytic  signal  and  Hilbert  transform  are  well  developed 
[22-26,43-54].  The  comparable  one-dimensional  discrete-time 
analytic  signal  and  the  associated  Hilbert  transform  can  be 
found  in  [55-57].  In  order  to  apply  the  analytic  signal 
concept  to  the  planar  array  antenna  system,  the  purposes  of 
this  chapter  are  to  deduce  the  two-dimensional  (hereafter 
written  as  2-D)  discrete-time  analytic  signals  and  to  derive 
the  associated  Hilbert  transforms.  First,  the  discrete-time 
analytic  signals  will  be  defined  for  both  1-D  and  2-D  cases. 
Then,  it  will  be  shown  that  there  are  two  types  of  2-D 
discrete-time  analytic  signals  and  Hilbert  transforms  which 
provide  both  theoretical  interest  and  practical  usefulness, 
and  thus  will  be  applied  in  later  chapters. 

From  the  application  point  of  view,  the  spectrum  property 
of  the  discrete-time  analytic  signal,  which  is  directly 
related  to  the  radiation  pattern  of  an  array  antenna  system, 
is  emphasized.  There  is  no  intention  to  exhaustively  discuss 
all  possible  properties  of  such  analytic  signals. 
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2 . 2 Review  of  the  One-Dimensional  Analytic  Signal  Theory  and 
Hilbert  Transform 

An  analytic  function1  is  defined  as  a function  where  all 
derivatives  exist  and  are  continuous  over  a region. 
Furthermore,  a complex  function  f (x+jy) =u (x,y) +jv (x,y)  is 
analytic  in  a complex  region  R,  if  the  partial  derivatives  of 
its  real  part  u(x,y)  and  imaginary  part  v(x,y)  are  continuous 
over  the  region  R and  satisfy  the  following  Cauchy-Riemann 
conditions  [58-59] 

3u  3v  3u  3v 

= and  = - (2.1) 

<9x  3y  dy  3x 

For  a real  continuous  signal  zr(t) , an  analog  analytic 
signal  z(t)  can  be  formed  by  [22] 


z(t)  = zr(t)  + j z 1 ( t ) 


(2.2) 


where  z^t)  is  the  analog  Hilbert  tranform  of  zr(t)  and  the 
Hilbert  transform  is  given  by 

A 1 

Zi(t)  = zr(t)  = zr(t)*  h(t)  = zr(t)  * 

7rt 


n 


r(£) 


(t-O 


de 


(2.3) 


where 


!Later  in  Definition  2.1  the  term  analytic  function  (note 
tilde  used)  will  be  defined.  The  analytic  function  is  not 
analytic  since  its  derivatives  do  not  exist  and,  conseguently , 
it  does  not  satisfy  the  Cauchy-Riemann  conditions  given  by 
(2.1) . 
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A 

z ( • ) represents  the  Hilbert  transform  of  z ( • ) . 

* indicates  the  convolution  operation. 
h(t)  = (l/7rt)  is  the  impulse  response  of  an  analog 
Hilbert  transform. 

pJc<->  denotes  the  Cauchy  principal  value. 


From  Equation  (2.3),  the  Hilbert  transform  of  a 
continuous  signal  can  be  obtained  from  the  residue  values 
evaluated  on  the  real  axis,  provided  that  the  signal  is 
analytic  (derivatives  exist)  over  upper  half  complex  plane. 
The  spectrum  H(w),  the  Fourier  transform  of  h(t),  is 


H(w) 


00 

h(t)e“^a,t  dt  = -jsgn(w) 
—oo 


where 


(2.4) 


sgn (w) 


1 

-1 


cd  > 0. 
cj  < 0 . 


(2.5) 


It  is  well  known  that  the  analytic  signal  z(t)  has 
positive-sided  spectrum  (i.e.  Z(u)=0  for  Ct><0)  . A dual  analytic 
signal  which  is  formed  as  zd (t)  =zr (t) -j z1  (t)  has  negative-sided 
spectrum  (i.e.  Z(w)=0  for  w>0) . 

In  the  discrete-time  signals  domain,  the  analyticity 
concept  defined  in  the  continuous  signals  will  be  meaningless. 
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since  the  discrete-time  signals  are  mathematically  represented 
as  sequences  of  numbers,  and  their  spectra  are  inherently 
periodic.  However,  in  order  to  equally  share  the  analytic 
terminology  in  the  complex  discrete-time  signal,  we  can  first 
generate  a complex  discrete-time  signal  by  sampling  an  analog 
analytic  signal,  then  from  the  property  of  the  spectrum  of 
such  a complex  sampled  analytic  signal,  a general  discrete- 
time analytic  signal  and  the  concept  of  a periodic  half-sided 
spectrum  could  be  appropriately  defined. 

Let's  assume  a bandlimited  analog  analytic  signal 
z (t)  =zr (t)  +j zi  (t)  , which  has  a maximum  spectrum  content  wm,  is 
sampled  with  an  impulse  (impulse  sampling  [60]);  the  sampling 
rate  is  ws=(27r/At)  and  ws  > 2wm.  At  is  the  sampling  interval. 
The  sampled  analog  analytic  signal,  denoted  as  zs(t),  is  [60] 

00 

zs(t)  = E [zr(kAt)+jzi  (kAt)  ] S (t-kAt)  (2.6) 

k=-°o 

Using  the  sampling  theory,  the  spectrum  Zs(w)  of  the  sampled 
z(t)  is 

00 

Zs(cd)  = (0>s/27T)  E Z(w-kws)  (2.7) 

k=-oo 


where 


Z(w) 


00 

z (t) dt 

-00 
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Next,  by  taking  the  Fourier  transformation  on  both  sides  of 
Equation  (2.6),  an  alternative  spectrum  representation  for 
zs(t)  is 


Zs(«) 


'oo  00 

[ s {[zr(kAt)+jz1(kAt)  ]«(t-kAt)  }]e-iut  dt 

—oo  k=-oo 


00 


= Z {[zr(kAt)+jz,(kAt)][ 

k=-oo 


'oo 

S (t-kAt)  e“^utdt] } 

—00 


00 

= Z [zr(kAt)+jzi(kAt)  ]e"jwkAt  (2.8) 

k=-<» 

That  is,  both  Equation  (2.7)  and  Equation  (2.8)  should 
represent  the  same  spectrum  of  the  sampled  z(t),  which  is 
periodic  with  periodicity  ios.  Let's  take  a linear 
transformation  of  the  variables  from  w to  u with  u=wAt;  that 
is,  we  map  the  spectrum  from  -( us/2 ) to  (ws/2)  in  u into  -n 
to  n in  u,  and  the  mapping  procedure  is  periodic  with  ws  in 
w and  2n  in  u.  After  changing  the  variable,  Equation  (2.8) 
becomes 


00 

zs(u)«  Z [zr(kAt)+jZl(kAt)  ]e"juk  (2.9) 

k=-oo 

Now,  the  spectrum  in  Equation  (2.9)  is  a periodic  funtion  of 
u with  periodicity  2n . If  the  spectrum  Z(w)  of  the  original 
analog  analytic  signal  is  as  depicted  in  Figure  2.1(a),  the 
spectrum  Zs(w)  of  the  sampled  analog  analytic  signal  in  either 
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Z(«) 


-w 


m 


W 


m 


(a) 


Zs(«) 


Zs(u) 


Figure  2.1  The  spectrum  transformation  between  the  analog 
and  discrete-time  analytic  signals. 

(a)  The  spectrum  of  an  analog  analytic  signal. 

(b)  The  spectrum  of  the  sampled  analytic 
signal  of  (a) . 

(c)  The  spectrum  of  a discrete-time  analytic 
signal  transformed  from  (b) . 


Equation  (2.7)  or  Equation  (2.8)  (because  both  of  them 
represent  the  same  spectrum)  and  the  variable  transformed 
spectrum  in  Equation  (2.9)  will  be  shown  as  in  Figure  2.1(b) 
and  Figure  2.1(c)  respectively.  As  seen  from  Figure  2.1,  the 
original  single-sided  bandlimited  spectrum  Z(w)  of  an  analog 


14 


analytic  signal  is  mapped  into  a periodic  half-sided  spectrum 
Zs(u),  by  ways  of  sampling  and  variable  transformation.  That 
is,  the  single-sided  spectrum  Z(w)=0  for  w<0  is  mapped  into 
the  periodic  half-sided  spectrum  Zs(u)=0  for  -7r<u<0  of  each 
period,  provided  that  ws  > 2wm. 

In  essence,  the  sampled  signal  [ zr  (kAt) +j zi  (kAt)  ] in 
Equation  (2.9)  will  be  defined  as  a discrete-time  analytic 
signal.  Consequently,  the  transformation  format  in  Equation 

(2.9)  can  be  regarded  as  the  Fourier  transform  for  a discrete- 
time signal,  which  is  formally  stated  as  follows. 

The  Fourier  transformation  of  an  absolutely  summable 
discrete-time  signal  x(n)  is  defined  as  [56,p.45] 

00 

X (u)  = £ x(n)e~jun  (2.10) 

n=-oo 

Hereafter,  X(u)  is  called  the  spectrum  of  the  discrete-time 
signal  x(n).  Note  that,  in  Equation  (2.10),  the  function 
variable  n is  an  integer  and  the  spectrum  variable  u is 
continuous.  The  spectrum  X(u)  of  a general  discrete-time 
signal  is  inherently  periodic  with  periodicity  2n . For  a 
finite  length  discrete-time  signal  x(n)  with  length  N,  the 
spectrum  calculated  using  the  discrete  Fourier  transform 
(DFT)  [60],  where  the  resultant  spectrum  is  both  discrete  and 
periodic,  is  the  sampled  version  of  the  spectrum  in  Equation 

(2.10)  with  sampling  period  (27r/N). 
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Since  the  spectrum  is  a periodic  function  of  u.  Equation 
(2.10)  can  be  viewed  as  a Fourier  series  representation  of 
X(u)  . According  to  the  theory  of  the  Fourier  series,  the 
discrete-time  signal  x(n)  can  be  derived  from  x(u)  as 


x(n) 


1 f7r 

X(u)e^un  du 

2n  J -n 

v 


(2.11) 


Equations  (2.10)  and  (2.11)  thus  form  a transform  pair  which 
is  commonly  defined  as  the  discrete-time  Fourier  transform 
pair,  and  denoted  as  x(n)«X(u) . Since  we  are  dealing  with  a 
discrete  system — array  antenna  system,  hereafter,  any 
statement  of  relations  between  the  spectrum  and  the  discrete- 
time signal  will  be  directly  related  to  the  Fourier  transform 
pair  of  Equations  (2.10)  and  (2.11). 

From  the  illustrated  periodic  half-sided  spectrum  of  the 
sampled  analog  analytic  signal  in  Figure  2.1(c) , we  then  could 
define  an  equivalent  discrete-time  analytic  signal  by  the 
notion  that  a complex  discrete-time  signal  which  has  the 
periodic  half-sided  spectrum  is  the  discrete-time  analytic 
signal,  and  neglect  the  analyticity  requirement  in  the 
continuous  case.  Thus,  let's  formally  define  the  periodic 
half-sided  spectrum  (due  to  the  periodic  property,  we  only 
need  to  explicitly  define  the  -n  to  ir  interval)  and  a 
discrete-time  analytic  signal  as  follows. 

Definition  2.1  A spectrum  X(u)  of  a discrete-time  signal 
is  said  to  be  periodic  half-sided,  if  X(u)=0  for  -n  < u < 0, 
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with  period  of  2ir  (refer  to  Figure  2.1(c)).  A complex 
discrete-time  signal  z(n)  is  defined  as  an  analytic  signal, 
if  its  spectrum  Z(u)  is  periodic  half-sided.2 

After  establishing  the  meaning  of  the  discrete-time 
analytic  signal,  the  impulse  response  h(n)  of  the  discrete- 
time Hilbert  transform  can  be  derived  from  the  spectrum 
property  of  the  analytic  signal  [56].  Let's  denote  the  defined 
half-sided  periodic  spectrum  as  Z(u),  and  its  corresponding 
analytic  signal  as  z(n).  Based  on  the  property  of  the 
discrete-time  Fourier  transform,  z(n)  should  be  complex  and 
can  be  alternatively  expressed  as 


z (n)  = zr(n)  + j z1  (n) 


(2.12) 


where  both  zr(n)  and  z^n)  are  real  and  represent  the  real  and 
imaginary  part  of  z(n),  respectively.  It  is  easily  shown  that 
zr(n)  = [ (1/2)  (z (n)+z* (n) ) ] and  zi  (n)  = [ (-j/2) (z (n) -z* (n) ) ] , where 
z*(-)  represents  the  complex  conjugate  of  z().  With  Zr(u)  and 
Zj(u)  denoting  the  spectra  of  zr(n)  and  z^n),  it  is 
straightforward  to  show  that 


Zr(u)  = (1/2) [ Z (u)  + Z*(-u) ] 

Zi(u)  = (~j/2) [Z (u)  - Z*(-u)]  (2.13) 


2Henceforth  we  will  use  the  "analytic  signal"  (with  a 
tilde  on  the  first  letter  of  the  word  "analytic")  to  denote 
the  discrete-time  analytic  signal  whenever  it  is  necessary 
to  clearly  distinguish  it  from  "analytic  signal"  which  is  used 
in  the  continuous  case.  Furthermore,  there  are  two  kinds  of 
analytic  signals  related  to  the  positive  (X(u)=0  for  -7r<u<0) 
and  negative  (where  X(u)=0  for  0<u<7r)  periodic  half-sided 
spectra. 
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Applying  the  half-sided  spectrum  property  of  Z(u),  the 
relations  among  Z(u),  Zr(u),  Z^u)  are 

Z(u)  = 2Zr(u)  = ZjZ^u)  for  0 < u < n 

Z*(-u)=  2Zr(u)  = -2jZi(u)  for  -n  < u < 0 (2.14) 

Thus,  from  Eguation  (2.14),  Zr(u)  and  Z^u)  are  related  as 
follows 


Zi(u)  = H(u)Zr(u) 

Zr(u)  = -H (u) Zi (u)  (2.15) 

where  H(u)  is  called  the  discrete-time  Hilbert  transformer  and 
defined  as 


H (u) = -jsgn(u) 


and 


sgn (u) = 


-1 

1 


-7T  < U < 0. 
0 < U < 7T 


(2.16) 


(2.17) 


Then,  the  impulse  response  h(n)  of  the  discrete-time 
Hilbert  transform  is  the  Inverse  Fourier  transform  of  H(u)  in 
Eguation  (2.16),  that  is 


h(n) 


1 IV 

H (u) e ^un  du 

2n  -7r 

v 


After  a straightforward  integration,  h(n)  becomes 
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h (n)  = 


1 

(l-cosn*)  , 

nn 

0 , 


if  n M 
if  n = 0 


(2.18) 


which  is  odd  symmetric  sequence  about  the  origin  and  depicted 
in  Figure  2.2. 


h(n) 


(2  /*) 


• -5  -4  -3  -2  -1 

* i * | ' r 

-(2/5*)  i 


1 


-(2/3*) 


(2/3*) 

■ ^ 


(2/5*) 

1 


2 3 4 5 


n 


- (2/*) 


Figure  2.2  The  impulse  response  of  one-dimensional 
discrete-time  Hilbert  transform. 


Using  Equation  (2.15),  the  relationship  between  the  real 
and  the  imaginary  parts  of  an  analytic  signal  z(n)  in  Equation 
(2.12)  can  be  represented  by  the  following  discrete-time 
convolution  pair, 


00 

zi  (n)  = E zr (k) h (n-k)  (2.19) 

k=-oo 


00 

zr(n)  = - E zi  (k) h (n-k) 

k=-oo 


(2.20) 
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Equations  (2.19)  and  (2.20)  form  the  desired  one-dimensional 

discrete-time  Hilbert  transform  relations  between  the  real  and 

the  imaginary  parts  of  an  analytic  signal.  It  should  be  noted 

that,  for  any  real  discrete-time  signal,  whether  finite  length 

or  not,  the  corresponding  analytic  signal  is  unrealizable 

since  it  always  has  an  infinite  sequence.  This  is  due  to  the 

fact  that  the  impulse  response  of  the  discrete-time  Hilbert 

transform  in  Equation  (2.18)  is  an  infinite  sequence.  Some  of 

the  properties  related  to  the  discrete-time  Hilbert  transform 

can  be  found  in  [56],  and  those  related  to  the  design  of  the 

radiation  pattern  of  an  array  antenna  will  be  stated  and 

proved  later  whenever  it  becomes  necessary  and  appropriate. 

2 . 3 The  2— D Discrete-Time  Analytic  Signal  and  2-D 
Discrete-Time  Hilbert  Transform 

In  the  2-D  signal  theory,  there  appears  to  be  no  defined 
analog  analytic  signal,  and  so  we  could  not  obtain  the 
meanning  of  a 2-D  discrete-time  analytic  signal  from  a 
comparable  2-D  continuous  analytic  signal.  We  will, 
nevertheless,  propose  a definition,  by  maintaining  the 
physical  meaning  and  the  terminologies  used  in  the  1-D  case 
and  meeting  the  application  requirement  of  finding  a 2-D 
complex  discrete-time  signal  whose  spectrum  is  periodic 
quadrant-plane  or  periodic  half-plane3  as  will  be  defined 


3Both  quadrant-plane  spectrum  and  half-plane  spectrum  are 
useful  concepts  and  properties  for  the  design  of  the  radiation 
pattern  of  a planar  array  antenna  system. 
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later.  Consequently,  we  will  define  a 2-D  complex  discrete- 
time signal,  which  possesses  the  spectrum  property,  as  an 
analytic  signal  and  derive  the  associated  Hilbert  transform 
which  will  generate  the  desired  2-D  analytic  signal. 

First,  the  Fourier  transform  of  an  absolutely  summable 
2-D  discrete-time  signal  x(m,n)  is  defined  as  [61,p.32] 

00  00 

X(u,v)=  2 2 x(m,n) e~i (um+vn)  (2.21) 
m=-°o  n=-oo 

X(u,v)  is  called  the  spectrum  of  the  2-D  discrete-time  signal 
x(m,n).  Note  that  the  spectrum  of  a 2-D  discrete  signal  is 
generally  continuous  and  periodic  in  both  spectrum  variables 
u and  v with  period  of  2n . 

Since  Equation  (2.21)  is  similar  to  the  Fourier  series 
representation  for  X(u,v) , the  2-D  discrete-time  signal  x(m,n) 
can  be  derived  from  X(u,v)  by  using  the  orthogonal  property 
of  the  transformation  kernel  (over  2n  period)  in  Equation 
(2.21) ; that  is, 

X(u,v) (um+vn)  dudv  (2.22) 

Thus,  Equations  (2.21)  and  (2.22)  form  a transformation  pair 
which  is  defined  as  the  2-D  discrete-time  Fourier  transform 
pair,  and  denoted  as  x (m,  n)  ++X  (u,  v)  . 


x (m, n) = 


An t‘ 


n 

— 7T 
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2.3.1  The  definition  of  2-D  discrete-time  analytic  signals 
As  mentioned  earlier,  by  maintaining  the  concepts  used 
in  a 1-D  case  and  meeting  the  application  reguirement,  a 2-D 
discrete-time  analytic  signal  can  be  defined.  However,  one 
generalization  is  formulated  by  requiring  that  the  spectrum 
be  zero  over  some  specific  region  within  each  period  of  the 
2-D  spectrum,  rather  than  half-side  interval  within  each 
period  of  the  spectrum  in  the  1-D  discrete-time  case.  Thus, 
there  are  two  2-D  discrete-time  analytic  signals.  Both  of  them 
and  their  associated  spectra  are  defined  as  follows. 

Definition  2.2  A spectrum  X(u,v)  of  a 2-D  discrete-time 
signal  is  said  to  be  periodic  quadrant-plane,  if  X(u,v)  is 
X(u,v)  = 0 for  ( — 7r  < u < 0 or  -7r  < v < 0) 

with  period  of  2n  in  both  u and  v.  A typical  rectangular 
periodic  quadrant-plane  spectrum  is  depicted  in  Figure  2.3(a) . 
Consequently,  a 2-D  complex  discrete-time  signal  z(m,n)  is 
defined  as  a 2-D  analytic  signal — Type  I.  if  its  spectrum 
Z(u,v)  is  periodic  quadrant-plane.4 

Definition  2.3  A spectrum  X(u,v)  of  a 2-D  discrete-time 
signal  is  said  to  be  periodic  half-plane,  if  X(u,v)  is 
X(u,v)  = 0 for  {-it  < u < 0) 

with  period  of  2n  in  both  u and  v.  A typical  rectangular 
periodic  half-plane  spectrum  is  depicted  in  Figure  2.3(b). 
Consequently,  a complex  2-D  discrete-time  signal  z(m,n)  is 

4There  are  four  possible  2-D  Type  I analytic  signals 
corresponding  to  the  four  periodic  quadrant-plane  spectra. 
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defined  as  a 2-D  analytic  signal — Type  II.  if  its  spectrum 
Z(u,v)  is  periodic  half-plane.5 
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Figure  2.3  The  2-D  quadrant-plane  and  half-plane  spectra. 

(a) The  2-D  quadrant-plane  spectrum  (Type  I). 

(b) The  2-D  half-plane  spectrum  (Type  II). 


2.3.2  The  definition  of  2-D  discrete-time  Hilbert  transforms 
A generalized  one-dimensional  discrete-time  Hilbert 
transform  which  is  capable  of  generating  an  analytic  signal 
with  half-sided  spectrum  is  indicated  in  Equation  (2.19). 
Using  the  same  token,  a generalized  2-D  discrete-time  Hilbert 
transform  of  a real  discrete-time  signal  zr(m,n)  can  be 
expressed  as  the  following  2-D  discrete-time  convolution. 


zi(m,n)  = zr(m,n)  **h(m,n) 


00  00 

= 2 2 zr  (€  , k)  h (m-£  , n-k)  , for  -oo  < m,n  < °o 

€=-oo  k=-oo 

(2.23) 


5There  are  four  possible  2-D  Type  II  analytic  signals 
corresponding  to  the  four  periodic  half-plane  spectra. 
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where  z^(m,n)  denotes  the  2-D  discrete-time  Hilbert  transform 
of  zr(m,n),  **  represents  the  2-D  discrete-time  convolution, 
and  h(m,n)  is  the  impulse  response  of  the  2-D  discrete-time 
Hilbert  transform. 

Applying  the  2-D  discrete-time  convolution  theorem  to 
Equation  (2.23),  the  generalized  spectrum  relation  of  the 
Hilbert  transformation  is 

Zf  (u,v)  = H(u,v)Zr(u,v)  (2.24) 

where  Zi(u,v),  Zr(u,v)  and  H(u,v)  are  the  spectrums  of  the 
2-D  discrete-time  signal  z^ir^n),  zr(m,n)  and  the  impulse 
response  h(m,n),  respectively. 

However,  in  order  to  be  consistent  with  the  physical 
meaning  of  a Hilbert  transformer,  the  2-D  discrete-time 
Hilbert  transforms  should  be  further  restrictedly  defined  as 
follows:  the  2-D  discrete-time  convolution  in  Equation  (2.23) 
can  be  defined  as  the  2-D  discrete-time  Hilbert  transform  of 
the  real  signal  zr(m,n),  if  the  spectrum,  H(u,v),  of  the 
impulse  response  h(m,n)  in  Equation  (2.23)  has  the  format 
as  H(u,v)=-jsgn(u,v) , where  sgn(u,v)  represents  the  2-D 
version  of  sgn  function,  that  is,  depending  on  u and  v,  its 
value  is  always  either  1 or  0 in  the  2-D  plane  (reference  to 
Equation  (2.40)).  Consequently,  zr(m,n)  and  zi(m,n)  in  Equation 
(2.23)  are  said  to  be  a 2-D  discrete-time  Hilbert  transform 
pair.  Based  on  the  definition,  it  will  be  shown  that  there  are 


24 


two  types  of  2-D  discrete-time  Hilbert  transforms  related  to 
the  defined  Type  I and  Type  II  analytic  signals  in  the  2-D 
case. 

Having  defined  the  2-D  analytic  signals  and  the  format 
of  2-D  discrete-time  Hilbert  transform,  we  first  show  that  no 
2-D  Type  I analytic  signal  (i.e.  the  periodic  quadrant-plane 
spectrum)  can  be  derived  in  general  using  H(u,v)=-jsgn(u,v) . 
However,  after  a complete  consideration  and  exploration,  there 
are  three  methods  that  can  be  used  to  generate  the  2-D 
analytic  signals  (either  Type  I or  Type  II)  which  have  these 
desired  quadrant-plane  and  half-plane  spectra.  These  three 
methods  are  stated  in  the  following  cases,  and  their 
associated  discrete-time  Hilbert  transforms  will  be  discussed 
and  derived  in  Section  2.3.4  through  Section  2.3.6. 

Case  la.  It  is  possible  to  obtain  a 2-D  Type  I analytic 
signal  using  two  stages  of  transformations. 

Case  lb.  It  is  possible  to  obtain  a 2-D  Type  I analytic 
signal  using  two  1-D  analytic  signals. 

Case  II.  It  is  possible  to  obtain  a 2-D  Type  II  analytic 

signal  using  a 2-D  discrete-time  Hilbert  transform. 

2.3.3  The  Properties  of  the  general  expression  of  2-D 
discrete-time  Hilbert  transform 

Formally,  we  can  directly  derive  the  2-D  impulse 
response  h(m,n)  from  the  general  expression  of  the  2-D 
discrete-time  Hilbert  transform  in  Equation  (2.23),  and 
explore  the  general  properties  of  the  2-D  discrete-time 
Hilbert  transform. 
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2. 3. 3.1  The  impulse  response  of  the  general  expression  of 
2-D  discrete-time  Hilbert  transform 

Let's  first  try  to  derive  a generalized  2-D  Hilbert 
transform,  hg(m,n),  and  then  examine  its  spectral  properties 
where  Hg(u,v)  is  the  2-D  Fourier  transform  of  hg(m,n)  . The 
derivation  of  Hg(u,v)  can  be  carried  out  as  follows:  let  the 
kernel  of  the  2-D  discrete-time  Hilbert  transform  in  Equation 
(2.23)  be  separable  (just  like  the  separability  of  the  kernel 
of  the  2-D  Fourier  transformation) ; that  is, 

hg  (m-£ , n-k) =h (m-£ ) h (n-k) . 

where  h(n)  is  defined  in  Equation  (2.18).  Thus,  the  discrete 
convolution  in  Equation  (2.23)  becomes 


z,  (m,  n)  =2  Z zr  ( £ , k)  h (m-£ ) h (n-k)  , 
£ k 


-oo  < m,n,£,k  < °o 

(2.25) 


By  taking  the  2-D  discrete-time  Fourier  transformation  on  both 
sides  of  Equation  (2.25),  the  spectrum  of  zi(m,n),  denoted  as 
Zi(u,v),  is 


Zi(u,v)=  Z Z [2  Z zr(£,k)h(m-£)h(n-k)  ]e“j  (um+vn) 
m n £ k 

= Z Z [Z  Z zr(£,k)h(m-£)h(n-k) (um+vn) 

£ k m n 

= Z Z zr ( £ , k) { [Z  h(m-£) e~jum] [Z  h (n-k) e~l vn] } 
£ k m n 


(2.26) 

From  Equation  (2.18),  the  impulse  response  h(n)  is  an  odd 
symmetric  function;  that  is, 
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h(n-k)  = h[- (k-n) ] = -h(k-n) 

Then,  one  of  the  discrete  convolutions  of  the  complex 
exponential  e-Jum  in  Equation  (2.26)  becomes 

2 h (m-£ ) e“juin  = - [ £ h(£-m)e~jum]  = - [h  (£ ) ] 

m m 

= - [2  h(m)^um] 

m m 

= [-jsgn(-u) ] = -jsgn(u) (2.27) 

where  sgn(u)  is  defined  in  Equation  (2.17).  Using  the  result 
of  Equation  (2.27),  the  spectrum  Zi(u,v)  in  Equation  (2.26) 
becomes 

zi(u,v)=  2 E zr(£  ,k)  [ jsgn(u)  ] [ jsgn(v)  e-^vk] 

£ k 

= [-sgn(u) sgn(v) ] E E zr (£  , k) e~i (u^+vk) 

£ k 

= [-sgn(u) sgn(v) ]Zr(u,v)  (2.28) 

where  Zr(u,v)  is  the  spectrum  of  zr(m,n)  . Comparing  Equations 
(2.24)  and  (2.28),  the  spectrum  of  the  generalized  2-D 
discrete-time  impulse  response,  denoted  as  Hg(u,v),  is 

Hg(u,v)  = [-sgn(u) sgn(v) ] (2.29) 

Then,  by  taking  the  2-D  discrete-time  inverse  Fourier 
transform  of  Hg(u,v)  in  Equation  (2.29),  the  2-D  discrete-time 
impulse  response,  denoted  as  h (m,n) , is 
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hg(m,n)=  

47T 


7T 

-7T 


7T 

-7T 


[-sgn(u)  sgn(v)  ] ei  (Ultl+vn) 


dudv 


After  some  mathematical  manipulation,  hg(m,n)  becomes 


hg  (m,  n) 


0,  for  m=n=0  or  m=o  or  n=0. 

1 

[cos  (mn)  cos  (nw)  -cos  (m7r)  -cos  ( rnr ) +1]  , 

(mn) n2 

elsewhere. 


(2.30) 


Examining  the  values  of  hg(m,n)  for  all  possible  combinations 
of  m and  n in  Equation  (2.30),  hg(m,n)  becomes 


hg  (m,  n) 


4 

7 ' 

(mn)  n 
0, 


both  m and  n are  odd 

(2.31) 

elsewhere 


From  Equation  (2.31),  the  2-D  impulse  response  derived 
from  the  general  expression  of  2-D  discrete-time  Hilbert 
transform  has  cross-quadrant  even  symmetry  and  adjacent- 
quadrant  odd  symmetry.  The  impulse  response  hg(m,n)  and  the 
spectrum  of  a uniformly  truncated  version  of  hg(m,n)  (where 
the  lengths  of  both  m and  n are  11)  are  plotted  in  Figure 
2.4(a)  and  (b) . 
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(a) 


Figure  2.4  The  impulse  response  and  spectrum  of  hg(m,n). 

(a)  The  impulse  response  hg(m,n). 

(b) The  spectrum  of  uniformly  truncated  version 

of  (a) . (M=N=11) 


29 


2. 3. 3. 2 The  properties  of  the  transformation  using  impulse 
response  hg(m.n) 

Using  the  impulse  response  hg(m,n)  in  Equation  (2.31),  a 
2-D  complex  discrete-time  signal  z(m,n)  can  be  formed  as 

z(m,n)  = zr(m,n)  + jz^n^n)  (2.32) 

where  zr(m,n)  is  a real  2-D  discrete-time  signal  and 
zi  (m,n)  = [zr(m,n)  **hg(m,n)  ] . From  Equation  (2.28),  the  spectrum 
of  zi(m,n),  denoted  as  Z^UjV),  is 

Zi(u,v)  = (-1) sgn(u) sgn(v) Zr(u,v) 

Thus,  the  spectrum  of  the  z(m,n)  in  Equation  (2.32),  denoted 
as  Z(u,v),  becomes 

Z(u,v)  = Zr(u,v)  - j [sgn(u) sgn(v) Zr(u,v)  ] 

= Zr(u,v)[l  - jsgn(u) sgn(v) ] (2.33) 

The  result  shown  in  Equation  (2.33)  indicates  that  the 
spectrum  of  the  2-D  complex  discrete-time  signal  in  Equation 
(2.32)  , which  is  formed  by  using  the  impulse  response  hg(m,n) , 
does  not  possess  the  previously  defined  quadrant-plane  or 
half-plane  spectrum,  and  consequently , the  complex  signal 
zfm.n)  in  Equation  (2.32)  is  not  a 2-D  analytic  signal. 

Again,  using  zr(m,n)  and  z^n^n)  in  Equation  (2.32),  let's 
form  another  2-D  discrete-time  signal  as 


z(m,n)  = zr(m,n)  + (j)2zi(m,n) 
= zr(m,n)  - zi(m,n) 


(2.34) 
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Note  that  z(m,n)  becomes  a real  2-D  discrete-time  signal,  and 
consequently,  is  not  an  analytic  signal.  However,  it  is  easy 
to  show  that  the  spectrum  of  z(m,n)  in  Equation  (2.34), 
denoted  as  Z(u,v),  is 

Z(u,v)  = Zr(u,v)  + sgn(u)  sgn(v)  Zr(u,v) 

= Zr  (u, v) [ 1+sgn (u) sgn (v) ] 

that  is. 


Z (u, v) 


2Zr(u,v) , 

Zr(u,v) , 
0, 


for  (-7r  < u < 0 and  -n  < v < 0)  or 
(0  < u < 7r  and  0 < v < n) 
for  (u=v=0  or  u=0  or  v=0) 

otherwise  (2.35) 


which  is  a cross-quadrant  spectrum  (spectrum  vanishes  in  both 
Quadrant  II  and  IV)  . We  will  define  this  cross-quadrant 
spectrum  as  a Type  III  spectrum. 

It  needs  to  be  pointed  out  that  the  2-D  signal  pair 
zr(m,n)  and  z^mjn)  in  either  Equation  (2.32)  or  Equation 
(2.34)  never  forms  a 2-D  discrete-time  Hilbert  transform  pair. 
This  is  because  the  spectrum,  Hg(u,v),  of  the  2-D  impulse 
response  hg(m,n)  in  Equation  (2.31)  is  not  consistent  with 
the  format  of  H(u, v) =-jsgn(u,v)  for  the  2-D  discrete-time 
Hilbert  transform  which  was  defined  in  Section  2.3.2.  Thus, 
based  on  the  separable  assumption,  derivation  and  discussion 
above,  we  can  establish  a conclusion  in  Lemma  2.1 
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Lemma  2 . 1 In  the  2-D  signal  domain,  there  is  no 
generalized  2-D  discrete-time  Hilbert  transform  which  can 
produce  a 2-D  analytic  signal,  whose  spectrum  is  quadrant- 
plane.  However,  the  2-D  impulse  response  which  is  derived  from 
the  general  expression  of  the  2-D  discrete-time  Hilbert 
transform  is  capable  of  generating  a 2-D  real  discrete-time 
signal,  whose  spectrum  is  cross-quadrant  (i.e.  Type  III) . 


2.3.4  The  generation  of  a 2-D  analytic  signal — Type  I using 
two  stages  of  transformations  (Case  la) 

According  to  Lemma  2.1,  given  a general  real  2-D 
discrete-time  signal,  the  2-D  impulse  response  in  Equation 
(2.31)  is  capable  of  generating  another  real  2-D  discrete-time 
signal  whose  spectrum  is  cross-quadrant  (Type  III) . The 
purposes  of  this  section  are  to  derive  a 2-D  discrete-time 
Hilbert  transform  (Type  I) . so  that  a 2-D  analytic  signal  with 
quadrant-plane  spectrum — Type  I can  be  obtained  from  a special 
real  2-D  discrete-time  signal  whose  spectrum  is  cross-quadrant 
(Type  III) , and  to  devise  a two-stage  transformation,  so  that 
a 2-D  analytic  signal  with  quadrant-plane  spectrum — Type  I can 
be  obtained  from  a general  real  2-D  discrete-time  signal. 

2. 3. 4.1  The  2-D  discrete-time  Hilbert  transform  (Type  I) 

Consider  a special  case  where  a 2-D  real  discrete-time 
signal  zrr(m,n)  has  a cross-quadrant  spectrum  (Type  III)  which 
is  denoted  as  Zrr(u,v)  ; that  is, 


Zrr(u,v)=2  2 zrr(m,n)  e~i  (um+vn)  = 0, 
m n 

for  (0  < u < 7r  and  -ir  < v < 0) 
or  (-7T  < u < 0 and  0 < v < n 
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Then,  let's  form  a complex  signal  z(m,n)  whose  real  part  is 
zrr(m,n)  ; that  is,  z(m,n)=  zrr(m,n)+j  z1{(m,n).  Now,  we  want  to 
deduce  a 2-D  discrete-time  Hilbert  transform  (Type  I)  which 
is  capable  of  generating  z(j(m,n)  from  zrr(m,n),  such  that  the 
spectrum  of  the  complex  signal  z(m,n),  denoted  as  Z(u,v),  is 
a quadrant-plane;  that  is, 

Z(u,v)  = 0,  for  (—7 r < u < 0 or  -n  < v < 0)  (2.36) 

Since  z(m,n)=zrr(m,n)+jzii(m,n)  , it  is  true  that 

zrr(m,n)  = (1/2)  [z(m,n)+z*(m,n)  ] 
jzii(m,n)  = (1/2)  [z  (m,n) -z*  (m,n)  ] 

where  z ( • , • ) is  the  complex  conjugate  of  z ( • , • ) . With 
Zrr(u,v)  and  Zif(u,v)  denoting  the  spectra  of  zrr(m,n)  and 
zii(m,n),  it  can  be  shown  that 

Zrr(u,v)  = (1/2)  [Z  (u,v)  + Z * ( -u , -v ) ] 
jZ^O^v)  = (1/2)  [Z(u,v)  - Z * ( -u , -v ) ] 

Let ' s presume  that  the  cross-quadrant  spectrum  property  of 
Z(u,v)  in  Equation  (2.36)  is  valid.  Then  the  relations  among 
Z(u,v),  Zrr(u,v)  and  Zii(u,v)  are 

Z(u,v)=  2Zrr(u,v)=  2jZ1i(u,v),  (0  < u < 7r  and  0 < v < n) 

(except  u=v=0) 

Z*  (-u, -v)  =2Zrr(u,  v)  =-2  jZi1  (u, v)  , (-7T  < u < 0 and  -7r<v<0) 

(except  u=v=0) 

(2.37) 
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Thus,  from  Equation  (2.37),  the  Zrr(u,v)  and  Zii(u,v)  should 
be  related  as  follows 


2^  (u,  v)  = Hj(u,v)Zrr(u,v) 

Zrr(u,v)  = -HI(u,v)Zii(u,v)  (2.38) 

where 


Hj(u,v)  = -jsgnj  (u,  v) 


(2.39) 


and 


sgnj(u,v) 


0 < u < n and  0 < v < n 
-7T  < u < 0 and  -7r  < v < 0 
u=v=0  and  elsewhere  (2.40) 


Then  the  impulse  response  of  the  2-D  discrete-time  Hilbert 
transform  (Type  I),  denoted  as  hj(m,n),  can  be  obtained  by 
taking  the  2-D  discrete-time  inverse  Fourier  transform  of 
Hj  (u, v)  in  Equation  (2.39), 


hj(m,n)  = — { 


47T 


'0 

-7T 


ro 


j e j (um+vn)  dudv 


-7T 


7T 


jej (um+vn)  dudv) 


After  taking  the  integration  with  all  possible  combinations 
of  m and  n,  the  impulse  response  hj(m,n)  is 
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hj  (m,  n) 


1-cos (nw) 
2nn 

1-cos  (m7r) 
2mn 
0, 


m = 0 and  n f 0 

(2.41) 

n = 0 and  in  t 0 
otherwise 


The  impulse  response  hj(m,n)  in  Equation  (2.41)  and  the 
spectrum  Hj(u,v)  of  a uniformly  truncated  version  (the  lengths 
of  both  m and  n are  11)  of  h,(m,n)  are  plotted  in  Figure 
2.5(a)  and  (b) . 

From  Equation  (2.38),  the  relations  between  zrr(m,n)  and 
zii(m,n)  are 

zii  (m, n)  = zrr(m,n)  **hj  (m,n) 

zrr(m,n)  = — zn  (m,n)  **h,  (m,n)  (2.42) 

It  is  straightforward  to  prove  that  if  the  real  signal 
zrr(m,n)  has  cross-quadrant  spectrum  (Type  III) . then  the  2-D 
complex  signal  z(m,n)=zrr(m(n)+jz1i(m,n)  has  quadrant-plane 
spectrum  (Type  I ) . Thus,  hj(m,n)  in  Equation  (2.41)  is  the 
desired  impulse  response  for  the  2-D  discrete-time  Hilbert 
transform  (Type  I)  . Furthermore,  by  definition,  the  2-D 
discrete-time  convolution  in  Equation  (2.42)  forms  the  desired 
2-D  discrete-time  Hilbert  transform  relations  (Type  I)  between 
the  signal  pair  zrr(m,n)  and  z1i(m,n). 


Spectrum  C dB) 
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Figure  2.5  The  2-D  discrete-time  Hilbert  transform  (Type  I). 

(a) The  impulse  response  hj(m,n). 

(b) The  spectrum  of  uniformly  truncated  version 
of  (a) . (M=N=11) 
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2. 3. 4. 2 The  generation  of  a 2-D  analytic  signal — Type  I 

After  deriving  the  two  required  2-D  transformations  in 
Section  2.3.3  and  Section  2. 3. 4.1,  the  two-stage 
transformation  and  the  procedures  used  to  produce  a 2-D 
analytic  signal  with  quadrant-plane  spectrum  (Type  I)  from  a 
real  2-D  signal  zr(m,n)  can  be  depicted  in  Figure  2.6. 


z (m,  n) 


2 : Summation 
X : Multiplication 

**:2-D  discrete-time  convolution 
Figure  2.6  A general  2-D  analytic  signal — Type  I generator. 

In  Figure  2.6,  hg(m,n)  and  hj(m,n)  denote  the  2-D  impulse 
responses  in  Equation  (2.31)  and  Equation  (2.41).  zr(m,n)  is 
a real  2-D  discrete-time  signal,  and  the  2-D  signal  zrr(m,n) 
is  formed  as 


zrr(m,n)  = zr(m,n)  - [ zr  (m,  n)  **hg  (m,  n)  ] (2.43) 

According  to  Equation  (2.34)  and  (2.35),  the  spectrum  of 
zrr(m,n),  denoted  as  Zrr(u,v),  is  cross-quadrant;  that  is. 


Zrr(u,v) 


2Zr(u,v)  , for  ( — 7r  < u <0  and  -n  < v < 0)  or 
(0  < u < 7r  and  0 < v < n) 
Zr(u,v) , for  (u=v=0  or  u=0  or  v=0) 

0,  otherwise 


(2.44) 
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where  Zr(u,v)  is  the  spectrum  of  zr(m,n).  Note  that  zrr(m,n) 
is  a real  2-D  signal  and  has  cross-quadrant  spectrum  Zrr(u,v)  . 
Finally,  z(m,n),  the  output  at  the  second  stage,  can  be 
expressed  as 

z(m,n)  = zrr(m,n)  + j [zrr(m,n)  **hj  (m,n)  ] (2.45) 

The  spectrum  of  z(m,n),  denoted  as  Z(u,v),  is 

Z(u,v)  = Zrr(u,v)  + j [Zrr(u,v)Hj(u,v)  ] (2.46) 

where  Hj(u,v)  is  the  spectrum  of  the  impulse  response  hj(m,n) 
in  Equation  (2.41).  After  applying  the  spectrum  property  of 
Hj(u,v)  in  Equation  (2.39),  the  spectrum  Z(u,v)  in  Equation 
(2.46)  becomes 


Z(u,v) 


4Zr(u,v)  , 
2Zr(u,v)  , 


(0  < u < 7r  and  0 < v < n) 

(u=v=0)  or  (u=0  and  0 < v < n) 
or  (v=0  and  0 < u < tt) 

otherwise  (2.47) 


Because  zr(m,n)  is  assumed  to  be  real,  zrr(m,n)  is  also  real, 
and  consequently  z(m,n)  in  Equation  (2.45)  is  complex.  Since 
z(m,n)  is  complex  and  has  the  quadrant-plane  spectrum 
(Type  I) , by  definition,  it  is  a 2-D  analytic  signal — Type  I 
corresponding  to  a general  real  signal  zr(m,n) . Thus,  we  have 
Lemma  2.2  which  is  stated  as  follows. 

Lemma  2 . 2 In  a 2-D  signal  domain,  a 2-D  analytic  signal, 
whose  spectrum  is  the  quadrant-plane  (Type  I)  of  a 2-D  real 
discrete-time  signal,  can  be  generated  by  using  two  stages  of 
transformations.  (Refer  to  Figure  2.6) 
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2.3.5  The  generation  of  a 2-D  analytic  signal — Type  I using 
two  1-D  analytic  signals  (Case  lb) 

It  is  also  possible  to  obtain  the  desired  2-D  analytic 
signal — Type  I using  two  one  dimensional  analytic  signals. 
Let's  assume  a 2-D  complex  discrete-time  signal  z(m,n)  is 
separable,  that  is  z (m, n) =z (m) z (n) , where  z(m)  and  z(n)  are 
two  1-D  complex  signals.  We  want  to  show  that  if  both  z (m)  and 
z(n)  are  analytic,  then  z(m,n)  is  also  a 2-D  analytic  signal 
with  the  defined  quadrant-plane  spectrum  (Type  I)  . The  2-D 
complex  signal  z(m,n)  can  be  alternatively  expressed  as 

z(m,n)  = z(m)  z(n) 

= [ zr  (m)  + jzi  (m)  ] [zr(n)  + jz^n)] 

= [zr(m)zr(n)  - zi(m)zi(n) ] 

+j[zr(m)zi(n)  +zr(n)z1(m)]  (2.48) 

where  zr(  ) and  zi(-)  represent  the  real  and  imaginary  parts 
of  z(  ),  both  of  them  are  real  and  z^-)  is  the  discrete-time 
Hilbert  transform  of  zr(  ).  Since  both  zr(  ) and  zi(  ) are 
real,  the  2-D  signal  z(m,n)  in  Equation  (2.48)  is  complex. 
Furthermore,  by  using  the  spectrum  relations  between  the  real 
and  imaginary  part  of  a 1-D  analytic  signal  in  Equation 
(2.15),  it  is  straightforward  to  show  that  the  spectrum  of 
z(m,n)  in  Equation  (2.48),  denoted  as  Z(u,v),  is 

Z(u,v)  = Zr(u)Zr(v)  - [-jsgn(u)  Zr(u)  ] [-jsgn(v)  Zr(v)  ] 

+j  [-jsgn(v)  Zr(v)  Zr(u)  ] + j [-jsgn(u)  Zr(u)  Zr(v)  ] 


= Zr(u)Zr(v)[l  + sgn(u)sgn(v)  + sgn(u)  + sgn(v)] 
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that  is, 


Z(u,v) 


4Zr(u)Zr(v)  , 

for 

2Zr(u)Zr(v)  , 

for 

Zr(u)Zr(v)  , 

for 

0, 

(0  < u < n and  0 < v < n) 
(u=0  and  0 < v < n)  or 
(v=0  and  0 < u < n) 
u=v=0 

otherwise  (2.49) 


which  is  quadrant-plane  spectrum.  Since  the  2-D  discrete-time 
signal  z(m,n)  in  Equation  (2.48)  is  complex  and  has  quadrant- 
plane  spectrum  (Type  I) , by  definition,  it  is  a 2-D  analytic 
signal — Type  I. 

The  above  derivation  also  shows  that  a 2-D  separable 
analytic  signal — Type  I can  be  generated  by  using  1-D 
discrete-time  Hilbert  transform,  and  thus  leads  to  Lemma  2.3, 
the  main  result  in  this  section. 

Lemma  2 . 3 A complex  2-D  discrete-time  signal  formed  by 
the  product  of  two  1-D  analytic  signals  is  also  analytic,  and 
consequently  has  a periodic  quadrant-plane  spectrum  (Type  I) . 

2.3.6  The  generation  of  a 2-D  analytic  signal — Type  II  using 
2-D  discrete-time  Hilbert  transform  (Type  II) 

(Case  III 

Theoretically,  the  concept  and  procedures  of  deducing 
the  2-D  discrete-time  Hilbert  transform  (Type  II)  used  to 
generate  the  desired  analytic  signal  with  a half-plane 
spectrum  (Type  II)  are  similar  to  those  used  to  deduce  the 
discrete-time  Hilbert  tansform  (Type  I)  in  Section  2. 3. 4.1. 
Thus  referring  to  Figure  2.3(b),  by  definition  the  spectrum 
of  the  2-D  analytic  signal — Type  II,  denoted  as  Z(u,v),  is 
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Z(u,v)  = 0,  for  ( -7r  < u < 0)  (2.50) 

Let's  denote  the  2-D  discrete-time  signal  corresponding  to 
Z(u,v)  as  z(m,n)  which  must  be  complex,  due  to  the  asymmetric 
property  of  Z(u,v).  Thus,  z(m,n)  can  be  represented  by  its 
real  and  imaginary  parts,  that  is,  z(m,n)  = zr(m,n)+jzj(m,n)  , 
where  both  zr(m,n)  and  z^mjn)  are  real.  It  can  be  shown  that 

zr(m,n)  = (1/2) [z (m,n)+z* (m,n) ] 
jzi(m,n)  = (1/2) [z(m,n)-z*(m,n) ] 

where  z*(-,)  is  the  complex  conjugate  of  z (•,•)•  With  Zr(u,v) 
and  Z^UjV)  denoting  the  spectra  of  zr(m,n)  and  zi(m,n),  it  is 
easy  to  show  that 

Zr(u,v)  = ( 1/ 2 ) [ Z ( u , v ) + Z* (-u, -v) ] 
jZ1(u,v)  = (1/2) (Z(u,v)  - Z*(-u,-v) ] 

Using  the  periodic  half-plane  spectrum  property  of  Z(u,v), 
the  relations  among  Z(u,v),  Zr(u,v)  and  Zi(u,v)  are 

Z(u,v)  = 2Zr(u,v)  = 2jZi(u,v),  for  (0  < u < n ) 
Z*(-u,-v)=2Zr(u,v)=-2jZ1(u,v)  , for  (-7T  < u < 0)  (2.51) 

Thus,  from  Eguation  (2.51),  Zr(u,v)  and  Zf(u,v)  are  related 
as  follows 


zi(u,v)  = Hn  (u,  v)  Zr  (u,  v) 
Zr(u,v)  = -HII(u,v)Zi(u,v) 


(2.52) 
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where 


H„(u,v)  = -jsgnH(u,v) 


(2.53) 


and 


sgnH (u,v) 


(0  < u < ir)  and  (-jr  < v < n) 

(2.54) 

(-7r  < u < 0)  and  (-tt  < v < n) 


The  impulse  response  hH(m,n)  corresponding  to  Hn(u,v)  can  be 
obtained  by  carrying  out  the  2-D  discrete-time  inverse  Fourier 
transformation  of  Hn(u,v)  in  Equation  (2.53);  that  is, 


h„(m,n)  = { 


47T 


'0 

-7T 


7T 

0 


V 

. -7r 

n 

-7T 


j ej(um+vn)  dudv 


j e j (um+vn)  dudv} 


After  taking  the  integration  with  all  possible  combinations 
of  m and  n,  hn(m,n)  is  found  to  be 


hn(m,n)  = 


1-cos (mn) 

mn 
0 , 


n = 0 and  m t 0 
otherwise 


(2.55) 


Again,  the  2-D  impulse  response  hn(m,n)  in  Equation  (2.55) 
and  the  spectrum  of  a uniformly  truncated  version  of  hn(m,n) 
(the  lengths  of  both  m and  n are  11)  are  shown  in  Figure 
2.7(a)  and  (b) . 
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-3  X 
-4  • 

-5  X 

2 • • 

hII(m)=  , | 

mxr  -xn 

o: positive  value,  x: negative  value,  otherwise,  zero  value. 

(a) 


Figure  2.7  The  2-D  discrete-time  Hilbert  transform  (Type  II)  . 

(a) The  impulse  response  hn(m,n). 

(b) The  spectrum  of  uniformly  truncated  version 
of  (a) . (M=N=11) 
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After  deriving  the  2-D  impulse  response  hn(m,n)  of  the 
2-D  discrete-time  Hilbert  transform  (Type  II) , it  is  easy  to 
show  that,  for  an  arbitrary  real  2-D  discrete-time  signal 
zr(m,n),  the  2-D  complex  signal  z (m, n)  =zr (m, n)  +j zi  (m, n)  has  a 
half-plane  spectrum.  Thus,  by  definition,  z(m,n)  is  a 2-D 
analytic  signal — Type  II  which  can  be  generated  by  using  the 
2-D  discrete-time  Hilbert  transform  (Type  II) . 

Furthermore,  according  to  Equation  (2.52),  the  desired 
discrete-time  Hilbert  transform  relations  (Type  II)  between 
the  real  and  imaginary  part  of  a Type  II  analytic  signal  can 
be  expressed  as 

zi  (m,n)  = zr(m,n)  **hn  (m,n) 

zr(m,n)  = -zi  (m,n)  **hjj  (m,n)  (2.56) 

It  should  be  noted  that  based  on  the  property  of  the 
2-D  impulse  response  hn(m,n)  in  Equation  (2.55),  the  discrete- 
time Hilbert  transform  relations  (Type  II)  in  Equation  (2.56) 
can  be  simplified  as  a 1-D  discrete-time  convolution,  that  is 

zi(m,n)  = £ zr(k,n)hn  (m-k) 
k 

zr(m,n)  = - £ z1  (k,n)hn  (m-k)  (2.57) 

k 

2 . 4 Summary 

From  the  application  point  of  view,  the  2-D  discrete-time 
analytic  signals  and  their  associated  discrete-time  Hilbert 
transforms  have  been  discussed  and  derived  in  this  chapter. 
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The  derivation  procedures  can  be  extended  to  continuous  as 
well  as  multi-dimensioned  signals.  However,  in  order  to 
compatible  with  the  planar  array  antenna  system,  only  the  2- 
D case  is  emphasized.  It  was  found  that  it  is  possible  to 
obtain  the  quadrant-plane  spectrum  (Type  I)  of  a real  2-D 
discrete-time  signal  by  using  two  stages  of  transformations. 
If  the  original  2-D  signal  is  separable,  it  was  found  that  the 
2-D  analytic  signal — Type  I can  be  obtained  from  the  product 
of  the  associated  1-D  analytic  signal.  Furthermore,  a 
technique  was  developed  to  obtain  2-D  analytic  signal  with 
half-plane  spectrum  (Type  II)  . As  far  as  application  is 
concerned,  these  concepts  of  the  analytic  signals  and  the 
associated  half-sided,  quadrant-plane  and  half-plane  spectra 
will  provide  an  important  contribution  to  the  area  in  the 
design  of  the  radiation  patterns  of  an  array  antenna  system. 


CHAPTER  III 


ANALYSIS  OF  BROADSIDE  NARROW  BEAMWIDTH  ARRAY 
RADIATION  PATTERNS  USING  ANALYTIC  EXCITATIONS 

3 . 1 Introduction 

Having  defined  the  discrete-time  analytic  signals  and 
derived  the  reguired  mathematical  tool  - discrete-time  Hilbert 
transforms,  we  will  fully  analyze  the  utilization  of  the 
analytic  excitation1  for  an  array  antenna  system.  To  do  so,  we 
first  show  the  discrete-time  Fourier  transform  relations 
between  the  radiation  pattern  and  its  element's  excitation  of 
an  array  antenna  system.  It  is  this  relationship  that  shows 
why  the  concept  of  the  analytic  signal  theory  is  applicable 
to  the  design  of  the  array  antenna  radiation  pattern.  Then, 
the  performance  of  a narrow  beamwidth  radiation  pattern  using 
an  analytic  excitation  which  is  directly  corresponding  to  an 
even  symmetric  broadside  array  excitation  is  analyzed.  The 
analysis  procedures  will  be  divided  into  two  steps.  In  the 
first  step,  the  attractive  properties  of  the  array  antenna 
radiation  pattern  using  an  analytic  excitation  are 
illustrated.  These  attractive  properties  stimulate  the  idea 


!The  "analytic  signal"  mentioned  in  Chapter  II  will  be 
replaced  by  "analytic  excitation"  whenever  the  signal  is  the 
excitation  of  an  array  antenna  system. 


45 


46 


using  the  analytic  excitation  in  the  design  of  the  array 
antenna  radiation  patterns.  In  the  second  step,  the 
performance  of  radiation  patterns  using  a realizable  analytic 
excitation  is  quantitatively  analyzed.  From  there,  several 
properties  related  to  this  new  realization  technique  are 
extracted.  The  properties  are  useful  and  essential  for 
designing  several  array  antenna  radiation  patterns  that  will 
be  developed  in  the  next  chapters. 

3 . 2 Fundamentals  of  an  Array  Antenna  Theory 

This  section  begins  with  a description  of  a generalized 
far  field  radiation  pattern  of  an  Array  Antenna  System.  Two 
basic  but  important  array  structures,  linear  and  rectangular 
grid  planar  array  antennas,  will  be  studied  in  this 
dissertation.  Consequently,  the  discrete-time  Fourier 
transform  relations  between  the  far  field  radiation  pattern 
and  its  element's  excitation  will  be  shown. 

3.2.1  The  generalized  structure  of  an  array  antenna  system 

Generally,  the  phasor  for  a spherical  wave  that 
propagates  along  a radially  outward  direction  can  be 
expressed  as  [3] 

E(r  ,0,0)  = [A(9  ,<p)/r]e~i*r  (3.1) 

Where  A (6  ,<p)  is  the  complex  amplitude  of  the  propagated  wave, 
k=(2?r/A)  is  the  wave  number,  \ is  the  wavelength  of  the 
propagated  wave,  r is  the  distance  between  the  radiation 
source  and  the  observation  point  in  the  space,  and  (1/r)  is 
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the  expression  for  the  law  of  energy  conservation.  The  actual 
waveform  is  e(t )=Re  [E  (r , 6 , <p)  (wc*-)  ] r this  representation 
assumes  that  the  propagated  waveform  is  a temporal  narrow  band 
signal,  and  to  shorten  the  notation  the  temporal  variation  is 
omitted  in  the  wave  function  of  Equation  (3.1) . Two  propagated 
waves  with  same  temporal  frequency  may  be  summed  so  that  they 
combine  or  interfere  with  each  other,  and  it  is  this  principle 
that  forms  the  array  antenna  radiation  pattern. 

An  antenna  array  is  a collection  of  a specific  spatial 
distribution  of  several  identical  and  similarly  oriented 
radiation  elements,  and  can  be  considered  as  the  sampled 
version  of  an  aperture  antenna.  An  array  antenna  with  N 
arbitrary  elements  is  described  in  Figure  3.1,  where  one  of 
the  elements,  the  nth  element,  is  explicitly  shown. 


Figure  3.1  A generalized  geometrical  structure  of  an 
array  antenna  system. 
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In  Figure  3.1,  a phase  reference  is  assumed  to  be  the 
coordinate  origin.  The  position  vector  of  the  nth  array  element 
with  respective  to  the  phase  origin  is  rn 

rn  = XnXo  + YnYo  + ZnZo 

where  the  boldface  letter  represents  a vector,2  x0,  y0  and  zD 
are  unit  vectors  in  the  x,  y,  and  z direction  respectively. 
The  far  electromagnetic  field  that  launched  from  this  nth 
element,  according  to  Equation  (3.1),  is  proportional  to  [3] 

En  = ( [ J (n)  fn(0n,<*>n)]/Rn}<? "jkRn  (3.2) 

where  J(n)  and  fn(0n,0n)  are  the  complex  current  excitation 
and  the  far  field  pattern  of  the  nth  array  element 
respectively,  and  Rn  is  the  distance  between  the  nth  array 
element  and  the  observation  point  P.  Let's  assume  that  all 
individual  elements  have  the  same  radiation  pattern,  denoted 
as  f(8,<p),  and  they  are  distributed  within  a volume  with 
radius  a,  which  is  much  smaller  than  the  distance,  r,  from  the 
phase  reference  to  the  observation  point  (a«r) . This 
assumption,  by  using  the  Taylor  series  expansion,  assures 
a good  approximation  [3],  (l/Rn)«(l/r)  , 6n~6,  <Pn~<P,  and 


2 

Hereafter,  the  boldface  denotes  a vector.  Thus,  rn  is  a 
vector,  and  rn  is  a scalar  quantity. 

For  the  far  field  radiation  pattern,  this  assumption  is 
always  valid.  Generally,  the  distance,  r,  shall  be  equal  to 
or  greater  than  (2d2/A)  to  be  considered  as  the  far  field 
region,  where  d is  the  size  of  measurement  device  in  the  far 
field  observation  point.  In  the  far  field  region,  a spherical 
wave  can  be  locally  represented  as  a plane  wave,  and  thus 
simplify  many  propagated  wave  analysis  problems. 
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e-jkRp  w e-^kr  e jk(xnsin0cos0  + ynsin0sin0  + zncos 6) 

* ^-jkr  ejk(ru-rn) 

where  ru  is  the  normalized  unit  vector  along  the  r direction, 
which  is  defined  as 

r 

ru=  [ ] = sinflcos  </>xo  + sinflsin<£y0  + cos0zo 

M 

( ru * rn)  represents  the  inner  product  of  ru  and  rn,  and  k(ru*rn) 
is  the  measurement  of  the  phase  difference  between  the  phase 
reference  and  the  nth  element  as  measured  along  the  ru 
direction.  Then,  by  virtue  of  the  superposition  principle, 
the  total  far  field  strength  of  an  array  antenna,  denoted  as 
Es(r ,0,<p),  is  proportional  to  [3] 

N 

Es(r,0,0)  = (e'3Kr/r)f (0,0)  [ e J(n)c^(ru'rn)]  (3.3) 

n=l 

where  J(n)  is  the  normalized  complex  excitation  of  nth  array 
element,  and  the  component  in  the  bracket  of  the  right-hand 
side,  denoted  as  A{6  ,<p)  , 

A = [ S J(n)  ^k(rn*rn)  ] (3.4) 

n 

is  usually  called  the  array  factor  or  space  factor,  because 
it  depends  on  the  geometrical  structure  in  the  array  system, 
for  given  a specified  set  of  excitation  (J(n)}  and  the  number 
of  array  elements  N. 
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As  mentioned  earlier,  the  factor  (e-3kr/r)  in  Equation 

(3.3)  is  the  expression  for  the  law  of  energy  conservation  and 
the  phase  reference  and  is  not  a function  of  direction.  In 
antenna  theory,  the  far  field  radiation  pattern  F (0  ,<p)  of  an 
array  antenna  system,  which  is  defined  as  the  function  of 
direction  only,  is  identified  with  the  multiplicative  of  the 
element  factor  f (6  ,<p)  and  the  array  factor  A (0  ,<p)  in  Equation 

(3.4)  ;4  that  is, 

N 

F (0,0)  = f(0,0)[  Z J(n)^k(ru-rn)  ] (3.5) 

n=l 

The  objective  of  this  research  is  to  apply  a complex 
excitation  so  that  the  far  field  radiation  patterns  are 
obtained.  We  shall  not  deal  with  the  problems  related  to 
designing  the  element  factor,  f ($,</>),  for  the  individual 
radiation  elements.  Hence,  without  loss  of  any  generality,  an 
isotropical  element  factor  is  assumed  and  consequently,  the 
array  factor  A (0,0)  and  the  array  radiation  pattern  F(6 ,$)  are 
assumed  to  be  interchangable  throughout  this  dissertation. 

In  general,  a linear  phase  gradient  is  fed  along  with 
the  individual  array  element,  such  that  the  mainbeam  of  the 
narrow  beamwidth  radiation  pattern  can  be  steered  at  a 
specific  direction  (0o,<po)  . Mathematically,  the  normalized 
steering  vector  in  the  ( 90,<P0 ) direction,  denoted  as  rs,  can  be 
expressed  as 

4This  is  the  well-known  principle  of  the  pattern 
multiplication  of  an  array  antenna  system. 
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rs  = sin0ocos0oxo  + sin0osin0oyo  + cos0ozo 


(3.6) 


Then,  the  required  phase  gradient  on  each  array  element  can 
be  obtained  by  k times  the  inner  product  of  rs  and  rn,  and  the 
generalized  far  field  radiation  pattern  (strictly  speaking, 
the  array  factor)  of  an  N elements  array  antenna  becomes 


Two  specific  array  antenna  structures,  the  linear  and  the 
planar  rectangular  grid  array  systems  are  used  to  illustrate 
the  principle  of  the  general  array  antenna  theory.  For  some 
other  structures,  the  reader  can  refer  to  [3]. 

3.2.2  The  linear  array  antenna  system 

A typical  linear  array  antenna  which  consists  of  N 
linearly  oriented  array  elements  along  z axis  is  depicted  in 
Figure  3.2. 


(3.7) 


n=l 


P 


0 


z 


2 


3 


N-l 


Figure  3.2  The  linear  array  antenna  system. 
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In  Figure  3.2,  the  phase  reference  of  the  array  system 
is  assumed  to  be  at  the  origin  of  the  coordinates.  dn  is  the 
distance  between  the  nth  element  and  the  coordinate  origin. 
Let's  assume  that  the  interelement  distance  is  constant 
(assume  d) . Then,  the  location  of  the  nth  element  with  respect 
to  the  coordinate  origin  can  be  expressed  as 

rn  =ndZo 

and  the  phase  difference  between  the  nth  array  element  and  the 
phase  origin  is 

k(ru*rn)  = kndcosfl 

Further  assume  that  a constant  linear  phase  gradient 
k(rs-rn)  = kndcos#0 

is  applied  along  with  the  array  elements,  such  that  the 
mainbeam  maximum  of  the  radiation  pattern  is  oriented  at  the 
angular  direction  6 0.  Then,  based  on  Eguation  (3.7),  the  far 
field  radiation  pattern  of  the  linear  array  antenna  system  is 

N-l 

F (6)  = S j(n)e^n(kdcos^-kdcos^o)  (3.8) 

n=0 

It  is  normally  more  convenient  to  represent  the  linear  array 
radiation  pattern  given  in  Eguation  (3.8)  as  function  of  u, 
where 


u(0)  = kd  (cos0-cosflo) 
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Then,  the  linear  array  far  field  radiation  pattern  is 
alternatively  expressed  as 

N-l 

F(u)  = S J (n) e1un  (3.9) 

n=0 

which  shows  the  linear  array  far  field  radiation  pattern  and 
its  element's  current  excitation  forming  a Fourier  transform 
relationship  (with  reversed  sign  in  u)  .5  Refering  to  Figure 
3.2,  for  0o=(7T/2),  the  linear  array  is  called  broadside 
array . and  for  0o=O,  is  called  the  endfire  array.  For  a 
broadside  array  system,  the  angle  (6)  of  the  radiation  pattern 
from  0 to  7r  corresponding  to  the  range  of  u from  +kd  to  -kd, 
is  called  the  visible  range.  The  far  field  radiation  pattern 
of  a linear  array  system  is  inherently  a periodic  function  of 
u with  periodicity  2n . However,  it  is  different  from  the 


5It  is  without  loss  of  generality  that  the  index  n of  the 
summation  in  Equation  (3.9)  can  be  extended  from  -a>  to  +«  so 
as  to  meet  the  format  of  the  Fourier  transform  in  Equation 
(2.10),  although  the  number  of  array  elements  is  finite. 

It  is  also  possible  to  transform  an  N elements  linear 
array  system  intq>  a (N-l) -order  polynomial  function  of  z,  by 
letting  z = e Ju,  and  the  array  function  in  Equation  (3.9) 
becomes 


F(z)  = E J (n)  z-n 
n 

Since  u is  always  real,  the  magnitude  of  z is  always  equal  to 
unity,  and  the  array  pattern  thus  is  evaluated  along  with  a 
Schelkunoff  unit  circle  in  the  complex  plane.  For  an  N 
elements  array  system,  there  are  (N-l)  roots  in  polynomial 
F(z),  and  the  locations  of  these  roots  determine  the  shape  of 
the  array  radiation  pattern.  For  the  analysis  of  the  array 
radiation  pattern  using  z transform,  the  reader  can  refer  to 
[5-7]  . 
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spectrum  of  a discrete-time  signal  in  that  only  the  pattern 
corresponding  to  the  visible  range  need  to  be  noteworthy.  It 
is  easy  to  show  that  if  the  interelement  distance  d is  larger 
than  ( A/2 ) , then  the  visible  range  will  be  greater  than  2n . 
Thus,  for  a steering  angle  of  6=  0 (on  the  endfire 
direction) , a second  radiation  peak  (called  grating  lobe)  will 
appear  within  the  visible  range  of  the  radiation  pattern. 
Commonly,  the  interelement  distance  d of  a linear  array  system 
is  assumed  to  be  smaller  than  (A/2),  so  as  to  eliminate 
grating  lobe  within  the  visible  range.6 
3*2*3  The  planar  rectangular  grid  array  antenna  system 

A planar  rectangular  grid  array  antenna  system,  whose 
elements  are  distributed  in  x-y  plane  with  constant 
interelement  distances  dx  and  dy  is  depicted  in  Figure  3.3. 


Figure  3.3  The  rectangular  grid  planar  array  antenna 
system. 


g 

The  interelement  distance  will  be  assumed  to  be  half  of 
the  wavelength  (d=A/2)  throughout  this  dissertation. 
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In  Figure  3.3,  it  is  assumed  that  there  are  2M+1  and  2N+1 
elements  in  each  column  (x  axis)  and  row  (y  axis)  . The 
reference  element  is  exactly  positioned  at  the  coordinate 
origin,  the  interelement  distances  in  x and  y axes  are 
constant  and  denoted  as  dx  and  dy  respectively.  Thus,  the 
location  of  mth  element  in  x axis  is  xm=rndx,  and  the  nth  element 
in  y axis  is  yn=ndy,  in  which  -M  < m < M,  -N  < n < N.  Again, 
both  dx  and  dy  are  smaller  or  equal  to  (A/2),  such  that  the 
effect  of  grating  lobe  is  eliminated.  Then,  the  location  of 
the  mnth  element  of  the  rectangular  grid  planar  array  system 
can  be  expressed  as 

r (m,n)  = mdxx0  + ndyy0 

and  the  phase  difference  between  the  mnth  element  and  the  phase 
origin  is 

k[ruT(m,n)  ] = k[mdxsinflcos0  + ndysinflsin0] 

Again,  a constant  phase  gradient  is  required  to  apply  to  each 
planar  array  element,  such  that  the  radiation  pattern  maximum 
can  be  steered  to  a specific  direction  ( 0O,0O ) , and  the  phase 
gradiant,  according  to  Equation  (3.6),  should  be 

k[rs*r(m,n)  ] = k[mdxsin0ocos0o  + ndysintfosin0o] 

Then,  based  on  Equation  (3.7),  the  far  field  radiation  pattern 
of  the  rectangular  grid  planar  array  antenna  system  is 
proportional  to 
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F(0,0)=  £ £ J(m,n)  exp{  j [kmdx (sinflcos0  -sin0Qcos0o)  ] 
m n 

+j  [kndy  (sin0sin0-sin0osin0o)  ] } (3 . 10) 

where  J(m,n)  is  the  complex  excitation  of  the  mnth  array 
element.  From  analysis  point  of  view,  once  an  excitation 
function  is  given,  the  characteristics  of  the  field  pattern 
can  be  evaluated  by  using  Eguation  (3.10).  However,  from  the 
design  point  of  view,  it  is  more  convenient  to  have  the 
following  variable  transformations.  Let's  assume 

u(0,0)  = kdx(sin0cos0  - sinflocos0o) 

v(0,0)  = kdy(sin0sin0  - sin0osin0o)  (3.11) 

Then,  the  planar  far  field  radiation  pattern  in  Equation 
(3.10)  becomes 

F (u,  v)  = E £ J (m,  n)  (um+vn)  (3.12) 

m n 

which  again  exactly  shows  the  2-D  discrete-time  Fourier 
transform  relations  (with  reversed  sign  in  both  u and  v) 
between  the  far  field  radiation  pattern  of  the  rectangular 
grid  planar  array  system  and  its  element's  current  excitation. 
Due  to  the  periodic  property  of  an  exponential  function,  the 
planar  far  field  radiation  pattern  in  Equation  (3.12)  is  a 
periodic  function  on  both  u and  v with  the  period  of  2n . Once 
a desired  field  pattern  F (5,0)  is  specified,  it  can  be  mapped 
into  F(u,v)  via  Equation  (3.11).  Thus,  any  available 
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techniques  of  2-D  waveform  design  can  be  applied  to  the  design 
of  the  radiation  pattern  of  this  planar  array  system. 

It  is  always  true  that  the  array  excitation  J(m,n)  is 
controllable,  such  that  the  following  separable  condition  is 
realized  [8] 

J(m,n)  = J(m,o)J(o,n)  (3.13) 

Thus,  the  radiation  pattern  of  the  planar  array  antenna  as 
given  by  Equation  (3.12)  can  be  mathematically  expressed  as 
the  product  of  two  linear  array  radiation  patterns;  that  is, 

F(u,v)=  [ E J(m,o) e3um] [2  j(o,n)<?jvn] 
m n 

= F(u)F(v)  (3.14) 

The  separable  consequence  shown  in  Equation  (3.14)  provides 
a useful  idea  in  designing  the  radiation  pattern  of  a planar 
array  antenna  by  applying  the  well-developed  design  technique 
in  the  linear  array  system.  This  is  the  reason  why  most 
efforts  in  dealing  with  the  radiation  pattern  of  an  array 
antenna  have  been  concentrated  on  the  linear  array  system. 

Having  shown  the  Fourier  transform  relations  of  the  far 
field  radiation  pattern  and  its  element's  current  excitation 
of  an  array  antenna  system,  we  can  now  apply  any  concepts  in 
the  traditional  waveform  problem  to  the  design  of  radiation 
pattern  for  an  array  antenna  system.  Since  several  properties 
of  an  analytic  signal  are  so  attractive,  the  potential 


58 


advantages  of  using  such  an  analytic  excitation  for  an  array 

antenna  will  be  formally  exploited  and  analyzed  in  the 

following  sections.  It  is  noted  that,  in  order  to  be 

compatible  with  a traditional  signal  analysis  representation, 

hereafter,  it  is  desirable  to  change  the  variables  u,  v to 

-u,  -v  in  Equation  (3.9)  and  (3.12),  which  reverse  the 

radiation  pattern.  However  it  does  not  change  the  shape  and 

characteristics  of  the  radiation  pattern.  First,  we  discuss 

and  analyze  the  application  of  the  analytic  excitation  for 

designing  narrow  beamwidth  radiation  patterns.  The  procedures 

will  be  extended  to  designing  patterns  for  monopulse  tracking 

radar  systems  which  is  the  subject  of  Chapter  VI.  Due  to  the 

similarity  of  the  error  analysis  between  the  linear  and  the 

planar  array  systems,  the  derivation  and  analysis  will  be 

emphasized  on  the  linear  array  antenna  system. 

3 . 3 Analysis  of  Broadside  Narrow  Beamwidth  Array  Radiation 
Patterns  using  Analytic  Excitations 

The  frequent  need  to  both  detect  the  long  range  radar 
targets  and  discriminate  among  closedly-spaced  radar  targets 
has  led  to  the  high  beamwidth  resolution,  and  high  antenna 
directivity  radiation  pattern  design.  Both  requirements  are 
describable  in  terms  of  the  shape  of  the  radiation  pattern. 
A high  beamwidth  resolution  pattern  implies  the  pattern  has 
a narrow  mainlobe  beamwidth  and  low  sidelobe,  while  a high 
antenna  directivity  implies  a large  amount  of  energy  radiated 
toward  the  designated  direction  relative  to  the  total  energy 
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radiated.  It  has  been  recognized  that  a uniformly  excited 
array  antenna  system  which  has  half  wavelength  interelement 
distance  provides  the  maximum  directivity  and  the  radiation 
efficiency.  However,  the  high  sidelobe  ratio  of  the  uniform 
array  antenna  is  unacceptable  in  most  of  the  radar 
applications.  The  optimal  linear  in-phase  excited  array 
antenna  system  in  trade-off  between  the  mainbeam  width  and  the 
sidelobe  ratio  has  been  highlighted  by  the  work  of  Dolph  [19]  . 
Since  it  will  be  extensively  applied  throughout  this 
dissertation,  the  procedures  to  derive  the  Dolph-Chebyshev 
array  system  [19]  is  attached  in  Appendix. 

In  the  following  sections,  we  present  and  analyze  a new 
concept  to  improve  the  characteristics  of  a typical  symmetric 
broadside  narrow  beamwidth  radiation  pattern  of  an  array 
antenna  system,  such  as  directivity,  the  radiation  efficiency 
and  the  beamwidth  resolutions.  For  the  sake  of  convenience, 
these  characteristic  indices  of  a narrow  beamwidth  radiation 
pattern  to  be  studied  are  collectively  defined  as  follows. 

a.  The  directivity  of  an  array  antenna  is  a dimensionless 
quantity  and  defined  as  the  ratio  of  the  power  radiated  in 
the  direction  of  mainbeam  maximum  (0o,<po)  to  the  average  power 
density  from  the  antenna.  It  is  mathematically  expressed  as 


D 
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|f(*o,0o)  I' 
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b.  The  main-beam  radiation  efficiency  (77)  of  an  N 
elements  array  antenna  is  a dimensionless  quantity  and  defined 
as  the  normalized  ratio  of  the  power  radiated  in  the  direction 
of  mainbeam  maximum  (0O,0O)  to  the  total  power  of  the  current 
excitation  of  the  individual  array  element.  It  is 
mathematically  expressed  as  [16] 7 

|F(*O,0O)  I' 

n = ; (3.16) 

N[£(J(n)2] 

c.  Half  power  (or  3dB  power)  beam-width  (HPBW)  of  a 
narrow  beamwidth  radiation  pattern  of  a linear  array  system 
is  defined  as  the  angular  width  of  the  mainbeam  between  the 
angular  points  half  power  (3dB)  lower  than  the  mainbeam 
maximum  in  the  direction  (60)  ; that  is, 

HPBW  = ez  - (3.17) 

where  F(^2)=F(^1)=0.707F(5o) 

d.  First  null  beam-width  (BWnul1)  of  a narrow  beamwidth 


7It  has  been  shown  [16]  that  77=1  which  is  the  possible 
maximum  value  of  77  for  uniform  cophasal  excitation. 
Furthermore,  for  an  N elements  broadside  linear  array  antenna 
system  (60=n/2  or  u=0) , based  on  Equation  (3.9),  the  numerator 
in  Equation  (3.16)  is  proportional  to  the  algebraic  sum  of  all 
the  current  excitation,  thus  the  radiation  efficiency  (77)  in 
Equation  (3.16)  can  be  simplified  as: 

tgJ(n)  ]2 


N[gJ(n)2] 
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radiation  pattern  of  a linear  array  system  is  defined  as  the 
angular  width  between  the  first  zero  crossing  on  either  side 
of  the  mainbeam  maximum,  which  is  used  to  measure  the  mainbeam 
width  of  the  pattern?  that  is, 

BWnull  = *2  - *1  (3.18) 

where  F (02)  =F  (6  :)  =0  . 

e.  The  peak  sidelobe  (SLpeak)  is  the  ratio  of  the  maximum 
strength  of  the  radiation  pattern  outside  the  mainbeam  to  the 
strength  of  the  mainbeam  maximum. 

The  directivity  (D)  is  used  to  measure  the  transmitting 

power  gain  of  the  radiation  pattern  of  an  antenna  system.  The 

main  beam  radiation  efficiency  (n)  is  used  to  measure  the 

radiation  intensity  in  the  main  beam  direction  and  the 

excitation  efficiency  of  the  current  excitation  of  an  array 

antenna  system.  The  half-power  beamwidth  (HPBW) , first-null 

beamwidth  (BWnull)  and  the  peak  sidelobe  ratio  (SLpeak)  are  the 

pattern  characteristics  which  are  associated  with  the 

capability  of  detection  and  tracking  in  the  radar  application. 

3.3.1  The  potential  advantages  of  broadside  narrow  beamwidth 
array  radiation  patterns  using  analytic  excitations 

The  potential  advantages  of  the  broadside  narrow 
beamwidth  array  radiation  pattern  using  analytic  excitation 
are  first  illustrated  in  this  section.  The  performance  of  the 
broadside  narrow  beamwidth  array  radiation  pattern  excited  by 
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a realizable  analytic  excitation  will  be  analyzed  in  the 
following  Section  3.3.2. 

3. 3. 1.1  The  generation  of  a half-sided  narrow  beamwidth 
radiation  pattern  using  an  analytic  excitation 

Consider  we  are  given  a conventional  real  excitation 

Jr(n)  which  produces  a symmetric  broadside  narrow  beamwidth 

radiation  pattern  denoted  as  Fr(u) . A desired  analytic 

excitation  Jd(n)  can  be  generated  from  Jr(n)  by  applying  the 

discrete-time  Hilbert  transform  in  Equation  (2.19);  that  is, 

Jd(n)  = Jr  (n)  + jJi(n)  (3.19) 

where  J^n)  is  the  discrete-time  Hilbert  transform  of  Jr(n)  . 
After  using  the  spectrum  relations  of  the  real  and  imaginary 
part  of  an  analytic  signal  in  Equation  (2.15),  the  radiation 
pattern,  denoted  as  Fh(u) , corresponding  to  the  analytic 
excitation  in  Equation  (3.19)  is 

Fh(u)  = Fr(u)  + j [-jsgn(u)Fr(u)  ] 

= Fr(u)  [1  + sgn (u) J 


That  is, 


Mu) 


2Fr(u)  , 
Fr(u) , 


0 < U < 7T 

u=0  (3.20) 

-7T  < U < 0 


Equation  (3.20)  shows  that  the  radiation  pattern  excited  by 
an  analytic  excitation  is  half-sided.  As  a result,  both  the 
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HPBW  and  the  BWnull  beamwidths  of  Fh(u)  are  exactly  half  of 
those  in  Fr(u) . 

Using  the  same  procedures  for  given  a typical  2-D  real 
symmetric  excitation  Jr(m,n)  which  produces  a narrow  beamwidth 
planar  radiation  pattern,  denoted  as  Fr(u,v) , with  mainbeam 
maximum  at  u=v=0,  a half-plane  radiation  pattern  of  a planar 
array  system,  denoted  as  Fh(u,v)  (assumed  in  u plane),  can  be 
produced  by  using  a 2-D  analytic  excitation — Type  II  which  is 
denoted  as  Jd(m,n)  and  is  formed  as 

Jd(m,n)  = Jr(m,n)  + jj^n^n)  (3.21) 

where  J^it^n)  is  the  Hilbert  transform  of  Jr(m,n)  using  the 
2-D  discrete-time  Hilbert  transform  (Type  II) . It  is  easy  to 
show  that  the  half-plane  pattern  Fh(u,v)  is 


Fh(u,v) 


2Fr(u,v)  , for  (0  < u < 0)  and  (-7 r < v < n) 

Fr(u,v),  for  (u=0)  and  (-n  < v < n) 

0,  for  ( — 7T  < u < 0)  and  (-n  < v < ir) 

- (3.22) 


It  is  obvious  that  both  HPBW  and  BWnu11  beamwidths  of  the  half- 
plane radiation  pattern  Fh(u,v)  as  given  by  Equation  (3.22) 
are  half  of  those  in  Fr(u,v)  in  u plane.  Furthermore,  a 
desired  quadrant-plane  radiation  pattern  can  be  produced  from 
the  real  excitation  Jr(m,n)  by  using  the  procedures  as  given 
in  Section  2. 3. 3. 2.  (the  procedures  will  not  be  repeated 
here).  Thus,  according  to  Equation  (2.47),  the  quadrant-plane 
radiation  pattern  Fq(u,v)  is  expressed  as 
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Fq(u,v) 


4Fr(u,v),  for  (0  < u < n and  0 < v < it) 
2Fr(u,v),  for  (u=v=0)  or  (u=0  and  0 < v < 0) 
or  (v=0  and  0 < u < 0) 

0,  otherwise  (3.23) 


As  a result,  both  HPBW  and  BWnu11  of  Fq(u,v)  in  Equation  (3.23) 
are  half  of  those  in  Fr(u,v)  in  both  u and  v planes. 

The  common  features  of  these  asymmetric  radiation 
patterns  as  given  by  Equation  (3.20),  (3.22)  and  (3.23)  are 
that  the  beamwidths  (both  HPBW  and  BWnull)  are  reduced  by  one 
half,  and  the  radiation  strengths  over  the  radiation  region 
are  two  times  uniformly  increased  in  the  half-sided  pattern 
Fh(u)  and  the  half-plane  pattern  Fh(u,v)  , and  four  times 
uniformly  increased  in  the  quadrant-plane  pattern  Fq(u,v) . 
Therefore,  the  beamwidth  resolutions  are  doubly  increased,  and 
the  directivity  and  the  peak  sidelobe  remain  the  same  compared 
to  the  conventional  broadside  narrow  beamwidth  symmetric 
radiation  patterns.  Furthermore,  the  undesired  sidelobes  in 
either  half-side,  half-plane  or  three  of  four  quadrants  of 
Fh(u) , Fh(u,v)  and  Fq(u,v)  respectively  are  totally  eliminated. 
Thus,  we  have  the  following  first  assertion  related  to  the 
utilization  of  an  analytic  excitation  for  an  array  system. 

Assertion  3 . 1 An  analytic  excitation  can  theoretically 
produce  a desired  asymmetric  radiation  pattern  whose  beamwidth 
resolutions  are  doubly  increased  and  half  of  the  sidelobe 
region  is  totally  eliminated,  without  sacrificing  the 
directivity  and  peak  sidelobe  compared  to  a conventional 
symmetric  broadside  radiation  pattern  of  an  array  antenna 
system. 
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Two  linear  broadside  array  radiation  patterns  excited  by 
a uniform  excitation  and  the  corresponding  theoretical  ideal 
analytic  excitation  are  plotted  in  Figure  3.4(a)  and  (b) . 


Figure  3.4  The  broadside  narrow  beamwidth  radiation 

patterns  excited  by  a uniform  excitation  and  the 
corresponding  ideal  analytic  excitation. 

(a)  The  uniform  excitation. 

(b) The  corresponding  analytic  excitation. 

3. 3. 1.2  The  symmetry  property  of  an  analytic  excitation 

One  of  the  desirable  excitation  characteristics  for  a 
large  array  antenna  system  is  the  symmetric  property  of  the 
magnitude  of  the  excitation.  This  is  because  the  symmetric 
property  reduces  the  required  categories  of  power  divider  by 
one  half  in  the  physical  realization  of  the  beamforming 
network.  Therefore,  most  of  the  excitations  related  to  a large 
array  antenna  system  have  the  symmetric  property. 


(a) 


(b) 


Assume  that  Jr(n)  is  a typical  real  even  symmetric 
excitation  (i.e.  Jr(-n)=Jr(n) ) , for  an  array  antenna  system. 
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Then  based  on  the  properties  of  discrete-time  Fourier 
transform,  the  corresponding  radiation  pattern,  denoted  as 
Fr(u)  , should  be  real  and  even  symmetry.  That  is,  Fr (-u)  =Fr (u)  . 
The  radiation  pattern  F^u)  generated  by  the  Hilbert  transform 
of  Fr(u),  that  is  Fi  (u)  = [- jsgn (u)  Fr(u)  ] , should  be  pure 
imaginary  and  odd  symmetry.  Thus,  the  excitation  J^n) 
corresponding  to  F^u)  should  be  real  and  odd  symmetry.  The 
above  deductions  show  the  fact  that  the  Hilbert  transform  of 
a real  even  symmetric  excitation  is  always  a real  odd 
symmetric  excitation.  Thus,  if  an  analytic  excitation  J(n)  is 
formed  as 

J(n)  = Jr(n)  + j J1  (n) 

, J,(n) 

= a (n)  tan  x ( ) (3.24) 

Jr(n) 

2 2 

where  a (n)  = [ Jr  (n) +J,  (n)  ] 1/2  is  the  magnitude  and  tan-1()  is 
the  phase  of  J(n),  respectively.  Then,  based  on  the  even  and 
odd  symmetric  property  of  Jr(n)  and  J^n),  it  is  easy  to  show 
that  J(-n)  is  complex  conjugate  of  J(n).  Due  to  the  complex 
conjugate  symmetric  property,  the  desirable  even  symmetric 
property  in  the  magnitude  of  the  excitation  is  maintained  in 
the  analytic  excitation.  Thus,  the  symmetric  property  of  an 
analytic  excitation  is  stated  as  follows: 

Assertion  3.2  The  analytic  excitation  generated  from  a 
real  even  symmetric  excitation  has  the  property  of  complex 
conjugate  symmetry,  that  is,  its  magnitude  is  even  symmetry 
while  its  phase  is  odd. 
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3. 3. 1.3  The  improvement  of  radiation  efficiency  (n)  of  a 
tapered  excitation  using  an  analytic  excitation 

The  radiation  efficiency  (r?)  plays  an  important  role  in 
the  radiation  pattern  design  of  an  antenna  system.  This  is 
because  the  transmitting  efficiency  of  an  array  antenna  system 
is  related  to  this  parameter.  It  will  be  shown  that  an 
analytic  excitation  which  is  generated  by  a real  tapered 
excitation  always  has  better  radiation  efficiency  than  the 
original  real  tapered  excitation,  especially  compared  to  a 
heavily  tapered  real  excitation. 

Assume  that  an  analytic  excitation  J(n)  is  generated  from 
a real  excitation  Jr(n)  ; that  is, 

J(n)  = Jr(n)  + j J i (n)  (3.25) 

where  J^n)  is  the  discrete-time  Hilbert  transform  of  Jr(n)  . 
As  deduced  in  the  preceeding  section,  J^n)  should  be  an  odd 
symmetric  excitation.  Now,  let's  examine  what  will  later  be 
defined  as  the  complemental  property  of  the  real  and  imaginary 
parts  of  an  analytic  excitation  by  using  a pictorial 
illustration.  If  a given  (N=9)  conventional  normalized  array 
excitation  Jr(n)  is  as  depicted  in  Figure  3.5(a),  then,  based 
on  the  operation  process  of  the  discrete-time  convolution  and 
the  sequence  of  the  impulse  response  of  the  discrete-time 
Hilbert  transform  h(n)  in  Equation  (2.18),  it  is  not  hard  to 
show  that  Ji  (n)  will  be  as  in  Figure  3.5(b). 
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Jr(n) 


(a) 

J i (n) 


Figure  3.5  The  pictorial  illustration  of  the  complement a 1 
property  of  the  real  and  the  imaginary  parts  of 
a typical  analytic  excitation. 

(a)Jr(n).  (b) Ji (n)  . 

The  meaning  of  the  pictorial  results  shown  in  Figure  3 . 5 
can  be  stated  as  follows:  whenever  Jr(n)  approaches  maximum 
value,  J^n)  always  approaches  a zero  value,  and  at  the  two 
tapered  terminals  of  Jr(n),  J^n)  normally  becomes  dilated. 
This  is  the  complemental  property  between  the  real  and 
imaginary  parts  of  an  analytic  excitation.  An  extremely 
evident  example  of  the  complemental  property  of  an  analytic 
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excitation  is  the  sinusoidal  excitation,  that  is,  an  analytic 
excitation  is  formed  as  a sum  of  a sine  and  cosine  excitation, 
that  is,  J (n)  =cos (n) +j sin (n)  (it  is  well  known  that  the  cos(n) 
and  sin(n)  form  a Hilbert  transform  pair) , and  it  is  easy  to 
show  that  the  magnitude  of  J(n)  is  always  unity  for  all  the 
array  elements.  Furthermore,  the  following  inequality  holds 

I J(n)  I = |Jr(n)  + j JJn)!  > |jr(n)|  (3.26) 

The  complemental  property  of  the  real  and  imaginary  parts 
of  an  analytic  excitation  illustrated  in  Figure  3.5  and  the 
inequality  equation  of  (3.26)  assure  the  improvement  of 
tapering  condition  as  well  as  radiation  efficiency  (rj)  of  the 
analytic  excitation  over  the  real  symmetric  excitation  (will 
be  seen  from  computer  simulations  later) . Therefore,  the 
assertion  3.3,  the  result  of  this  section,  can  be  stated  as 
follows: 

Assertion  3.3  Due  to  the  complemental  property  of  the 
real  and  the  imaginary  parts,  an  analytic  excitation  always 
provides  higher  radiation  efficiency  than  a tapered  real 
excitation  for  an  array  antenna  system. 

3.3.2  The  analysis  of  broadside  narrow  beamwidth  array 

radiation  patterns  using  realizable  analytic  excitations 

In  the  previous  Section  3.3.1,  we  have  illustrated  the 
outstanding  and  special  characteristics  of  the  broadside 
narrow  beamwidth  array  radiation  pattern  using  a theoretical 
analytic  excitation.  However,  such  a theoretical  analytic 
excitation  is  unrealizable,  because  it  can  be  achieved  only 
with  an  infinite  number  of  array  elements.  A physically 
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realizable  array  antenna  system  consists  of  a finite  number 
of  array  elements.  Thus,  the  realization  of  an  analytic 
excitation  in  the  array  antenna  system  becomes  an 
approximation  problem. 

Consider  an  N=2M+1  elements  linear  broadside  array 
antenna  system,  where  M is  an  integer,  and  a typical 
conventional  real  excitation  Jr(n),  where  -M  < n < M. 
Accordingly,  the  desired  analytic  excitation  Jd(n)  can  be 
generated  by 

Jd(n)  = Jr(n)  + j[Jr(n)*h(n)  ] (3.27) 

where  * denotes  the  discrete-time  convolution,  and  h(n)  is  the 
impulse  response  of  1-D  discrete-time  Hilbert  transform 
defined  in  Equation  (2.18).  Then,  the  desired  half-sided 
radiation  pattern  Fd(u)  corresponding  to  Jd(n)  is 

Fd(u)  = Mu)  + j <Fr(u)  [-jsgn(u)  ] } 

= Fr(u)  + [Fr (u)  sgn (u)  ] (3.28) 

where  Fr(u)  is  the  radiation  pattern  corresponding  to  Jr(n). 
It  is  obvious  that  the  analysis  of  the  half-sided  radiation 
pattern  Fd(u)  using  a realizable  analytic  excitation  is 
equivalent  to  the  analysis  of  the  approximation  to  the  ideal 
transformed  pattern  [ -j  Fr  (u) sgn (u) ] . Based  on  the  property  of 
the  discrete-time  Hilbert  transformation,  the  excitation  of 
[Jr(n)*h(n)]  corresponding  to  the  pattern  of  [ -j Fr (u) sgn (u) ] 
is  an  infinite  sequence,  due  to  the  infinite  length  of  h(n). 
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For  a finite  elements  array  antenna  system,  the  most 
straightforward  approach  to  obtain  a realizable  analytic 
excitation  is  to  truncate  the  infinite  sequence  of 
[Jr(n)*h(n)]  into  a length  of  N with  some  suitable  window 
function;8  that  is,  the  truncated  analytic  excitation  Jt(n) 
can  be  represented  as 

Jt(n)  = Jr(n)  + j [Jr(n)  *h(n)  ]w(n)  (3.29) 

where  w(n)  represents  a typical  N points  truncation  window, 

w(n)  =0  M < | n | 

By  applying  the  window  theorem  [17,pp.59],  the  realizable 
radiation  pattern  Ft(u)  corresponding  to  Jt(n)  in  Equation 
(3.29),  becomes 

Ft(u)  = Fr(u)  + j{Fr(u)  [-jsgn(u)  ] }*W(u) 

= Fr(u)  + [Fr(u)sgn(u)  ]*W(u)  (3.30) 

where  W(u)  is  the  spectrum  of  the  selected  window  function 
w(n),  and  { [Fr (u) sgn (u) ] *W(u) } represents  the  periodic 
convolution  of  the  ideal  transformed  pattern  [Fr (u) sgn (u) ] 
and  the  spectrum,  W(u),  of  the  window  function.  With  F^u) 
denoting  the  resulting  spectrum  after  periodic  convolution, 
the  periodic  convolution  is  defined  as  [56] 

8 

In  this  dissertation,  our  principal  interest  is  to 
convey  a new  concept  to  the  design  of  the  radiation  pattern 
of  an  array  antenna  system,  and  the  window  method  will  be 
adopted  in  the  analysis  and  design  of  the  radiation  pattern 
using  analytic  excitation. 
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Fi(u) 
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[Fr(0sgn(0  ]W(u-£)  d£ 

2n  J -7T 

V 


(3.31) 


Actually,  the  realization  of  the  ideal  pattern 
[Fr (u) sgn (u) ] which  has  discontinuity  in  the  u=0  direction  by 
using  windowing  method  is  similar  to  the  realization  of  an 
ideal  finite  impulse  response  (FIR)  discrete  filter  by  using 
windowing  method  [56-57],  a subject  which  has  received  a great 
deal  of  study,  and  a particular  set  of  tapering  windows  which 
moderated  the  Gibbs  phenomenon  [56]  at  the  discontinuity  has 
been  found  [12],  Thus,  it  is  well  known  that  due  to  the 
effects  of  the  window  function,  the  realized  pattern  Fi  (u)  in 
Equation  (3.31)  becomes  a smeared  version  of  [ Fr (u) sgn (u) ] . 
Using  a typical  even  symmetric  pattern  Fr(u)  and  a typical 
spectrum  W(u)  of  a window  function,  the  periodic  convolution 
in  Equation  (3.30)  is  illustrated  in  Figure  3.6.  Figure  3.6(a) 
is  the  typical  symmetric  broadside  radiation  pattern  Fr(u) , 
Figure  3.6(b)  shows  the  convolution  process  at  u=0,  and  the 
pattern  in  Figure  3.6(c)  represents  the  relizable  transformed 
pattern  F^u)  which  is  the  smeared  version  of  the  desired 
pattern  [Fr(u)  sgn(u)  ] . As  seen  from  Figure  3.6,  the  most 
familiar  effects  of  the  window  method  are  that  the 
discontinuity  of  the  pattern  [Fr (u) sgn (u) ] at  u=0  becomes  a 
roll-off  region,  and  the  pattern  of  [Fr(u)(sgn  u)  ] is  also 
leaked  by  the  spectrum,  W(u),  of  the  window  function  for  all 
range  of  u.  However,  a quantitative  estimation  of  the  width 
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Figure  3.6  The  illustration  of  periodic  convolution. 

(a) The  symmetric  broadside  pattern  Fr(u). 

(b) The  convolution  process  at  u=0. 

(c) The  realized  transformed  pattern  F^u). 

of  the  roll-off  region  and  the  extent  of  the  leakage  effect 
of  the  windowing  method  is  unobtainable.  This  is  because  not 
only  they  depend  on  both  the  overall  shapes  of  the  original 
symmetric  pattern  Fr(u)  and  the  spectrum  W(u)  of  the 
truncation  window,  but  also  the  symmetric  pattern  Fr(u)  and 
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the  spectrum  W(u)  always  have  approximately  the  same  mainbeam 
width  and  mainlobe  width  under  this  realization  situation. 
Nevertheless,  it  will  be  shown  by  computer  simulations  that 
the  performance  of  the  realized  asymmetric  radiation  pattern 
using  a realizable  analytic  excitation  is  found  to  be 
proportional  to  the  parameter  Main-Lobe  Width  Ratio  (MLWR) , 
which  is  defined  as 


MLWR 


BW, 


wnull 


BW; 


inul  1 


2BW, 


wnull 


BW 


rnul  1 


(3.32) 


where  BWwnull  and  BWrnull  are  the  mainlobe  width  of  the  spectrum 
W(u)  of  the  window  function  and  the  mainbeam  width  of  the 
original  even  symmetric  pattern  Fr(u)  respectively,  and  BW1null 
is  the  half-side  mainbeam  width  of  the  pattern  [ Fr (u) sgn (u) ] 
which  is  half  of  BWrnull.  For  the  definitions  of  BWwnul1,  BWrnull 
and  BWinul1,  please  refer  to  Figure  3.6. 

It  is  clear  that  the  performance  of  the  realizable 
asymmetric  radiation  pattern  Ft(u)  in  Equation  (3.30)  is 
determined  by  the  extent  of  the  pattern  Fi  (u)  approximation  to 
the  ideal  pattern  [Fr(u)  sgn(u)  ] . As  far  as  the  distortion 
which  results  from  the  truncation  window  is  concerned,  it 
should  be  equivalent  to  the  extent  of  the  distortion  of  the 
ideal  pattern  [Fr  (u)  sgn (u)  ] around  the  u=0  results  from  a 
realizable  Hilbert  transformer  that  is  also  assumed  to  be 
implemented  by  the  windowing  method.  Consequently,  the 
selection  of  the  truncation  window  becomes  the  narrower  the 
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mainlobe  width  the  better.  That  is,  the  spectrum  of  the  window 
function  should  possess  the  narrowest  mainlobe  width  under  a 
specified  sidelobe  ratio,  or  vice  versa.  It  is  well  known  that 
the  optimal  window  function  which  possesses  the  required 
properties  is  the  Dolph-Chebyshev  window  [12].  Actually,  the 
Dolph-Chebyshev  window  is  exactly  the  same  as  the  Dolph- 
Chebyshev  linear  array  system  where  the  interelement  distance 
is  set  to  half  of  the  wavelength.9  For  the  derivation  of  the 
optimal  Dolph-Chebyshev  window,  refer  to  Appendix.  The 
property  of  the  Dolph-Chebyshev  window  function  is  that 
whenever  the  sidelobe  ratio  R is  specified,  the  mainlobe  width 
is  automatically  minimized. 

After  acquiring  the  optimal  truncation  window,  the 
performance  of  an  array  antenna  radiation  pattern  using 
realizable  analytic  excitation  can  be  thoroughly  analyzed. 
However,  a generalized  error  analysis  is  unobtainable.  By 
using  computer  simulation,  two  specific  linear  array  antenna 
systems  are  used  to  demonstrate  the  achievements  of  this  new 
realization  technique.  The  procedures  used  to  analyze  and  to 
simulate  the  realizable  analytic  excitation  and  its 
corresponding  asymmetric  radiation  pattern  are  described  in 
Figure  3.7. 

9Refer  to  Equation  (3.9),  when  d=(A/2),  the  visible  range 
of  the  linear  array  system  is  from  -n  to  n (total  2n  range), 
which  is  exactly  the  same  as  the  periodicity  of  the  spectrum 
of  the  discrete-time  signal.  Thus  the  Dolph-Chebyshev  linear 
array  radiation  pattern  is  equivalent  to  the  spectrum  of  a 
discrete-time  window  function  which  is  optimal  in  the 
narrowest  mainlobe  width. 
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X : mult ipl icat ion . 
E: summation. 


Figure  3.7  The  generations  of  the  realizable  analytic 
excitation  and  the  corresponding  asymmetric 
radiation  pattern. 

From  Figure  3.7,  once  a real  excitation  Jr(n)  is 
specified,  the  imaginary  part  of  a realizable  analytic 
excitation  can  be  obtained  from  Jr(n)  by  carrying  out  the 
discrete-time  Hilbert  transform,  and  multiplying  it  with  a 
selected  truncation  window  w(n) . The  analytic  excitation  Jt(n) 
is  then  the  quadrature  summation  of  the  real  excitation  Jr(n) 
and  the  truncated  version  of  Hilbert  transform  of  Jr(n). 
Accordingly,  the  corresponding  asymmetric  array  radiation 
pattern  Ft(u)  is  the  Fourier  transform  of  the  finite  length 
analytic  excitation  Jt(n).  In  doing  the  computer  simulation, 
there  are  two  steps  that  need  to  be  emphasized.  The  first  one 
is  that  based  on  the  discrete-time  convolution  operation,  the 
length  of  the  impulse  response  h(n)  of  the  discrete-time 
Hilbert  transfrom  which  is  required  to  compute  the  N points 
imaginary  part  of  an  analytic  excitation  is  reduced  to  the 
points  from  -(N-l)  to  (N-l) , for  total  (2N-1)  points.  The 
second  one  is  that  the  Romberg  integration  [62-63],  which 
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applies  the  Richardson  extrapolation  process  to  improve  the 
accuracy  of  the  trapezoidal  integration  rule,  is  used  to 
compute  the  average  radiation  power  in  the  denominator  of 
Equation  (3.15)  that  is  required  to  evaluate  the  antenna 
directivity  given  a particular  set  of  excitation. 

3.3.3  Analysis  examples 

Example  1 In  this  first  analysis  example,  we  use  the 
heavily  tapered  Hamming  function  as  an  original  real  even 
symmetric  excitation  to  generate  the  desired  analytic 
excitation.  Then,  the  performance  of  the  radiation  patterns 
excited  by  the  Hamming  and  the  corresponding  realizable 
analytic  excitations  are  analyzed  and  compared. 

The  Hamming  excitation  of  an  N (N=2M+1,  M is  an  integer) 
elements  array  system  can  be  expressed  as  [12] 


Jh(n) 


— 2n7r 

0.54  + 0.46  cos  ( ), 

N-l 

0, 


-M  < n < M, 
otherwise 


(3.33) 


and  its  radiation  pattern,  denoted  as  Fh,m(u)  , is 


sin[ (-f-)u] 


Fham(U)=  0.54  (- 


-)  + 0 . 23  (- 
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sin[(^-)  (u+f^y)] 
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Equation  (3.34)  shows  that  the  pattern  excited  by  the  Hamming 
function  is  the  summation  of  three  weighted  and  shifted 
[sin(Nx)/sin(x)  ] waveforms  that  effectively  eliminate  the  peak 
sidelobe  of  the  [sin(Nx)/sin(x) ] waveform,  thus  its  normalized 
peak  sidelobe  can  be  as  low  as  -41  dB.  However,  the  mainbeam 
width  (BWnull)  is  wider  than  that  of  uniformly  excited  array 
system  ( about  lOw/N)  [12].  Due  to  being  heavily  tapered,  the 
radiation  efficiency  of  the  Hamming  excitation  is  poor.  The 
truncated  (or  realizable)  analytic  excitation  related  to  the 
Hamming  excitation  can  be  obtained  by  implementing  the 
transformation  and  truncation  processes  in  Figure  3.6.  Thus, 
the  truncated  analytic  excitation,  denoted  as  Jth(n)  and  the 
corresponding  asymmetric  radiation  pattern,  denoted  as  Fth(u), 
are  expressed  as 

Jth (n)  = Jh(n)  + j [Jh(n)  *h(n)  ]w(n) 

Fth(u)  = Fham(u)  + [Fham(u)sgn(u)  ]*W(U)  (3.35) 

where  W(u)  is  the  spectrum  of  the  Dolph-Chebyshev  truncation 
window  mentioned  earlier.  The  various  radiation  patterns  used 
to  analyze  and  to  compare  the  performance  of  the  Hamming 
excitation  and  the  corresponding  realizable  analytic 
excitation  are  plotted  in  Figure  3.8(a)  through  Figure  3.8(d), 
where  the  number  of  array  elements  is  assumed  to  be  21  (N=21)  . 
Figure  3.8(a)  overlaps  the  even  symmetric  radiation  pattern 
Fham(u)  excited  by  the  Hamming  excitation  and  the  transformed 
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pattern  { [Fham(u)  (sgn  u)  ]*W(u)  } corresponding  to  the  excitation 
{ [ Jh (n) *h (n) ] w(n) > , where  the  truncation  window  is  the  Doloh- 
Chebyshev  window  with  sidelobe  ratio  R=30dB.  The  Hamming 
pattern  Fham(u)  and  the  realized  asymmetric  radiation  pattern, 
let  it  be  denoted  as  Fth(u)  R=3  0,  corresponding  to  the 
summation  of  the  two  patterns  in  Figure  3.8(a)  are  plotted 
together  in  Figure  3.8(b).  Both  Figure  3.8(c)  and  Figure 
3.8(d)  have  the  same  interpretations  as  Figure  3.8(a)  and 
Figure  3.8(b)  respectively.  However,  the  truncation  window  is 
the  Dolph-Chebyshev  window  with  sidelobe  ratio  R=20dB.  thus, 
let  the  realized  asymmetric  pattern  in  Figure  3.8(d)  be 
denoted  as  Fth(u)  R=20.  Let's  first  examine  the  shapes  of  the 
realized  asymmetric  radiation  pattern  Fth(u)  R=30  in  Figure 
3.8(b)  . The  most  noticeable  phenomenon  is  a high  peak  sidelobe 
that  appears  within  the  sidelobe  region  of  u < 0.  The  cause 
°f  this  high  peak  sidelobe  can  be  examined  from  Figure  3.8(a)  , 
which  shows  the  peak  sidelobe  is  the  result  of  the  mis-match 
of  the  mainbeams  between  the  Hamming  symmetric  pattern  and  its 
Hilbert  transformed  pattern.  It  is  the  consequence  of  the 
periodic  convolution  between  the  idealized  Hilbert  transformed 
pattern  and  the  spectrum  of  the  truncation  window,  that 
induces  the  mainbeam  mis— match  phenomenon.  Next,  by  examining 
the  realized  asymmetric  pattern  Fth(u)  R=20  in  Figure  3.8(d), 
we  find  that  the  mainbeam  mis-match  phenomenon  is  moderated. 
This  is  because  the  Dolph-Chebyshev  truncation  window  used  to 
realize  the  asymmetric  pattern  Fth(u)  R=20  has  a higher 
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Figure  3.8  The  Hamming  and  the  corresponding  realized 

asymmetric  patterns. 

(a)  The  Hamming  and  its  Hilbert  transformed 
patterns. (w(n)=Dolph-Chebyshev  with  R=30  dB) 

(b) The  Hamming  and  the  realized  asymmetric 
pattern  Fth(u)  R=30,  corresponding  to  (a)  . 

(c) The  Hamming  and  its  Hilbert  transformed 
patterns. (w(n)=Dolph-Chebyshev  with  R=20  dB) 

(d) The  Hamming  and  the  realized  asymmetric 
pattern  Fth(u)  R=20,  corresponding  to  (c)  . 
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(d) 


Figure  3 . 8 


continued 
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sidelobe  ratio  and  consequently  narrower  mainlobe  width 
compared  to  that  of  the  window  used  to  generate  the  pattern 
Fth(u)  R=30.  This  is  an  increase  in  the  parameter  the  Main- 
Lobe  Width  Ratio,  (MLWR) . The  improvement  of  the  mainbeam  mis- 
match phenomenon  can  be  seen  from  Figure  3.8(c).  As  a result, 
the  peak  sidelobe  produced  by  the  mainbeam  mis-match  is 
dwindled  in  the  asymmetric  pattern  Fth(u)  R=20  of  Figure 
3.8(d).  However,  theaccomplishment  of  the  dwindled  peak 
sidelobe  is  achieved  by  sacrificing  the  overall  peak  sidelobe 
of  the  realized  asymmetric  pattern.  The  quantitative 
characteristics  of  the  patterns  of  Fham(u)  , Fth(u)  R=30  and 
Fth(u)  R=20  are  listed  in  Table  3.1.  In  the  characteristic  of 
peak  sidelobe,  the  overall  peak  sidelobe  scenario  SLpeak  is 
further  decomposed  into  left  SLLpeak  and  right  SLRpeak  sidelobe 
regions  with  respect  to  the  mainbeam  of  the  asymmetric 
pattern,  in  order  to  clearly  describe  the  inherent  property 
of  the  asymmetric  radiation  pattern. 


Table  3 . 1 The  characteristics  of  the  broadside  narrow 

beamwidth  radiation  patterns  excited  by  Hamming 
and  the  two  corresponding  analytic  excitations. 


type  of 
pattern 

HPBW 

(Rad) 

BWnull 

(Rad) 

MLWR  r? 

(None) (None) 

D SLL  k SLR  k 

(None)  (dB)  (dB) 

Fham(U) 

.401 

1.359 

— 

.709 

14.885  -40.64  -40.64 

Fth(u)  R=3 0 

.361 

1.268 

1.55 

.789 

16.575  -21.64  -34.63 

Fth(u)  R=20 

.322 

.907 

2.01 

.877 

18.419  -24.85  -24.79 

By  comparing  the  quantitative  analysis  results  of  these 
three  radiation  patterns  tabulated  in  Table  3.1,  we  find  that 
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the  patterns  excited  by  the  analytic  excitations  have  the 
advantages  stated  in  Section  3.3.1,  that  is,  they  have  higher 
beamwidth  resolutions  in  both  HPBW  and  BWnull,  and  have  higher 
antenna  directivity  and  radiation  efficiency.  However,  the 
higher  peak  sidelobes  that  appear  in  the  realized  asymmetric 
patterns,  especially  the  undesired  peak  sidelobe  produced  by 
the  convolution  process,  has  greatly  demolished  these 
advantages  gained  from  the  analytic  excitation.  By  examining 
the  tabulated  characterictics  of  the  two  realized  asymmetric 
radiation  patterns  Fth(u)  R=30  and  Fth(u)  R=20,  we  find  that 
the  shapes  of  both  realized  asymmetric  radiation  patterns  are, 
as  expected,  very  much  determined  by  the  spectrum 
characteristics  of  the  Dolph-Chebyshev  truncation  window.  The 
sidelobes  of  both  realized  asymmetric  radiation  patterns  are 
leaked  by  the  sidelobe  of  the  truncation  window  which  is 
normally  higher  than  that  of  the  original  symmetric  pattern. 
It  also  indicates  that  the  higher  the  Main-Lobe  Width  Ratio 
(MLWR) , the  more  advantages  obtained  from  the  analytic 
excitation,  which  include  the  HPBW,  BWnull,  directivity  (D)  , 
radiation  efficiency  (rj)  , and  the  lower  peak  sidelobe  which 
is  produced  by  the  periodic  convolution  process.  However, 
those  advantages  are  limited  by  the  required  sidelobe  ratio 
of  the  realized  asymmetric  radiation  pattern. 

The  excitations  for  the  symmetric  radiation  patterns 
Fham(u)  and  Fth(u)  R=2°  are  listed  in  Table  3.2,  where  only  half 
of  the  element's  excitations  are  tabulated,  and  the  another 
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half  of  the  element's  excitations  are  just  the  complex 
conjugate  of  the  tabulated  results,  that  is,  J(-n)=J*(n). 
Comparing  the  two  excitation  sets  listed  in  Table  3.2  show 
that  the  heavily  tapered  amplitudes  of  the  original  Hamming 
excitation  are  greatly  improved  in  the  corresponding  analytic 
excitation,  which  is  another  important  advantage  of  the 
analytic  excitation. 


Table  3.2  The  normalized  Hamming  and  the  corresponding 
analytic  excitations. 

(N  =21,  J(-n)=J*(n) ) 


n 

J(n) 

Fham(U) 

Fth(u) 

R=20 

(amplitude) 

(amplitude) 

(phase) 

10 

. 080 

.504 

80.87 

9 

.103 

.254 

66.17 

8 

. 168 

.361 

62 .26 

7 

.270 

.480 

55.84 

6 

.398 

.614 

49.64 

5 

.540 

.737 

42.85 

4 

.682 

.842 

35.90 

3 

.810 

.918 

28.02 

2 

.912 

.967 

19.43 

1 

.978 

.992 

9.94 

0 

1.000 

1.000 

0.00 

Example  2 The  purpose  of  this  example  is  to  analyze  the 
performance  of  a Dolph-Chebyshev  excited  linear  array  system 
and  its  corresponding  analytic  excitations  (instead  of  Hamming 
excitation  as  used  in  Example  1)  , so  that  a deeper  insight 
into  the  pattern  characteristics  of  an  analytic  excitation  can 


be  obtained. 
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The  radiation  pattern  of  an  N elements  Dolph-Chebyshev 
linear  array  system  can  be  expressed  as  (refer  to  Appendix) 

Fc(u)  = tn-i(xocos(u/2)  ) (3.36) 

where  Tn_j(x)  is  the  (N-l)th  Chebyshev  polynomial,  x0  is  chosen 
such  that  Tn_1(x0)=R,  and  R is  the  required  peak  sidelobe 
ratio.  Let's  assume  that  the  sidelobe  ratio  R=30dB  Dolph- 
Chebyshev  array  system  is  the  original  real  excitation,  and 
the  two  truncation  windows  used  to  realize  the  analytic 
excitations  are  Dolph-Chebvshev  windows  with  R=25  dB  and 
R=18dB.  The  number  of  array  elements  is  the  same  as  that  in 
Example  1 (N=21) . The  symmetric  radiation  pattern  Fc(u)  of  a 
Dolph-Chebyshev  array  system  and  the  corresponding  two 
realized  asymmetric  radiation  patterns,  denoted  as  Ftc(u)  R=25 
and  Ftc(u)  R=18 , are  plotted  together  in  Figure  3.9(a)  and 
Figure  3.9(b)  respectively,  and  their  quantitative  pattern 
characteristics  are  listed  in  the  first  three  rows  of  Table 
3.3,  where  they  are  indicated  as  Fc(u)(30),  Ftc(u)  R=25  and 
Ftc(u)  R=18.  By  examining  and  comparing  the  shapes  of  the 
radiation  patterns  in  Figure  3.9,  and  the  quantitative 
characteristics  of  these  radiation  patterns  in  Table  3.3,  we 
can  conclude  that  the  analysis  results  related  to  the  pattern 
characteristics  of  an  analytic  excitation  are  very  much 
similar  to  those  of  Example  1.  That  is,  the  advantages  gained 
from  the  analytic  excitation  include  HPBW,  BWnu11,  antenna 
directivity  and  radiation  efficiency,  and  the  disadvantage  of 
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Chebyshev  (ROOdS) 


u Chadian) 


asymmetric  pattern 


(a) 


□ Chebyshev  CR=30d8^ 


u Chadian) 


+ asymmetric  pattern 


(b) 

Figure  3.9  The  Chebyshev  (R=30dB)  and  the  corresponding 
realized  asymmetric  radiation  patterns. 

(a)  The  Chebyshev  and  realized  asymmetric 
array  pattern  Ftc(u)  R=25. 

(b) The  Chebyshev  and  realized  asymmetric 
array  pattern  Ftc(u)  R=18 . 
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Table  3 . 3 The  characteristics  of  the  broadside  narrow 
beamwidth  radiation  patterns  excited  by  the 
Dolph-Chebyshev  and  the  two  corresponding 
analytic  excitations. 


type  of 
pattern 

HPBW 

(Rad) 

BWnul! 

(Rad) 

MLWR 

(None) 

ri 

(None) 

D SLL  k SLKeak 

(none)  (dB)  (dB) 

Fc(u)  {30} 

.330 

.880 

— 

.869 

18.242  -30.00  -30.00 

Ftc(u)R=25 

.310 

.789 

1.13 

.912 

19.165  -18.06  -25.96 

Ftc(u)R=18 

.287 

.691 

1.38 

.940 

19.758  -19.61  -19.84 

Fc(u)  {26} 

.314 

.793 

.913 

19.163  -26.00  -26.00 

the  realized  asymmetric  patterns  is  the  higher  sidelobe 
scenario.  Furthermore,  the  advantages  gained  from  the  analytic 
excitation  really  depend  on  the  MLWR  parameter,  the  Main-Lobe 
Width  Ratio,  and  the  advantages  are  limited  by  the  required 
sidelobe  ratio  in  the  realized  asymmetric  radiation  pattern. 

Very  often,  an  asymmetric  sidelobe  scenario,  such  as 
( 18dB/26dB)  sidelobes  of  Ftc(u)  R=25  in  Table  3.3,  can  be 
acceptable.  Thus,  a fair  performance  comparison  is  assumed  to 
be  made  between  the  Ftc(u)  R=25  and  the  symmetric  broadside 
Dolph-Chebyshev  pattern  with  the  sidelobe  ratio  R=26  dB,  whose 
characteristics  are  listed  in  the  last  row  of  Table  3.3  as 
Fc(u){26}.  Comparing  the  performance  of  these  two  patterns 
indicats  that  the  advantages  gained  from  an  analytic 
excitation  seem  to  be  modest.  However,  the  practical 
consideration  of  the  applications  of  an  asymmetric  sidelobe 
radiation  pattern  is  not  merely  based  on  these  modest 
advantages. 

From  the  simulation  results  of  Example  1 and  2,  the  most 
unfavorable  characteristic  in  the  asymmetric  radiation  pattern 
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corresponding  to  an  analytic  excitation  is  the  peak  sidelobe 
generated  by  the  mainbeam  mis-match  which  results  from  the 
periodic  convolution  in  the  realization  of  an  analytic 
excitation.  Theoretically,  the  larger  the  number  of  array 
elements  N,  the  narrower  the  mainlobe  width  of  the  spectrum 
of  the  truncation  window.  Therefore,  the  number  of  array 
elements  N could  be  used  to  improve  this  mainbeam  mis-match 
phenomenon.  It  is  interesting  and  worthwhile  to  check  this 
theoretical  concept  by  letting  N be  an  independent  variable. 
Let's  again  assume  that  the  sidelobe  ratio  R=30dB  Dolph- 
Chebyshev  array  system  is  the  original  real  excitation,  and 
the  truncation  window  used  to  realize  the  corresponding 
analytic  excitation  is  Dolph-Chebyshev  window  with  R=25dB, 
and  the  number  of  array  elements  N has  changed  from  11  to  25. 
The  two  sides  of  the  peak  sidelobe  of  the  realized  asymmetric 
radiation  pattern  as  function  of  the  number  of  array  elements 
N are  plotted  in  Figure  3.10. 

In  Figure  3.10,  the  upper  curve  represents  the  peak 
sidelobe  in  the  sidelobe  region  of  u < 0,  which  depends  on  the 
extent  of  the  mainbeam  mis-match  mentioned  earlier,  and  the 
lower  curve  represents  the  peak  sidelobe  in  the  sidelobe 
region  of  u > 0,  which  depends  on  the  peak  sidelobes  of  both 
original  symmetric  pattern  and  realizable  Hilbert  transformed 
pattern.  It  can  be  seen  from  Figure  3.10  that  both  sides  of 
the  peak  sidelobe  are  almost  independent  of  the  number  of  the 
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Nurrber  of  Array  Elements  (fQ 
□ SLLpeak  + SLRpeak 

Figure  3.10  The  peak  sidelobe  ratios  of  the  asymmetric 

radiation  patterns  as  function  of  the  number 
of  array  elements. 

array  elements  N.  The  reason  for  the  independence  of  peak 
sidelobe  with  respect  to  N for  the  case  of  u>0  (lower  curve) 
is  simply  that  we  keep  the  peak  sidelobe  constant  in  both 
original  symmetric  pattern  and  the  spectrum  of  the  truncation 
window  and,  consequently,  in  the  realized  Hilbert  transformed 
pattern  as  well.  However,  for  the  case  of  u<0  (upper  curve) 
the  undesired  peak  sidelobe  is  independent  of  N because  as 
we  increase  the  number  of  the  array  elements,  the  mainlobe 
width  of  the  spectrum  of  the  truncation  window  becomes 
narrower.  However,  the  mainbeam  width  of  the  pattern  of  the 
original  real  excitation  is  narrower  as  well,  that  is,  the 
MLWR  parameter  almost  remains  constant  with  the  increase  in 
the  number  of  the  array  elements.  Therefore,  the  increase  in 
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the  number  of  the  array  elements  does  not  improve  the  mainbeam 
mis-match  phenomenon  in  the  realization  of  the  analytic 
excitation  for  an  array  antenna  system. 

3 . 4 Discussion  and  Summary 

Based  on  above  qualitative  and  quantitative  analysis 
results,  we  could  reach  several  statements  related  to 
radiation  pattern  of  an  array  antenna  system  using  analytic 
excitation.  They  are  stated  as  follows: 

1.  The  far  field  radiation  pattern  produced  by  an 
analytic  excitation  is  inherently  asymmetric,  that  is,  the 
analytic  signal  theory  can  provide  a generalized  analytic 
method  to  design  a practical  desired  asymmetric  far  field 
radiation  pattern  of  an  array  antenna  system. 

2.  In  general,  the  asymmetric  radiation  pattern  of  an 
analytic  excitation  has  the  advantages  on  every  performance 
indices  over  the  corresponding  original  broadside  narrow 
beamwidth  symmetric  radiation  pattern,  except  the  peak 
sidelobe  scenario.  Therefore,  a required  low  peak  sidelobe  of 
radiation  pattern  is  implausible  using  a realizable  analytic 
excitation. 

3.  By  adjusting  the  Main-Lobe  Width  Ratio  (MLWR) , and  the 
peak  sidelobe  of  the  truncation  window,  we  can  obtain  a design 
result.  The  larger  the  value  of  MLWR,  the  more  advantages  we 
obtain  in  the  characteristics  of  HPBW,  BWnu11,  directivity, 
radiation  efficiency,  and  the  lower  the  undesired  peak 
sidelobe  generated  by  the  mainbeam  mis-match  which  is 
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consequence  of  the  periodic  convolution  in  the  realization  of 
an  analytic  excitation.  However,  the  larger  the  value  of  MLWR, 
the  higher  the  overall  sidelobe  level  of  the  truncation  window 
as  well  as  the  realized  asymmetric  radiation  pattern. 
Therefore,  the  advantages  are  taken  in  the  compromise  between 
those  pattern  characteristics. 

4.  The  increase  in  the  number  of  array  elements  (N)  does 
not  affect  the  peak  sidelobe  scenario  of  the  asymmetric 
radiation  pattern  of  an  analytic  excitation,  because  the 
parameter  MLWR  is  normally  independent  from  the  number  of 
array  elements  N. 

5.  The  cause  for  producing  the  undesired  peak  sidelobe 
in  the  realization  of  an  analytic  excitation  explained  earlier 
can  be  attributed  to  the  realizable  Hilbert  transform  that  can 
not  faithfully  transform  a radiation  pattern.  For  the  cases 
where  the  design  of  a desired  radiation  pattern  can  be  started 
from  a symmetric  radiation  pattern  which  has  small  amount  of 
components  in  the  u=0  direction,  the  utilization  of  the 
analytic  excitation  will  be  found  to  be  valuable. 

The  information  extracted  from  the  analysis  results  in 
this  chapter  provides  both  essential  and  valuable  foundation 
for  designing  several  practical  array  radiation  patterns  using 
the  analytic  excitation,  which  will  be  developed  in  the 
following  chapters. 


CHAPTER  IV 


THE  BROADSIDE  NARROW  BEAMWIDTH  ASYMMETRIC  SIDELOBE  ARRAY 

RADIATION  PATTERN  DESIGN  USING  ANALYTIC  EXCITATION 

4 . 1 Introduction 

To  overcome  the  disadvantageous  property  of  the  Hilbert 
transformation  in  the  realization  of  an  analytic  excitation 
in  the  preceeding  chapter,  in  this  chapter,  we  exploit  the 
utilization  of  the  analytic  excitation  which  will  generate 
useful  broadside  narrow  beamwidth  asymmetric  sidelobe 
radiation  patterns  for  an  array  antenna  system. 

The  required  shape  of  an  antenna  radiation  pattern  may 
very  much  depend  on  applications.  There  are  many  environmental 
situations  where  the  asymmetric  sidelobe  radiation  pattern  is 
desired  [28-34].  However,  the  design  of  an  asymmetric  sidelobe 
narrow  beamwidth  radiation  pattern  is  not  as  straightforward 
as  a symmetric  radiation  pattern.  An  iterative  try  and  error 
numerical  method  is  commonly  adopted  in  dealing  with  such 
problems.  This  chapter  will  first  present  a new  gneralized 
design  concept,  which  combines  the  analytic  signal  concept  and 
modulation  theory,  for  designing  the  desired  broadside  narrow 
beamwidth  asymmetric  sidelobe  radiation  patterns  of  an  array 
antenna  system.  Then,  the  design  criteria  and  error  analysis 
for  the  physical  implementation  of  this  generalized  design 
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concept  will  be  carried  out.  Several  practical  design  examples 
will  show  that  the  asymmetric  sidelobe  radiation  patterns  have 
various  advantages  over  the  symmetric  radiation  patterns,  and 
the  features  of  this  design  technique  will  be  generalized. 
Finally,  we  will  show  that,  by  using  the  2-D  analytic 
excitations  and  this  design  concept,  an  asymmetric  sidelobe 
planar  array  radiation  pattern  can  be  obtained. 

4.2.  The  Fundamental  Theories  for  Designing  a Broadside  Narrow 

Beamwidth  Asymmetric  Sidelobe  Array  Radiation  Pattern 

Knowing  that  the  realizable  discrete-time  Hilbert 
transformer  can  not  handle  the  patterns  which  have  a large 
amount  of  components  around  u=0  direction,  we  will  derive  a 
generalized  modulated  analytic  excitation  which  is  generated 
by  combining  analytic  excitation  and  modulation  theory,  so 
as  to  shift  the  typical  broadside  narrow  beamwidth  radiation 
pattern  off  the  u=0  and  to  reduce  the  distortion  caused  by 
the  realizable  Hilbert  transfomation.  As  a result,  we  can 
generate  a broadside  narrow  beamwidth  radiation  pattern  which 
has  a practical  desired  asymmetric  sidelobe  scenario. 

Let's  ideally  assume  that  Jr(n)  is  a real  excitation  for 
an  N elements  broadside  linear  array  system,  and  its 
corresponding  narrow  beamwidth  radiation  pattern,  denoted  as 
Fr(u),  has  a negligible  response  in  the  range  um<|u|<7r 
(approximation  to  a beamlimited  pattern) . Then,  in  order  to 
shift  the  pattern  Fr(u)  off  the  u=0  direction,  a modulated 
real  excitation  Jmr(n)  can  be  formed  as 
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Jmr(n)=Jr(n)  s(n)  = Jr  (n)  COS  (ucn)  (4.1) 

where  s (n) =cos (ucn)  represents  a generalized  sinusoidal 
excitation,  where  uc  is  the  modulation  angle  in  radians.  Due 
to  the  periodicity  of  the  array  pattern,  the  following 
condition  is  also  assumed  to  be  satisfied 

[uc  + uj  < n (4.2) 

Taking  the  discrete-time  Fourier  transform  on  both  sides  of 
Eguation  (4.1),  and  applying  the  modulation  theory,  the 
radiation  pattern,  denoted  as  Fmr(u),  of  Jmr(n)  is 

Fmr(u)  = (1/2)  [Fr(u-uc)+Fr(u+uc)  ] (4.3) 

Equation  (4.3)  shows  that  by  applying  modulation  theory,  we 
can  shift  a typical  broadside  narrow  beamwidth  radiation 
pattern  away  from  u=0  direction.  However,  in  order  to 
eliminate  one  of  the  shifted  patterns  [Fr(u-uc)]  or  [Fr(u+uc)] 
presented  in  the  visible  range,  we  can  apply  the  analytic 
signal  theory.  Let's  form  a generalized  modulated  complex 
analytic  excitation  Jma(n)  as 

Jma(n)  = Jmr(n)  + jJmi(n)  (4.4) 

where  Jm1(n)  is  the  discrete-time  Hilbert  transform  of  Jmr(n)  . 
Thus,  Fmi(u)  denotes  the  pattern  of  Jm1(n),  and  according  to 
Equation  (2.15) , 


-3Fmr(u)  , 
jFmr(u)  / 


0 < U < 7T 
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Fm1(u)  = 


-n  < u < 0 


Consider  the  facts  that  the  pattern  Fr(u)  corresponding 
to  Jr(n)  is  assumed  to  be  beamlimited,  and  the  value  of  uc  is 
set  as  Equation  (4.2).  Those  assumptions  assure  that  the 
pattern  Fmi(u)  can  be  expressed  as1 

Fmi(u)  = (1/2)  [-jFr(u-uc)  + j Fr (u+uc)  ] (4.5) 

By  taking  Inverse  Fourier  transform  of  Equation  (4.5),  the 
excitation  Jmi(n)  becomes 

Jmi(n)  = (1/2)  [-jjr(n)  <?jucn  + jjr(n)e-3ucn 

= Jr (n)  sin (ucn)  (4.6) 

Then,  the  complex  analytic  excitation  in  Equation  (4.4) 
becomes 

Jma(n)  = Jmr(n)  + jJmi(n) 

= Jr (n)  [cos (ucn)  + jsin(ucn)] 

= Jr(n)  e^ucn  (4.7) 

By  taking  the  discrete-time  Fourier  transformation,  the  narrow 
beamwidth  radiation  pattern,  Fma(u),  of  the  modulated  analytic 
excitation  Jma(n)  in  Equation  (4.7)  is 


The  error  caused  by  the  non-beamlimited  (sidelobe 
component)  of  a typical  array  radiation  pattern  will  be 
analyzed  later. 
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Fma(U)  = Fr(U"Uc)  (4.8) 

It  is  a single  narrow  beamwidth  radiation  pattern  which 
remains  within  the  visible  range.  Then,  by  using  the  phase 
shifting  property  of  an  array  antenna  system,  the  shifted 
pattern  in  Eguation  (4.8)  can  be  shifted  back,  such  that  the 
boresight  of  the  realized  asymmetric  sidelobe  radiation 
pattern  is  pointed  to  the  u=0  direction.  That  means,  by 
controlling  the  parameter  uc,  we  can  obtain  a half-sided 
sidelobe  radiation  pattern  (which  is  eguivalent  to  saying  that 
one  of  the  two  sides  sidelobes  can  be  eliminated)  without 
sacrificing  the  system  performance  of  an  original  broadside 
symmetric  radiation  pattern,  provided  that  the  sidelobe  ratio 
of  the  original  radiation  pattern  is  reasonably  low  enough. 
Thus,  the  analytic  signal  theory  really  provides  a generalized 
theory  for  designing  desired  broadside  narrow  beamwidth 
asymmetric  sidelobe  radiation  patterns  of  an  array  antenna 
system. 

A generalized  theory  has  been  found  to  design  narrow 
beamwidth  asymmetric  sidelobe  radiation  patterns,  and  will  be 
explained  in  the  following  sections. 

4 . 3 System  Design  Constraint 

Using  a modulated  analytic  excitation,  a portion  of  the 
half-side  sidelobe  of  a broadside  symmetric  radiation  pattern 
can  be  eliminated.  The  extent  of  sidelobe  elimination  depends 
on  the  value  of  the  modulation  angle  uc. 


However,  from 
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practical  application  points  of  view,  the  purpose  of  using 
such  a modulated  analytic  excitation  is  to  reduce  one  of  the 
two  sides  sidelobes  in  an  array  radiation  pattern  without 
sacrificing  the  system  characteristics.  Thus,  the  system 
design  constraint  is  that  the  value  of  uc  is  restricted  to  be 
exactly  equal  to  half  of  the  mainbeam  width  of  the  original 
broadside  symmetric  pattern  (uc=BWnull/2 ) . Thus  the  basic 
concept  of  this  design  methodology  is  simply  to  shift  the 
mainbeam  of  the  radiation  pattern  to  the  pass-band  of  the 
spectrum  of  the  realizable  Hilbert  transformer,  and  to 
position  the  first  null  of  the  pattern  on  the  u=0,  so  as  to 
reduce  the  distortion  caused  by  the  realizable  Hilbert 
transformer  about  u=0,  and  to  eliminate  the  sidelobe  left  on 
the  reverse. 

4 . 4 System  Analysis  and  Performance  Evaluation 

With  the  constraint  value  of  uc,  there  are  two  factors 
that  will  determine  the  achievements  of  this  design 
methodology.  They  are  presented  and  analyzed  as  follows: 

Consider  a typical  real  even  symmetric  excitation  Jr(n) 
and  its  corresponding  broadside  symmetric  narrow  beamwidth 
radiation  pattern  Fr(u)  which  is  depicted  in  Figure  4.1(a). 
The  modulated  real  excitation  is  then  Jmr  (n)  =Jr  (n) cos (ucn)  and 
its  radiation  pattern  becomes  Fmr(u)  = (l/2)  [Fr  (u-uc) +Fr  (u+uc)  ] . 
The  pattern  of  Fmr(u)  with  the  modulation  angle  uc=  (BWnu11/2 ) 
can  be  further  decomposed  into  mainbeam  pattern  and  sidelobe 
pattern  as  depicted  in  Figure  4.1(b)  and  (c) . 
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Figure  4.1  The  decomposed  mainbeam  and  sidelobe  patterns 
of  a real  modulated  excitation. 

(a)  The  real  symmetric  pattern  Fr(u) 

(b) The  Fmrl(u)  and  its  decomposed  patterns. 

(c) The  Fmr2(u)  and  its  decomposed  patterns. 
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Thus, 


Fmr(u)  = (1/2)  [Fr(u-uc)  + Fr(u+uc)] 

= (1/2)  [Fmrl  (u)  + Fmr2  ( u ) ] 

= (l/2)[Fj;(u)  + Fm~(u)  + fJ2(u)  + Fm~(u)] 

(4.9) 


Where  Fmrl(U)=Fr(U“Uc)  ' Fmr2  (u)  =Fr  (U+Uc)  , and 

0,  0<U<7T 

Fmrl  (U)  ' -7r^u<0 


Fil<u>- 


Fmrl(u),  0<U<7T 

0,  -n<u<0 


F»rl(u>  = 


Ftr2<U) 


Fmr2(u),  0<U<7T 
0,  -7T<U<0 


0<U<7T 

-n<u<0 


Then,  the  radiation  pattern  Fmi(u)  corresponding  to  the 
excitation  which  is  the  Hilbert  transform  of  Fmr(u)  is 

Fmi  (u)  = -jsgn(u)Fmr(u) 

- M/2 ) [ Kn  (u>  +Ft2  (u)  ] + ( j/2 ) [ F~rl  (u)  +F“  2 (u)  ] 

(4.10) 

From  Equation  (4.9)  and  (4.10),  the  resulting  asymmetric 
sidelobe  radiation  for  the  corresponding  realizable  modulated 
analytic  excitation,  Fma(u),  becomes 
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Fma(u)  = Fmr(u)  + jFmi(u) 

= (l/2)[Ft1(u)  + F;rl(u)  + Ftr2(u)  + F;r2(u)] 

+ (l/2)[F+1(u)+Ft2(u)]  - (1/2)[F"1(U)+F"2(U)] 

= F+rl(u)  + F+r2(u)  (4.11) 

The  first  term  on  the  right-hand  side  of  Equation  (4.11) 
represents  the  desired  asymmetric  sidelobe  radiation  pattern, 
and  the  second  term  is  the  undesired  sidelobe  which  leaks  into 
the  desired  asymmetric  sidelobe  radiation  pattern.  Thus, 
because  of  the  sidelobe  leakage  effect,  the  first  factor  that 
influences  the  accomplishments  of  this  design  method  is  the 
sidelobe  level  of  the  original  broadside  symmetric  pattern. 
The  distortion  caused  by  this  sidelobe  leakage  effect 
indicates  that  the  lower  the  sidelobe  in  the  original 
symmetric  radiation  pattern  Fmr(u)  , the  more  advantages  we  gain 
in  the  resulting  asymmetric  sidelobe  radiation  pattern.  The 
second  factor  is  that  the  small  amount  of  components  in  the 
shifted  pattern  Fmr(u)  around  the  u=0  will  still  be  distorted 
by  a realizable  Hilbert  transformer.  However,  due  to  the 
shifting  of  Fmr(u),  the  amount  of  distortion  is  greatly 
reduced.  Although  the  distortion  appears  to  be  unavoidable, 
it  will  be  shown  that  the  generation  of  an  asymmetric  sidelobe 
radiation  pattern  using  a realizable  modulated  analytic 
excitation  is  both  practical  and  valuable. 
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In  evaluation  of  the  performances  of  the  asymmetric 
sidelobe  radiation  pattern,  in  addition  to  these 
characteristic  indices  defined  in  Section  3.3,  Chapter  III, 
there  are  two  more  characteristic  indices  which  are  required 
to  assess  the  figure  of  merits  of  the  asymmetric  sidelobe 
radiation  pattern  compared  to  the  original  broadside  symmetric 
radiation  pattern  of  a linear  array  antenna  system.  These  are 
defined  as  follows  and  the  related  parameters  can  be  refered 
to  Figure  4.2. 

a.  The  peak  sidelobe  suppression  gain  (SSGpeak)  is  defined 
as  the  difference  (dB)  of  the  peak  sidelobes  between  the 
original  symmetric  and  the  asymmetric  sidelobe  radiation 
patterns  in  the  supressed  side  (u<0) . Let  u<0  represent  the 
left-hand  side,  the  SSGpeak  can  be  expressed  as 

SSGpeak(dB)  = SLLspegk(dB)  - SLLapeak(dB) 
where  SLLspeak  and  SLLapeak  denote  the  peak  sidelobe  of  the 
symmetric  and  asymmetric  patterns  at  the  left-hand  side  (u<0) . 

b.  The  supressed  radiation  power  ratio  (SRPR)  is  defined 
as  the  ratio  of  the  radiation  powers  of  the  asymmetric 
sidelobe  radiation  pattern  to  the  original  symmetric  radiation 
pattern  in  the  suppressed  side  (u<0) . Mathematically,  it  can 
be  expressed  as 

'anullj 

|Fma(u)  |2  du 

J ~TT 

SRPR  = — (4.12) 

'snullj 

I Fr  (u)  |2  du 

— 7T 

V 
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where  the  Fma(u)  and  Fr(u)  represent  the  asymmetric  sidelobe 
and  the  symmetric  radiation  patterns,  and  anullj  and  snullj  are 
the  first  null  in  the  supressed  side  of  F^fu)  and  Fr(u)  with 
respect  to  the  mainbeam. 


Figure  4.2  The  parameters  used  to  define  the  SSG  k and 
SRPR.  (a) the  SLL  k and  snullj. 

(b)the  SLLapeak  and  anullj 

4 . 5 System  Design  Procedure 

The  design  procedure  of  a broadside  narrow  beamwidth 
asymmetric  sidelobe  radiation  pattern  of  an  array  antenna 
system  using  modulated  analytic  excitation  starts  from  a well- 
developed  real  excitation,  and  can  be  formally  summaried  as: 

(a)  Select  a suitable  real  excitation  Jr(n)  and 
evaluate  its  mainbeam  width  (BWnu11)  . 

(b)  Form  the  real  modulated  excitation  Jmr(n)  as  given 
by  Eguation  (4.1)  with  a modulating  angle 

uc=(BWnull/2). 

(c)  Acguire  the  required  modulated  analytic  excitation 
J^n)  as  given  by  Equation  (4.4)  by  carrying  out  the 
discrete-time  Hilbert  transform. 

(d)  Evaluate  the  system  performance. 
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Since  the  ideal  modulated  analytic  excitation  in 
Equation  (4.7)  is  unobtainable,  as  long  as  the  original 
symmetric  radiation  Fmr(u)  is  not  beamlimited,  the  acquisition 
of  the  modulated  analytic  excitation  Jma(n)  in  step  (c)  should 
be  carried  out  as  follows 

Jma(n)  = Jr(n)cos(ucn)+j  { [Jr(n)cos(ucn)  ]*h(n)  }w(n)  (4.13) 

where  h(n)  is  the  discrete-time  Hilbert  transform  defined  in 
Equation  (2.18),  and  w(n)  is  a truncation  window  function. 
Thus,  the  realized  asymmetric  sidelobe  radiation  pattern, 
Fma(u)/  corresponding  to  the  modulated  analytic  excitation  of 
Jma(n)  in  Equation  (4.13)  becomes 

Fma(u)  = Fnu-(U)  + j [ Fmr  (u)  H (u)  ] *W  (u)  (4.14) 

where  Fmr(u)  is  the  pattern  of  Jr (n)  cos  (ucn)  , and  H(u)  and  W(u) 
are  the  spectrums  of  h(n)  and  w(n)  respectively,  and  * 
indicats  the  discrete-time  convolution.  The  purpose  of  the 
window  function  w(n)  in  realization  of  this  asymmetric 
sidelobe  radiation  pattern  as  given  by  Equation  (4.13)  is  the 
same  as  that  in  Chapter  III.  However,  the  effects  of  the 
truncation  window  on  the  asymmetric  sidelobe  radiation  pattern 
are  different  from  those  discussed  in  Chapter  III.  The  major 
difference  is  that  there  is  no  discontinuity  in  the  shifted 
pattern  Fmr(u) , and  consequently,  there  is  no  roll-off  region 
in  the  pattern  ( [Fmr (u) H (u) ] *W(u) } . However,  because  of  the 
convolution  process,  the  spectrum  of  the  window  function  will 
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utimately  leak  into  the  pattern  [Fmr(u)H(u)  ] , and  result  in  a 
ripple  on  the  waveform  of  [Fmr(u)H(u)  ] . 

After  recalling  the  effects  of  the  truncation  window  in 
the  realization  of  a modulated  analytic  excitation,  we  can 
consider  the  question  of  what  is  the  best  selection  of  the 
truncation  window  under  this  circumstance.  Here,  the  major 
concern  is  to  eliminate  one  of  two  sides  sidelobes  of  an  array 
radiation  pattern.  Furthermore,  the  truncation  window  should 
have  the  negligible  influence  on  the  shape  of  [Fmr (u) *H (u) ] , 
not  only  on  the  mainbeam  but  also  on  the  sidelobe.  In  order 
to  achieve  this  result,  the  more  important  characteristic  of 
the  spectrum  of  the  desired  truncation  window  (in  addition  to 
the  narrow  mainlobe  width)  is  that  the  mainlobe  should  contain 
as  much  of  the  energy  as  possible  and  the  sidelobe  should 
decrease  rapidly  as  u moves  toward  to  n . 2 Thus,  the  solution 
of  the  best  chosen  truncation  window  in  implementation  of  a 
realizable  modulated  analytic  excitation,  in  countrary  to 
Chapter  III,3  converges  to  the  Kaiser  window  [12,57,64]  which 
is  claimed  to  be  an  approximation  to  a class  of  functions 
called  the  prolate  spheroidal  wave  functions  [36]  which 
possess  the  desired  properties  mentioned  above.  The  Kaiser 

2 • • 

This  is  because  a window  function  that  has  this  property 
produces  the  least  amount  of  ripple  on  the  convolved  waveform. 

3The  major  concern  in  Chapter  III  is  to  seek  a window 
function  which  has  minimun  mainlobe  width  given  a peak 
sidelobe  ratio,  in  order  to  reduce  the  roll-off  region  at  the 
discontinuity  of  a spectrum.  That  leads  to  the  Dolph-Chebyshev 
window  function  as  the  best  selection. 


window  is  expressed  as  [57] 
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wk(n) 


I0[6(l-[2n/(N-l)]2)1/2] 

lot*] 


-((N-l)/2)  < n < ( (N-l)/2) 

(4.15) 


where  N is  odd,  and  I0[  • ] represents  the  modified  Bessel 
function  of  first  kind  of  order  zero,  and  B is  a constant  that 
determines  the  tradeoff  between  the  peak  sidelobe  and  the 
mainlobe  width.  The  larger  the  6,  the  lower  the  sidelobe,  and 
the  wider  the  mainlobe  width  of  the  spectrum.  It  should  be 
noted  that  when  6=0,  the  Kaiser  window  becomes  the  uniform 
rectangular  window. 

There  is  an  additional  point  that  needs  to  be  considered. 
If  the  value  of  uc  is  given  as  the  first  null  of  the  original 
selected  symmetric  pattern  with  some  tolerance,  then,  the 
shifted  pattern  Fmr(u)  corresponding  to  Jr(n)  cos  (ucn)  will  have 
component  exactly  at  the  u=0  that  will  produce  a noticable 
distortion  due  to  the  realizable  Hilbert  transformer. 
Therefore,  in  order  to  allow  for  some  error  in  the  value  of 
uc  i-n  the  physical  calculation  of  the  modulated  analytic 
excitation,  the  procedure  used  to  compute  the  modulated 
analytic  excitation  in  Equation  (4.13)  can  be  judiciously 
modified  as 


Jma(n)=Jr(n)  sin(ucn)  + j { [ Jr ( n)  sin (ucn)  ] *h (n)  } w (n)  (4.16) 
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That  is,  instead  of  cos(ucn),  the  real  excitation  Jr(n)  is 
modulated  by  sin(ucn)  . Then,  with  a typical  real  even 
symmetric  excitation  Jr(n)  , Jr(n)  sin(ucn)  is  always  an  odd 
symmetric  excitation  for  any  specific  value  of  uc,  and  its 
corresponding  pattern  will  exactly  have  no  component  at  u=0. 
That  is,  this  modification  assures  that  the  pattern  which  will 
be  transformed  by  a realizable  Hilbert  transformer  has  zero 
component  at  u=0,  and  consequently,  minimize  the  distortion 
caused  by  the  realizable  Hilbert  transformer.  It  is  trivial 
to  show  that  the  modification  of  the  modulated  analytic 
excitation  in  Equation  (4.16)  does  not  change  the  shape  or 
characterictics  of  the  asymmetric  sidelobe  radiation  pattern 
as  indicated  in  Equation  (4.14),  except  a phase  constant  that 
does  not  affect  any  shape  or  performances  of  a desired 
pattern. 

There  are  three  parameters  that  allow  us  to  consider  to 
generate  the  asymmetric  sidelobe  radiation  pattern.  They  are 
the  number  of  the  array  elements  (N) , the  peak  sidelobe  of  the 
original  symmetric  pattern,  and  the  characteristics  of  the 
truncation  window.  Normally,  the  number  of  array  elements  does 
not  affect  the  sidelobe  level,  provided  that  N is  large 
enough.  Thus,  once  the  best  choice  of  window  function  is 
determined,  the  levels  of  the  sidelobes  in  the  realized 
asymmetric  sidelobe  pattern  will  be  controlled  by  the  original 
selected  symmetric  patterns.  A rule  of  thumb  design  strategy 
will  be  demonstrated  in  the  following  design  Example  2. 
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4 . 6 Design  Examples 

Example  1 For  demonstration  of  the  generation  of  an 
asymmetric  sidelobe  radiation  pattern  utilizing  the  modulated 
analytic  excitation,  we  use  the  Hamming  excitation  Jh(n) 
(refer  to  Example  1,  Chapter  III)  as  an  original  real 
symmetric  excitation,  the  number  of  array  elements  is  kept  the 
same  as  that  used  in  Chapter  III,  N=21,  and  the  truncation 
window  is  Kaiser  wk(n)  with  8=0.  First,  the  mainbeam  width  of 
the  array  radiation  pattern  excited  by  Jham(n)  is  computed  to 
be  BWnuii=1* 359  (rad)  • Then,  the  modulated  real  excitation  is 
Jmr  (n)  = [ Jham(n)  sin  (°  • 679n)  ] . Following  the  design  procedures 
described  in  previous  Section  4.5,  the  corresponding 
realizable  modulated  analytic  excitation  Jma(n)  as  given  by 
Eguation  (4.16)  can  be  obtained  by  carrying  out  the  Hilbert 
transform  and  truncation  process.  The  original  symmetric 
Hamming  pattern  and  the  shifted  radiation  pattern  Fma(u+uc) 
(uc=0.679)  corresponding  to  the  realized  modulated  analytic 
excitation  Jma(n)  are  plotted  together  in  Figure  4.3.  By 
examining  the  shape  of  the  pattern  Fma(u+uc)  in  Figure  4.3,  we 
find  that  the  sidelobe  in  the  range  of  u<0  has  been  greatly 
suppressed.  Furthermore,  it  is  interesting  to  compare  the 
asymmetric  pattern  in  Figure  4.3  with  the  asymmetric  pattern 
in  Figure  3.8(b).  Both  patterns  are  excited  by  an  analytic 
excitation  generated  by  using  Hamming  excitation  as  an 
original  real  excitation.  We  find  that  the  undesired  high  peak 
sidelobe  in  the  u<0  range  in  Figure  3.8(b)  disappears  in 
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Figure  4.3.  The  improvement  is  actually  attributed  to  the 
modulation  theory  which  has  been  applied  to  the  analytic 
excitation  that  produces  the  asymmetric  sidelobe  radiation 
pattern  of  Figure  4.3.  In  other  words,  the  disadvantageous 
property  of  an  asymmetric  radiation  pattern  occured  in  Chapter 
III  can  be  overcome  by  applying  modulation  theory. 


Figure  4.3  The  Hamming  Fham(u)  and  the  corresponding 

realized  asymmetric  sidelobe  Fma(u+uc)  radiation 
patterns. 

(N=21,  and  tuncation  window  is  Kaiser  with  6=0) 


Table  4 . 1 The  characteristics  of  the  radiation  patterns 

excited  by  Hamming  and  the  corresponding  modulated 
analytic  excitations.  (N=21) 


Type  of 
pattern 

HPBW 

(Rad) 

BWnun 

(Rad) 

D T)  SLL  k SLR  k 

(None)  (None)  (dlf)  (dfe) 

SSG 

(dB) 

SRPR 

(dB) 

Fham  ( u ) 
Fn,a(U+Uc) 

.401 

.404 

1.359 

1.312 

14.885  .709  -40.64  -40.64 

14.697  .700  -53.07  -40.29 

-12.4 

12.4 

41.6 

-41.6 
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Table  4.1  lists  the  characteristics  of  the  radiation 
patterns  for  which  Fham(u)  represents  the  pattern  excited  by 
Hamming  function  and  Fma(u+uc)  represents  the  pattern  of  the 
realized  asymmetric  sidelobe  radiation  pattern.  By  comparing 
the  quantitative  results  of  the  two  patterns  listed  in  Table 
4.1,  we  find  that  those  system  performances  of  the  original 
symmetric  pattern  are  almost  maintained  in  the  realized 
asymmetric  sidelobe  pattern,  except  the  peak  sidelobe  at  the 
left-hand  side  region  (SLLpeak)  , where  the  asymmetric  sidelobe 
radiation  pattern  has  tremendous  advantage  over  the  original 
symmetric  pattern  by  lower  than  12.4dB  (SSGpeak=12 . 4dB)  . 
Moreover,  the  suppressed  radiation  power  (SRPR)  in  the 
supressed  side  (u<0)  of  the  asymmetric  sidelobe  pattern  is 
41dB  lower  than  that  of  the  original  symmetric  pattern.  This 
characteristic  shows  one  of  the  reasons  why  the  asymmetric 
sidelobe  radiation  pattern  is  attractive  for  the  areas  where 
interferences  come  from  wide  range  of  angle  in  the  space. 
Moreover,  this  example  demonstrates  the  fact  that  the 
foregoing  analytic  signal  theory  really  can  provide  an 
effective  method  to  design  a desirable  asymmetric  sidelobe 
radiation  pattern.  In  the  next  design  example,  we  will  show 
a design  strategy  for  generating  a practical  and  useful 
asymmetric  sidelobe  radiation  pattern. 

Example  2 In  this  design  example,  by  using  the  sampled 
Taylor  aperture  distribution,  we  will  show  how  to  design  an 
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asymmetric  sidelobe  radiation  pattern  using  modulated  analytic 
excitations  that  are  practical  and  useful.  Then,  the 
performance  comparisons  are  made  between  the  symmetric 
radiation  pattern  excited  by  the  sampled  Taylor  aperture 
distribution  and  the  realized  aymmetric  sidelobe  radiation 
pattern.  The  reasons  for  choosing  the  sampled  Taylor  aperture 
distribution  in  this  design  example  are  that  it  not  only 
provides  a flexible  method  to  obtain  various  sidelobe  levels 
of  array  radiation  patterns  but  also  is  widely  adopted  in  the 
practical  realization  of  an  array  antenna  system,  because  it 
has  higher  directivity  as  well  as  excitation  (radiation) 
efficiency  than  the  Dolph-Chebyshev  array  system. 

It  was  recognized  that  the  continuous  counterpart 
(aperture  distribution)  of  the  Dolph-Chebyshev  array  system 
is  physically  unrealizable  [20].  By  changing  the  locations  of 
the  zeros  of  the  radiation  pattern,  Taylor  was  able  to  modify 
the  uniform  aperture  distribution  to  approximate  the  Dolph- 
Chebyshev  aperture  distribution.  The  Taylor  pattern  can  be 
expressed  as  [20] 


F(z , A, n) 


sin  7rz 

( ) 

nz 


n-1  z2 

II  (1 } 

n=l  ct2[A2  + (n-  (1/2 ) ) 2] 


n-1 

n 

n=l 


(1- 


z 


2 


) 


n 


2 


(4.17) 
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where  z=(L/A)cos0,  with  L the  length  of  the  aperture  in  units 
of  wavelength,  6 the  angle  measured  from  the  endfire 
direction,  and  \ is  the  wavelength.  The  shape  of  the  Taylor 
pattern  is  clearly  understandable  by  examining  Equation 
(4.16).  The  term  (sin  nz/nz)  in  the  numerator  is  the  typical 
uniformly  excited  aperture  pattern,  which  has  zeros  at  ±n, 

n=l,  2 , 3 , • • • ,oo.  Then,  the  first  (n  -1)  pairs  of  inner  nulls 

of  (sin  nz/nz)  are  to  be  modified,  so  that  the  first  (n  -1) 
pairs  of  inner  peak  sidelobes  of  the  pattern  (sin  nz/nz)  are 
converted  to  the  desired  Dolph-Chebyshev  equal  peak  sidelobe 
scenario.  Thus,  the  product  term  in  the  denominator 


n-1 

n 

n=l 


(1- 


z 


2 


is  used  to  remove  the  first  (n  - 1)  pairs  of  inner  zeros  of 
function  (sin  nz/nz) . These  removed  nulls  are  replaced  by  the 
zeros  of  the  product  term  in  the  numerator 


n-1  z2 

n {1 } 

n=l  a2[A2+  (n-(l/2) ) 2] 


(4.18) 


which  is  derived  from  the  Dolph-Chebyshev  array  pattern  by 
letting  the  number  of  array  elements  to  be  indefinitely 
increased,  where  A is  an  adjustable  real  parameter  such  that 
(cosh  7rA)=b,  with  b the  desired  voltage  sidelobe  ratio,  and 
a is  defined  by 
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n 

a = 

[A2  + (n-  ( 1/2 ) ) 2] 1/2 

The  aperture  distribution  corresponding  to  the  Taylor 
pattern  as  given  by  Eguation  (4.17)  can  be  obtained  by  using 
the  Fourier  transform  relations  between  the  aperture 
distribution  and  the  far  field  pattern.  After  doing  so,  the 
aperture  distribution,  J(x) , is  found  to  be 

n-1 

J(x)  = ( 1/L)  { 1 + 2 E [F(m,  A,n)  cos  ( 2?rmx/L)  ] ) , |x|  < (L/2) 

m=l 

(4.19) 

where  F(m,A,n)  is  defined  in  Equation  (4.17),  with  z replaced 
by  m.  From  Equation  (4.17),  the  shape  of  the  Taylor  pattern 

is  fully  specified  by  two  parameters,  n and  A.  The  required 

value  of  n for  different  levels  of  sidelobes  can  be  obtained 

from  Taylor  [20,  p.23.  Table  II].  In  this  design  example,  n 
is  set  to  4,  which  is  chosen  to  cover  the  sidelobe  levels  from 
OdB  to  30dB,  unless  otherwise  specified. 

Applying  the  aperture  sampling  method  [65-66],  the 
continuous  Taylor  aperture  distribution  can  be  applied  to  an 
array  antenna  as  well.  Consider  an  array  of  (N+l)  (N  is  even) 
elements  with  interelement  distance  d whose  unit  is 
wavelength,  then  the  length  of  an  aperture  is  assumed  to  be 
L=(N+l)d,  so  that  the  length  of  the  aperture  extends  a 
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distance  (d/2)  over  the  two  ends  of  the  array.  Thus,  the  array 
excitation  Jr(n)  corresponding  to  the  Taylor  aperture 
distribution  are  the  sampled  values  of  J(x),  that  is 
Jr(n)=J(xn),  where 

xn=  nd,  n = 0,±1,±2,±3,  • • • , ± (N/2 ) (4.20) 

Using  the  sampled  Taylor  aperture  distribution,  and  with  the 
number  of  array  elements  N=21,  a Taylor  pattern  for  the  case 

of  25dB  peak  sidelobe  with  n=4  is  plotted  in  Figure  4.4(a). 

With  the  same  parameters  of  N and  n,  the  BWnull  of  the  Taylor 
patterns  with  peak  sidelobes  from  20dB  to  30dB  are  plotted  in 
Figure  4 . 4 (b) . 

Having  acquired  the  Taylor  excitation,* * 4  we  are  able  to 
design  various  desired  asymmetric  sidelobe  radiation  patterns. 
Let's  keep  the  number  of  array  elements  to  be  21.  The  design 
goal  here  is  assumed  to  meet  the  practical  applications.  That 
is,  the  levels  of  the  peak  sidelobes  on  both  sides  of  the 
radiation  pattern  vary  from  20dB  to  35dB.  However,  it  allows 
the  two  sides  of  the  sidelobes  to  be  asymmetric.5  After  an 
experimental  judgement,  we  could  start  to  establish  the  design 


Hereafter,  the  Taylor  excitation  will  be  used  to 

represent  the  sampled  Taylor  aperture  distribution  for  an 

array  antenna  system. 

5 

The  most  practical  case  is  that  one  of  the  sidelobes  be 
about  20dB,  and  the  other  side  to  be  variable.  However,  the 
lower  the  sidelobe,  the  poorer  the  excitation  (radiation) 
efficiency.  Therefore,  a very  low  sidelobe  level  is  also 
impractical . 
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u (radian) 


(a) 


(b) 

Figure  4.4  The  25  dB  Taylor  pattern  and  mainbeam  widths. 

(a)  The  Taylor  pattern  with  N=21,  n=4 , SLJeak=25dB. 

(b)  The  BWnull  of  the  Taylor  pattern  with  N=21, 

n=4,  and  SLpeak  from  20dB  to  30  dB. 
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information  which  contains  the  relations  between  the  levels 
of  the  asymmetric  sidelobes  (SLLpeak  or  SLRpeak)  and  the 
symmetric  Taylor  patterns  with  various  peak  sidelobes.  That 
is,  by  changing  the  sidelobes  of  the  Taylor  patterns  from 
21dBto  27dB,  and  applying  the  procedures  for  generating  the 
corresponding  modulated  analytic  excitation  and  the  asymmetric 
sidelobe  patterns  described  in  Section  4.5,  for  which  the 
required  values  of  uc  can  be  obtained  from  Figure  4.4(b),  and 
Kaiser  truncation  window  with  B=0,6  the  various  asymmetric 
sidelobe  levels  (SLLpeak  and  SLRpeak)  corresponding  to  the 
symmetric  Taylor  patterns  can  be  obtained  and  are  plotted  in 
Figure  4.5,  where  SLLpeak  and  SLLpeak  symbolize  the  left  and 
right  side  of  the  peak  sidelobes  in  unit  of  dB  of  the 
asymmetric  sidelobe  radiation  pattern. 

From  the  system  design  point  of  view,  we  should  be  able 
to  acquire  the  array  excitation  which  produces  the  desired 
asymmetric  sidelobe  radiation  pattern.  Let's  use  a practical 
requirement  as  an  example  to  show  how  to  obtain  the  array 
excitation  of  an  asymmetric  sidelobe  radiation  pattern.  Assume 
that  the  requirements  of  the  sidelobe  levels  of  the  radiation 
pattern  of  an  array  antenna  system  are  (SLLpeak=35dB)  and 
(SLRpeak>2  0dB)  . Then,  from  the  results  plotted  in  Figure  4.5, 


c 

We  can  also  vary  R to  control  the  level  of  the 
asymmetric  sidelobe.  However,  the  larger  the  B,  the  more 
distortion  the  asymmetric  pattern.  B=0  shows  the  best  result 
in  the  realization  of  a practical  useful  modulated  analytic 
excitation. 
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Side  lobe  ([Taylor  patter  nXdET) 

□ SLRpeak  + SLLpaak 

Figure  4.5  The  relations  of  the  asymmetric  sidelobes 

SLL  k and  SLRpe;ak  and  the  Taylor  patterns  with 
various  peak  sidelobes. 

we  find  that  using  27dB  peak  sidelobe  of  Taylor  pattern  as  an 
original  symmetric  pattern  and  the  Kaiser  truncation  window 
with  6=0  can  generate  an  asymmetric  radiation  pattern  whose 
sidelobe  levels  satisfy  the  system  reguirements . Thus,  the 
Taylor  excitation  Jr(n)  can  be  calculated  by  using  Eguation 
(4.19)  and  (4.20).  The  corresponding  modulated  analytic 
excitation  Jma(n)  is  then  computed  from  Equation  (4.16)  with 
uc=0 . 424 (radian)  which  is  obtained  from  Figure  4.4(b)  (half 
of  the  BWnull  of  27dB  Taylor  pattern)  , and  the  truncation 
window  is  Kaiser  with  6=0.  The  resulting  shifted  (35/23) 
asymmetric  sidelobe  radiation  pattern7  and  the  35  dB  Taylor 


7Let  (SLLpeak/SLRpeijk)  be  used  to  symbolize  the  left-hand 
side  and  right-hand  side  peak  sidelobes  with  respect  to  the 
main  beam  of  an  arbitrary  radiation  pattern.  Thus,  (35/23) 
indicates  the  35dB  and  23dB  peak  sidelobes  at  the  left-hand 
and  right-hand  sides  of  an  asymmetric  sidelobe  radiation 
pattern  respectivity . 


117 


pattern  (N=21,  n=5)  are  plotted  together  in  Figure  4.6  for 
comparison  (will  be  discussed  later) , and  the  analytic 
excitation  for  the  (35/23)  asymmetric  sidelobe  radiation 
pattern  is  tabulated  in  Table  4.2.  For  the  shifted  (35/23) 
asymmetric  sidelobe  radiation  pattern  as  plotted  in  Figure 
4.6,  the  phase  in  each  array  excitation  is  the  sum  of  the 
tabulated  value  and  nuc,  where  uc=0 . 424 (radian) . The 
quantitative  characteristics  of  the  realized  shifted  (35/23) 
asymmetric  sidelobe  radiation  pattern  and  the  35  dB  Taylor 
pattern  are  listed  in  Table  4.3. 


u (radian) 

□ Taylor  pattern  + ASYMMETRIC  PATTERN 


Figure  4.6  The  realized  shifted  (35/23)  asymmetric  sidelobe 
and  the  Taylor  symmetric  sidelobe  (35dB) 
radiation  patterns. (N=21) . 
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Table  4.2  The  analytic  excitation  for  (35/23) 

asymmetric  sidelobe  radiation  pattern. 
[N=2 1 , Jma(-n)=J 


array  number 

^ma 

(n) 

n 

amplitude 

phase 

0 

1.019 

90.000 

1 

.988 

-65.715 

2 

.964 

-42.242 

3 

.891 

-17.021 

4 

.809 

5.891 

5 

.718 

31.796 

6 

.594 

54.641 

7 

.515 

80.217 

8 

.381 

-74.164 

9 

.349 

-50.966 

10 

.282 

4.958 

Table  4.3  The  characteristics  of  the  shifted  (35/23) 

asymmetric  sidelobe  and  (35/35)  Taylor  symmetric 
sidelobe  radiation  patterns. 


type  of 
pattern 

HPBW 

(Rad) 

BWnull 

(Rad) 

D 

(None) 

n 

(None) 

SLIw 

(dB) 

SLRpeak 

(dB) 

(35/23) 
Taylor (35) 

.334 
. 353 

. 884 
.997 

18.114 

16.971 

.863 

.808 

-35.045 

-35.046 

-23 . 629 
-35.046 

Let's  assume  that  the  requirement  of  a broadside  array 
antenna  system  on  the  sidelobes  can  be  stated  as  follows:  one 
side  of  the  sidelobes  is  required  to  achieve  the  level 
indicated  as  SLLpeak,  while  the  other  side  of  the  peak 
sidelobes  can  be  a specific  level  which  is  higher  than  SLLpeak. 
That  is,  the  peak  sidelobe  in  left-hand  side  with  respect  to 
the  mainbeam  is  required  to  be  lower  than  that  of  right-hand 
side.  Further  assume  that  the  (SLLpeak/SLRpeak)  asymmetric 
radiation  pattern  can  satisfy  the  system  requirement.  Thus, 
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in  order  to  satisfy  the  system  peak  sidelobe  requirements,  the 
symmetric  Taylor  pattern  should  be  designed  with  the  symmetric 
peak  sidelobe  of  SLLpeak.  The  fair  comparison  is  then  made 
between  the  symmetric  Taylor  pattern  with  peak  sidelobe  of 
SLLpeak  and  an  asymmetric  sidelobe  pattern  with  sidelobe 
(SLLpeak/SLRpeak)  . Using  the  patterns  plotted  in  Figure  4.6  and 
the  results  listed  in  Table  4.3,  let's  compare  the  (35/23) 
asymmetric  sidelobe  pattern  to  the  Taylor  pattern  with  peak 
sidelobe  35dB.  It  is  obvious  that  the  asymmetric  sidelobe 
radiation  pattern,  as  expected,  has  a substantial  amount  of 
advantages  in  all  characteristic  indices  over  the  symmetric 
Taylor  pattern.  The  incremental  beamwidth  resolutions, 
directivity,  and  excitation  (radiation)  efficiency  of  the 
asymmetric  sidelobe  radiation  pattern  can  be  attributed  to 
having  let  the  sidelobes  on  the  right  side  rise  from  the  level 

SLLpeak  to  the  level  SLRpeak* 

The  consequences  of  the  design  examples  show  that  the 
proposed  design  method  does  not  allow  a strong  control  on  the 
right-hand  side  sidelobe  (SLRpeak)  of  an  asymmetric  sidelobe 
pattern  due  to  the  convolution  process.  However,  the  result 
is  close  enough  to  achieve  the  most  practical  utilization  of 
the  asymmetric  sidelobe  radiation  pattern,  and  it  is 
reasonable  to  have  a conclusion  that  the  design  methodology 
using  a modulated  analytic  excitation  provides  an  effective 
and  straightforward  method  to  produce  a desirable  broadside 
narrow  beamwidth  asymmetric  sidelobe  radiation  pattern.  As 
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long  as  a practical  asymmetric  radiation  pattern  is  concerned, 

the  design  procedures  described  above  are  generalized.  For 

instance,  those  procedures  can  be  applied  to  a Dolph-Chebyshev 

array  system  and  the  same  results  can  be  obtained.  In  the 

following  sections,  we  extend  this  design  concept  in  the 

linear  array  system  to  the  2-D  planar  array  system. 

4 . 7 The  Generation  of  Narrow  Beamwidth  Asymmetric 
Sidelobe  Planar  Array  Radiation  Pattern 

The  2-D  analytic  signals  whose  spectrums  are  either  half- 
plane or  quadrant-plane  and  their  associated  Hilbert 
transformers  have  been  defined  and  derived  in  Chapter  II. 
Once,  the  key  transformation  tool — Hilbert  transformer  has 
been  obtained,  the  general  theory,  system  analysis  and  design 
procedures  involved  in  the  generation  of  an  asymmetric 
sidelobe  planar  array  pattern  using  the  concept  of  the 
modulated  analytic  excitation  are  theoretically  the  same  as 
those  applied  in  the  linear  array  system,  because  the 
properties  of  a realizable  2-D  Hilbert  transformer  and  the 
modulation  theory  remain  the  same  in  the  2-D  case. 
Nevertheless,  we  want  to  show  that  this  design  methodology  can 
be  extended  to  the  plannar  array  system  as  well.  However,  we 
will  omit  the  details  of  the  fundamental  theory  and  system 
analysis,8  and  directly  get  into  the  generation  of  the  planar 
asymmetric  sidelobe  patterns. 


Since  they  are  the  same  as  the  linear  case,  except  some 
symbolic  manipulation. 
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4-7.1  The  generation  of  half-plane  asymmetric  sidelobe 

radiation  pattern  using  Type-II  analytic  excitation 

Let's  first  consider  the  case  that  we  want  to  generate 
a half-plane  asymmetric  sidelobe  radiation  pattern  in  the  u 
plane.  Thus,  let's  assume  that  we  are  given  a 2-D  real 
excitation  Jr(m,n)  whose  pattern  is  denoted  as  Fr(u,v) . Using 
the  same  concept  in  the  linear  array  system,  in  order  to 
minimize  the  distortion  caused  by  the  realizable  2-D  Hilbert 
transformer  (Type  II)  around  the  line  of  u=0,  and  to  eliminate 
the  half-plane  sidelobe,  we  can  generate  a 2-D  modulated  real 
excitation  as 

Jmr  (m, n)  = Jr(m,n)sin(ucn)  (4.21) 

where  uc  is  the  angle (in  radians)  of  the  half  mainbeam  width 
of  the  Fr(u,v)  at  u plane.  Thus,  applying  the  modulation 
property  of  the  Fourier  transform,  the  pattern  corresponding 
to  the  excitation  in  Equation  (4.21),  denoted  as  Fmr(u,v),  is 

Fn,r(u'v)  = (j/2)  [Fr(U+Uc,V)  - Fr(u-Uc,v)]  (4.22) 

That  is,  the  original  pattern  Fr(u,v)  has  been  shifted  off  the 
line  of  u=0  by  the  modulation  process.  Then,  by  using  the 
Type-II  Hilbert  transformer  and  a suitable  window  function, 
we  can  obtain  a 2-D  modulated  analytic  excitation  Jma(rn,n)  as 

Jma(m'n)  = Jmr(m'n)  + j ( [ Jmr  (m  / n)  **hn  (m , n)  ] W (m , n)  } (4.23) 
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where  **  indicates  the  2-D  discrete-time  convolution,  hn(m,n) 
is  the  impulse  response  of  the  type  II  Hilbert  transformer, 
and  w(m,n)  is  the  2-D  truncation  window.  Note  that  because 
hn(m,n)  as  given  by  Equation  (2.55)  is  equivalent  to  h(m) , the 
convolution  in  Equation  (4.23)  can  be  carried  out  as  a 1-D 
convolution  and  w(m,n)  becomes  w(m) . Thus,  the  best  choice  of 
the  truncation  window  w(m)  should  be  the  same  as  that 
discussed  in  the  linear  array  system,  that  is,  it  is  the 
Kaiser  window  with  6=0.  Accordingly,  the  resulting  radiation 
pattern,  denoted  as  Fma(u,v),  corresponding  to  the  2-D 
modulated  analytic  excitation  Jma(m,n)  in  Equation  (4.23) 
should  have  a suppressed  sidelobe  in  u<0  half-plane,  and 
should  gain  various  advantages  compared  to  the  conventional 
symmetric  planar  array  pattern.  Let's  first  use  a nonseparable 
2-D  circular  Kaiser  window  function  as  a symmetric  planar 
array  excitation  to  illustrate  the  generation  of  a desired 
half-plane  asymmetric  sidelobe  array  pattern  using  the 
methodology  discussed  above.9  The  equivalent  2-D  nonseparable 
circular  Kaiser  excitation  can  be  expressed  as  [61] 


Jr(m,n)  = wk(m,n)  = 


I0[a(l-[  (m2+n2)  /25  ] ) 1/2 ] 
Io2[“] 

0, 


-,  m2+n2<25 


otherwise. 

(4.24) 


9 • 

This  example  is  for  demonstration  only,  due  to  the  very 
poor  excitation  (radiation)  efficiency  of  the  heavily  tapered 
Kaiser  window. 
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where  I0[]  is  the  modified  Bessel  function  of  the  first  kind 
of  order  zero,  and  a is  the  parameter  used  to  adjust  the 
sidelobe  level  and  the  mainbeam  width.  With  a=4  and  m=n=5,  the 
symmetric  planar  radiation  pattern  is  plotted  in  Figure 
4.7(a),  where  the  sidelobe  level  is  about  30dB.  The  contour 
of  the  3dB  beamwidth  (HPBW)  of  the  symmetric  Kaiser  planar 
pattern  is  ploted  in  Figure  4.7(b) . The  corresponding  realized 
half-plane  asymmetric  sidelobe  pattern  Fma(u,v)  and  its  3dB 
beamwidth  (HPBW)  contour  are  plotted  in  Figure  4.7(c)  and  (d) 
respectively.  As  seen  from  Figure  4.7(c),  the  sidelobe  of  the 
planar  pattern  Fma(u,v)  in  the  u<0  half-plane  has  been  greatly 
suppressed,  and  its  peak  sidelobe  levels  become  38dB  and  27dB 
in  the  u<0  and  u>0  half-planes  respectively.  However,  by 
comparing  the  contour  plots  in  Figure  4.6(b)  and  (d) , we  find 
that  the  3dB  beamwidth  (HPBW)  of  the  realized  asymmetric 
sidelobe  planar  pattern  remains  the  same  as  (at  least  very 
close  to)  that  of  the  original  symmetric  planar  array  pattern, 
except  the  shifted  mainbeam  maximum  due  to  the  modulation 
process . 

Next,  by  applying  the  principle  of  separability  of  the 
planar  array  excitation  in  Equation  (3.14),  we  can  form  a 
symmetric  planar  pattern  as 

Fr(u,v)  = [Fr(u)Fr(v)]  (4.25) 

where  Fr(u)  and  Fr(v)  are  assumed  to  be  the  Taylor  pattern 
with  n=4  and  27dB  sidelobe  level.  Let  the  number  of  the  array 
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elements  be  (11x11) . The  symmetric  Taylor  planar  pattern 
generated  by  using  Equation  (4.25)  and  its  3dB  beamwidth 
contour  plot  are  plotted  in  Figure  4.8(a)  and  (b)  , and  the 
corresponding  asymmetric  sidelobe  planar  array  pattern  and  its 
3dB  beamwidth  (HPBW)  contour  plot  are  plotted  in  Figure  4.8(c) 
and  (d) . 

By  comparing  Figure  4.8(a),  (b)  to  Figure  4.8(c),  (d) 
respectively,  we  can  find  the  same  results  as  shown  in 
previous  Kaiser  circular  planar  array  patterns.  That  is,  the 
planar  array  pattern  excited  by  the  corresponding  modulated 
analytic  excitation  has  a suppressed  sidelobe  in  the  u<0  half- 
plane and  has  about  the  same  3dB  beamwidth  as  the  original 
symmetric  Taylor  planar  pattern.  It  is  interesting  to  note 
that  the  asymmetric  sidelobe  levels  which  appear  in  Figure 
4.8(c)  are  34.5dB  and  22.5dB  in  the  u<0  and  u>0  half-planes, 
which  is  close  enough  to  the  results  of  the  linear  array 
system  (refer  to  Example  2 in  Section  4.6).  Thus,  it  is 
reasonable  to  conclude  that  the  analysis  results  and  the 
design  procedures  of  an  asymmetric  sidelobe  pattern  in  the 
linear  array  system  can  be  applied  to  the  half-plane 
asymmetric  sidelobe  planar  array  pattern,  provided  that  the 
excitation  for  a planar  array  system  is  separable. 

4.7.2  The  generation  of  cruadrant-plane  asymmetric  sidelobe 
radiation  pattern  using  Type  I analytic  excitation 

There  are  two  kinds  of  Type  I 2-D  analytic  excitations 
which  correspond  to  the  analytic  signals  of  Case  la  and  lb  in 
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Chapter  II  respectively.  The  generation  of  the  Type  I analytic 
excitation  of  Case  la  involved  two  stages  of  transformations 
(refer  to  Figure  2.6).  However,  as  far  as  the  excitation 
(radiation)  efficiency  of  an  array  excitation  is  concerned, 
we  found  that  the  analytic  excitation  of  Case  la  is  physically 
inapplicable  in  the  realization  of  the  quadrant -plane 
asymmetric  sidelobe  planar  array  radiation  patterns  since  the 
excitation  at  the  output  of  the  first  transformation  in  Figure 
2.6  is  formed  by  directly  substracting  two  real  excitations 
(not  quadratually)  that  will  greatly  reduce  the  efficiency  of 
the  realized  array  excitation.  However,  by  using  the  Case  lb 
analytic  excitation  which  is  equivalent  to  applying  the 
pattern  multiplicative  theory  of  the  array  antenna  system,  a 
quadrant -plane  asymmetric  sidelobe  planar  array  radiation 
pattern  can  be  easily  realized,  and  consequently,  the 
generation  of  a quadrant-plane  asymmetric  sidelobe  planar 
pattern  is  nothing  more  than  the  generation  of  two  asymmetric 
sidelobe  linear  array  patterns.  For  purpose  of  the  conceptual 
simplicity,  we  will  not  further  simulate  the  quadrant-plane 
asymmetric  sidelobe  planar  array  pattern,  as  long  as  the 
pattern  multiplicative  concept  is  applied. 

4.8.  Discussion  and  Summary 

There  are  two  major  contributions  in  this  chapter.  First, 
the  potential  advantages  of  the  asymmetric  sidelobe  array 
patterns  compared  to  the  conventional  symmetric  array  patterns 
have  been  quantitatively  shown.  Second,  an  effective  and 
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analytic  method  has  been  introduced  to  the  design  of  the 
desired  broadside  narrow  beamwidth  asymmetric  sidelobe 
radiation  patterns.  It  is  the  philosophy  of  associating  the 
modulation  theory  and  the  analytic  signal  concept  that  makes 
the  design  of  an  asymmetric  sidelobe  broadside  array  radiation 
pattern  so  straightforward.  Finally,  it  is  found  that  the 
design  concept  and  procedures  for  an  asymmetric  sidelobe 
linear  array  radiation  pattern  can  be  extended  to  the  design 
of  a half-plane  asymmetric  sidelobe  planar  array  radiation 
pattern. 


/r’ad  la  t ton  pa  t tarn  C d B) 
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(b) 

Figure  4.7  The  nonseparable  Kaiser  circular  and  the 

corresponding  half-plane  asymmetric  sidelobe  planar 
array  patterns. 

(a)  The  symmetric  Kaiser  planar  array  pattern. 

(b) The  3dB  beamwidth  contour  plot  of  (a). 

(c) The  half-plane  asymmetric  sidelobe  array  pattern. 

(d) The  3dB  beamwidth  contour  plot  of  (c) . 
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Figure  4.7 


continued 
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(b) 

Figure  4.8  The  11x11  separable  Taylor  and  the  corresponding 
half-plane  asymmetric  sidelobe  planar  array 
patterns . 

(a)  The  symmetric  Taylor  planar  array  pattern. 

(b) The  3dB  beamwidth  contour  plot  of  (a). 

(c) The  half-plane  asymmetric  sidelobe  array  pattern. 

(d) The  3dB  beamwidth  contour  plot  of  (c) . 
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(d) 


Figure  4 . 8 


continued 


CHAPTER  V 


THE  SHAPED  BEAM  ARRAY  ANTENNA  RADIATION  PATTERN 
DESIGN  USING  ANALYTIC  EXCITATION 

5 . 1 Introduction 

By  extending  the  design  concept  and  technology  in  the 
preceeding  chapter,  we  proceed  to  design  the  shaped  beam 
radiation  patterns  of  an  array  antenna  system  using  an 
analytic  excitation. 

In  addition  to  the  narrow  beamwidth  radiation  pattern, 
the  design  of  the  wide  angle  shaped  beam  radiation  pattern  is 
another  important  subject  in  the  area  of  antenna  system.  For 
instance,  a cosecant-squared  (esc 2d)  power  radiation  pattern 
is  widely  applied  in  the  ground-mapping  and  air-search  radar 
systems.  The  most  popular  methods  used  in  realization  of  a 
shaped  beam  pattern  were  the  Fourier  series  expansion  and  the 
Woodward's  voltage  pattern  sampling  method  [7,21,35]  which 
result  in  a minimum  mean-square-error  (MMSE)  approximation  to 
the  desired  shaped  beam  pattern.  Several  design  methods  for 
improving  the  shape  of  either  a voltage  or  a power  shaped  beam 
pattern  using  iterative  numerical  methods  have  been  found  in 
the  literature  [37-42],  Since  the  polynomial  representation 
of  an  asymmetric  shaped  beam  pattern  is  normally  complex,  the 
design  of  a shaped  beam  pattern  using  numerical  iterative 


131 


132 


method  is  not  as  effective  and  straightforward  as  the 
symmetric  pattern  which  can  always  be  represented  as  a real 
polynomial.  For  a large  array  system,  the  procedures  become 
very  ineffective,  cumbersome  and  sometimes  impossible. 
Moreover,  the  usual  wide  shaped  beamwidth  further  exacerbates 
the  fundamental  difficulties  in  designing  the  shaped  beam 
pattern  using  iterative  numerical  method,  and  skill  and  luck 
are  often  required  in  the  numerical  manipulation. 

In  this  chapter,  a totally  new  method  is  introduced  for 
designing  a typical  shaped  beam  pattern  using  an  analytic 
excitation.  It  will  be  shown  that  this  design  method  has 
greatly  alleviated  the  computational  effort,  and  has  provided 
an  idea  to  effectively  apply  the  well-developed  equiripple  (or 
minimax)  design  criterion  in  the  realization  of  an  asymmetric 
shaped  beam  pattern.  Furthermore,  the  inherent  distortion 
which  results  from  a realizable  Hilbert  transformer  can  be 
minimized,  and  a desirable  asymmetric  sidelobe  shaped  patttern 
can  be  easily  generated.  Thus,  this  design  methodology  is 
increasingly  attractive  from  the  points  of  view  of  practical 
utilization  of  the  shaped  beam  pattern  of  an  array  antenna 
system. 

5 . 2 The  Basic  Concepts  for  Shaped  Beam  Radiation  Pattern 
Design  using  Analytic  Excitation 

A desired  shaped  beam  radiation  pattern  is  always 
ill-defined,  and  the  design  of  a shaped  beam  pattern  becomes 
an  approximation  problem.  The  concepts  of  designing  a shaped 
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beam  radiation  pattern  using  an  analytic  excitation  are  first 
realizing  a symmetrical i zed  desired  shaped  beam  pattern, 
secondly  applying  the  analytic  signal  theory  to  eliminate  the 
half  side  of  the  symmetrical ized  shaped  beam  pattern,  and 
thirdly,  following  the  phase  shifting  property  of  the  array 
antenna  system  to  shift  the  realized  shaped  beam  pattern  back 
to  the  reguired  angular  position. 

Assume  a typical  desired  shaped  beam  pattern  of  an  N=2M+1 
elements  array  system,  denoted  as  Fd(u),  with  shaped  beamwidth 
us  depicted  as  in  Figure  5.1(a),  where  without  loss  of  any 
generality,  Fd(u)  is  assumed  to  be  real,  and  its  starting 
angle  of  the  shaped  region  is  assumed  to  be  exactly  positioned 
at  u— 0 . Then,  the  symmetricalized  shaped  beam  pattern,  denoted 
as  Fds(u),  can  be  expressed  as 


Fds(u) 


Fd(u"uc)  ' 

0, 

Fd(-U-Uc)  ' 


0 < u < 7 r 
u=0 

-n  < u < 0 


(5.1) 


where  the  value  of  uc  is  reguired  to  allow  a roll-off  region1 
which  occurs  in  the  realization  of  Fds(u)  and  to  eliminate  the 
distortion  caused  by  the  realizable  Hilbert  transformer  (will 
be  discussed  in  next  section) . The  symmetricalized  shaped  beam 
pattern  as  given  by  Equation  (5.1)  can  be  depicted  as  in 
Figure  5. 1 (b) , 


The  roll-off  region  here  is  defined  as  the  angular  width 
between  the  first  null  and  the  desired  magnitude  of  the 
realized  pattern  (or  spectrum) . 
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Figure  5.1  Sketches  of  the  desired  shaped  beam,  the 

symmetrical ized  shaped  beam  pattern  and  the 
realized  shaped  beam  pattern. 

(a) The  ideal  desired  shaped  beam  pattern  Fd(u) . 

(b) The  symmetricalized  shaped  beam  pattern  Fds(u) . 

(c) The  realized  symmetricalized  shaped  beam 
pattern  Fr(u). 

(d) The  realized  shaped  beam  pattern  Fta(u+uc). 
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After  the  symmetrical ization  process  and  according  to  the 
property  of  the  discrete-time  Fourier  transform,  the 
excitation  for  the  symmetrical ized  shaped  beam  pattern, 
denoted  as  Jds(n),  should  be  real  and  symmetric.  Thus,  the 
array  shaped  beam  pattern  Fds(u)  in  Equation  (5.1)  can  be 
expressed  as 

00 

Fds(u)  = £ a (n)  cos  (un)  (5.2) 

n=0 

where  a(0)=Jds(0)  and  a(n)=2Jds(n)  for  n^O.  Note  that  a desired 
shaped  beam  pattern  may  only  be  realized  with  an  infinite 
number  of  array  elements.  From  the  mathematical  representation 
of  the  symmetricalized  shaped  beam  pattern  as  given  by 
Equation  (5.2),  it  is  easy  to  see  that  the  symmetricalization 
process  not  only  converts  the  complex  polynomial 
representation  of  an  asymmetric  shaped  beam  pattern  into  a 
real  polynomial  representation,  but  also  allows  us  to  consider 
only  half  of  the  array  excitations  and  pattern,  and  thus 
greatly  simplifies  the  design  complexity  and  reduces  the 
computational  effort.  Furthermore,  because  of  the  form  of 
(5.2),  the  well-developed  equiripple  (minmax)  design  criterion 
[1,13,56-57]  can  then  be  applied  to  the  design  of  the 
symmetricalized  shaped  beam  pattern  Fds(u). 

It  is  the  analytic  signal  theory  that  makes  the  above 
symmetricalization  process  possible.  Let's  assume  that  a 
design  method  has  been  applied  to  realize  the  ideal 
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symmetricalized  shaped  beam  pattern  Fds(u),  and  let  the 
realized  symmetricalized  shaped  beam  pattern  and  its 
corresponding  excitation  be  represented  as  Fr(u)  and  Jr(n) 
respectively.  Note  that  Jr(n)  becomes  an  N elements  array 
excitation.  Let  the  realized  symmetricalized  shaped  beam 
pattern  Fr(u)  be  depicted  as  Figure  5.1(c).  Then,  by  applying 
the  analytic  signal  theory,  we  can  obtain  a realizable 
analytic  excitation  such  that  to  eliminate  one  of  the  shaped 
beams  of  Fr(u)  presented  within  the  visible  range.  That  is, 
the  analytic  excitation,  denoted  as  Jta(n),  can  be  formed  as 

Jta(n)  = Jr(n)  + jJti(n)  (5.3) 

where  Jti(n)  denotes  a truncated  N elements  array  excitation 
generated  from  Jr(n)  by  using  the  realizable  discrete-time 
Hilbert  transformer.  The  radiation  pattern  Fta(u)  corresponding 
to  the  analytic  excitation  of  Jta(n)  in  Eguation  (5.3)  is 

Fta(u)  = Mu)  + jFti  (u)  (5.4) 

where  Ft1(u)  is  the  pattern  corresponding  to  Jti(n). 
Accordingly,  Fta(u)  represents  the  realized  pattern  which  has 
single  shaped  beam  remaining  within  the  visible  range  of  an 
array  antenna  system.  Because  of  the  symmetricalization 
process,  the  realized  shaped  beam  pattern  Fta(u)  is  shifted 
with  angle  uc.  Then,  using  the  shifting  property  of  an  array 
antenna  system  with  a linear  phase  gradient  which  is  the  value 
of  Uc'  Fta(  u)  can  be  shifted  back  to  the  original  required 
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angular  position.  Thus,  the  ultimate  realized  shaped  beam 
pattern,  denoted  as  Fta(u+uc),  is  depicted  as  in  Figure  5.1(d). 
Since  the  symmetrical ization  process  is  a linear  operation, 
the  shape  of  the  pattern  in  the  to-be  shaped  region  remains 
the  same  as  the  one  which  is  synthesized  without  shifting 
process.  However,  based  on  the  property  of  an  analytic 
excitation,  at  least,  parts  of  half-side  sidelobe  of  the 
realized  shaped  beam  pattern  can  be  eliminated  (at  least 
suppressed) . As  discussed  in  Chapter  IV,  the  phenomenon  of  the 
eliminated  sidelobe  is  always  a desirable  property  in  the 
practical  application  of  an  antenna  system.  This  is  an 
important  advantage  obtained  from  this  design  methodology. 

5 . 3 System  Design  Constraint 

In  order  to  minimize  the  distortion  caused  by  a 
realizable  Hilbert  transformer  in  the  physical  implementation 
of  the  shaped  beam  pattern  using  analytic  excitation,  the 
quantity  of  uc  in  Equation  (5.1)  is  required  to  shift  the 
realized  symmetrical ized  shaped  beam  pattern  Fr(u)  to  the 
pass-band  of  the  realizable  Hilbert  transformer.  Let's 
delineate  the  roll-off  regions  (at  the  u=0  end)  of  both 
realized  symmetricalized  shaped  beam  pattern  Fr(u)  and 
spectrum  H(u)  of  the  Hilbert  transformer  as  in  Figure  5.2, 
where  Auh  denotes  the  roll-off  region  of  H(u) , Aur  represents 
the  roll-off  region  of  Fr(u)  and  Aunull  represents  the  angular 
distance  between  the  first  null  of  the  Fr(u)  and  the  starting 
angle  of  the  shaped  region  of  the  desired  symmetricalized  beam 
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pattern  Fds(u).  It  is  obvious  that  the  value  of  uc  should  at 
least  satisfy  the  following  condition: 

(Aunu11+Auh)  ^ uc  and  (Uc+Us+Aunull+Auh)  < n (5.5) 

and  Fr(u)  should  have  a small  amount  of  component  around  the 
u=0  direction  to  assure  that  the  distortion  caused  by  a 
realizable  Hilbert  transform  can  be  minimized. 

However,  in  the  practical  application,  it  is  advantageous 
to  eliminate  (at  least  suppress)  one  of  two  sides  sidelobes 
of  the  shaped  beam  pattern.  Thus,  the  value  of  uc  is  the 
possible  value  of  Aunu11.2 


Figure  5.2  The  definitions  of  the  roll-off  regions  of  the 
realized  symmetricalized  shaped  beam  pattern 
and  the  spectrum  of  the  realizable  Hilbert 
transformer. 


It  allows  a negligible  incremental  of  distortion  in  the 
shaped  beam  pattern  caused  by  the  realizable  Hilbert 
transformer. 
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The  extent  of  the  system  constraint  in  Equation  (5.5)  is 
actually  very  much  dependent  on  the  number  of  array  elements. 
The  larger  the  number  of  array  elements,  the  smaller  the 
limitation  on  the  angular  width  of  the  shaped  beamwidth  which 
can  be  realized  by  using  this  proposed  design  method. 

5 • 4 System  Design  Example  using  Windowing  Method 

Due  to  the  conceptual  simplicity,  it  is  easy  to  visualize 
the  design  process  in  the  realization  of  a shaped  beam  pattern 
using  an  analytic  excitation.  It  is  obvious  that  an  equiripple 
criteria  of  numerical  design  method  is  possible  and  preferred 
in  realizing  such  a symmetrical ized  shaped  beam  pattern  as 
given  by  Equation  (5.2).  However,  by  using  the  windowing 
method,  we  are  able  to  analytically  generate  a shaped  beam 
pattern,  and  to  demonstrate  the  simplicity  and  feasibility  of 
this  new  design  methodology.  That  is  the  windowing  method  is 
more  appropriate  for  the  purpose  of  illustration  of  this 
design  method  introduced. 

The  most  popular  shaped  beam  pattern  is  the  cosec 20  power 
pattern.  Since  it  gives  a constant  echo-signal  power  which  is 
independent  from  6 for  a target  moving  at  the  same  altitude, 
it  was  widely  used  in  an  air-search  and  the  ground-mapping 
radar  system.  A typical  requirement  of  a shaped  beam  pattern 
in  the  radar  application  is  that  from  1°  to  5°  is  the  vertical 
main  beam  (that  is  the  beamwidth  in  the  elevation  direction 
is  4°)  and  with  cosecant-squared  shaped  beam  from  5°  to  30°. 
Using  the  linear  array  system  as  discussed  in  Section  3.2.2 
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Chapter  II,  let  the  phase  reference  be  shifted  to  the  array 
center  point  and  the  interelement  distance  to  be  half  of  the 
wavelength  (d=A/2) , the  normalized  desired  cosecant  voltage 
shaped  pattern  is3 


Fd(*> 


is3 

1, 

sin(5°) 

85°  < e 

60°  < e 

< 89 

< 85' 

sin  (9O°-0 ) 

0, 

elsewhere 

(5.6) 


where  sin (5°)  is  used  to  normalize  the  shaped  beam  pattern. 
Note  that  u=7r cos  {6 ) for  half  wavelength  interelement  distance. 
After  changing  the  variable  of  6 to  u,  the  normalized  shaped 
voltage  beam  pattern  becomes 


Fd(u) 


U1  < u < u2 

u2  < u < u3  (5.7) 

elsewhere 


where  U]=7rcos  (89°)  , u0=u2=7rcos  (85°)  , and  u3=7rcos  ( 60°)  . The 
desired  shaped  voltage  beam  pattern  as  described  by  Eguation 
(5.7)  is  plotted  in  Figure  5.3. 


3Note  that  the  angle  6 in  the  Figure  3.2  of  the  linear 
array  system  is  measured  from  the  endfire  direction. 
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u (radian) 


Figure  5.3  The  desired  cosecant-squared  shaped  beam  power 
pattern. 

Having  defined  the  desired  shaped  beam  pattern,  we  can 
symmetrical ize  the  desired  shaped  beam  pattern  by  determining 
the  value  of  uc  in  Equation  (5.1).  According  to  Equation 
(5.5)  , the  value  of  uc  is  actually  determined  by  the  number  of 
array  elements  N as  well  as  the  characteristics  of  the 
windowing  functions  which  are  also  functions  of  the  number  of 
array  elements.  The  number  of  array  elements  determines  the 
degree  of  the  approximation  of  the  realized  pattern  to  the 
desired  shaped  pattern.  There  are  two  window  functions 
involved  in  using  the  windowing  method  to  realize  this  design 
technology.  The  first  one  is  used  to  realize  the  ideal 
symmetrical ized  shaped  beam  pattern,  and  the  other  one  is  used 
to  implement  the  finite  elements  of  the  analytic  excitation. 
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Since  the  mainlobe  widths  of  the  window  functions  are  to  be 
less  than  the  pass  band  of  the  desired  shaped  beam  pattern 
and  the  spectrum  of  the  ideal  Hilbert  transformer,  based  on 
the  convolution  process,  a fairly  accurate  qualitative 
estimation  of  the  angular  width  of  Aunull  and  the  roll-off 
region  Auh  is  possible.  That  is,  the  angular  width  Aunull  which 
is  the  angular  width  between  the  first  null  point  in  the 
realized  symmetricalized  shaped  beam  pattern  and  the 
discontinuity  in  the  ideal  symmetricalized  shaped  beam  pattern 
is  to  be  half  of  the  mainlobe  width  of  the  first  window 
function,4  and  the  roll-off  region  Auh  at  u=0  end  of  the 
realized  Hilbert  transformer  is  to  be  the  half  of  the  mainlobe 
of  the  second  window  function.  The  first  window  function  will 
determine  the  angular  width  Aunull  and  the  level  of  the 
sidelobe  in  the  realized  symmetricalized  shaped  pattern.  In 
order  to  have  a narrower  roll-off  region  (it  is  a desirable 
requirement  in  designing  a shaped  beam  pattern)  the  first 
window  should  be  Dolph-Chebyshev  window.  However,  in  order  to 
possibly  suppress  one  of  two  sides  sidelobes,  the  best  choice 
for  the  second  window  becomes  the  Kaiser  window  with  6=0.  The 
reasons  to  choose  such  window  functions  are  actually  the  same 
as  those  stated  in  the  preceeding  chapters. 

After  determining  the  type  of  the  windows,  let's  assume 
that  the  number  of  array  elements  is  17  [40],  and  a sidelobe 

4The  value  of  uc  is  assumed  to  be  larger  or  equal  to  the 
half  of  the  mainlobe  width  of  the  first  window  function. 
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ratio  R=15dB  Dolph-Chebyshev  window  is  used  to  realize  the 
symmetricalized  shaped  beam  pattern.5  Then,  the  values  of  Auh 
and  Aunul1  can  be  estimated  from  the  mainlobe  widths  of  the 
chosen  window  functions,  that  is  Auh=(27r/17)  which  is  equal 
to  half  of  the  mainlobe  width  of  a uniform  window  function  and 
Aunul1=  (1 . 9^/17)  which  is  equal  to  half  of  the  mainlobe  width 
of  the  Dolph-Chebyshev  window  with  R=15dB.  Accordingly,  the 
value  of  uc=0. 7 (radian)  can  satisfy  the  constraint  in  Equation 
(5.5),  and  the  symmetricalized  cosecant  shaped  voltage  beam 
pattern  can  be  represented  as 


Fds(u) 


1» 

u„ 


(U-Uc+Uj) 

0, 


uc  < I U I < (U^Ug-Ui) 

(Uc+U2"Ul)<lUl^  (Uc+U3“Ul) 

elsewhere 


(5.8) 


According  to  Equation  (5.2),  the  corresponding  excitation  for 
Fds(u)  in  Equation  (5.8)  can  be  obtained  by  using  the  inverse 
Fourier  transform.  After  doing  the  Inverse  transformation  and 
truncation  process,  the  excitation  for  the  realized 
symmetricalized  shaped  beam  voltage  pattern  can  be  expressed 
as 


Jr(n)  = a(n)w1(n) 


-M  < n < M 


(5.9) 


5The  mainlobe  width  of  the  Dolph-Chebyshev  window  can  be 
obtained  from  Figure  A.l  in  Appendix,  where  the  mainlobe  width 
is  equal  to  the  mainbeam  width  ratio  times  (4tt/N)  . 
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where  M=8,  w^n)  is  the  Dolph-Chebyshev  window  with  sidelobe 
ratio  R=15dB,  and  a(n)  is  the  excitation  for  the  ideal 
symmetricalized  shaped  beam  pattern  Fds(u)  as  given  by  Equation 
(5.8),  which  is  calculated  as 


a (n) 


1 


7T 


n 

0 


Fds  (u)  cos  (un)  du 


(5.10) 


It  can  be  seen  that  because  of  the  symmetric  property  of 
Fds(u),  the  calculation  of  the  excitation  is  simplified. 

The  resulting  symmetricalized  shaped  beam  pattern 
corresponding  to  the  excitation  Jr(n)  is  plotted  in  Figure 
5.4(a).  Then  applying  the  Hilbert  transform  and  truncation 
process,  the  analytic  excitation  Ja(n)  corresponding  to  the 
real  excitation  Jr(n)  in  Equation  (5.9)  can  be  obtained  as 

Ja(n)  = Jr(n)  + j [Jr(n)  *h(n) ]w2(n)  (5.11) 

where  * indicated  the  discrete-time  convolution,  h(n)  is  the 
impulse  response  of  the  discrete-time  Hilbert  transformer,  and 
w2(n)  is  the  uniform  window  function.  The  realized  cosecant 
shaped  beam  pattern  corresponding  to  the  analytic  excitation 
Ja(n)  in  Equation  (5.11)  after  shifting  back  to  the  original 
required  position  is  plotted  in  Figure  5.4(b).  The  realized 
pattern  in  Figure  5.4(b)  shows  that  the  levels  of  the  sidelobe 
in  the  non-shaped  region  and  the  ripple  in  the  shaped  region 
are  high.  By  changing  the  sidelobe  ratio  of  the  Dolph- 
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Chebyshev  window  function  in  the  realization  of  the 
symmetrical ized  shaped  pattern  in  Equation  (5.8),  we  can 
control  the  levels  of  the  sidelobe  and  the  pass  band  ripple. 
Thus,  by  applying  the  same  procedures  and  using  the  Dolph- 
Chebyshev  window  with  R=20dB  in  realization  of  the 
symmetricalized  shaped  beam  pattern,  the  corresponding 
realized  cosecant  shaped  beam  pattern  is  plotted  in  Figure 
5.4(c),  where  the  value  of  uc  is  set  to  0.5(radian).  When  we 
use  a Dolph-Chebyshev  window  with  R=30dB  to  implement  the 
symmetricalized  shaped  voltage  pattern,  and  adjust  the  first 
null  point  of  the  symmetricalized  pattern  to  be  at  u=0,  where 
the  value  of  uc  can  easily  be  estimated  to  be  uc=0 . 54  (radian) 
(uc=Aunu-|1=half  of  the  mainlobe  width  of  the  Dolph-Chebyshev 
window  with  R=30dB) , the  resulting  shaped  cosecant-squared 
power  pattern  is  plotted  in  Figure  5.4(d).  Comparing  those 
patterns  in  Figure  5.4(b),  5.4(c)  and  5.4(d),  it  is  easy  to 
see  that  the  lower  the  sidelobe  of  the  Dolph-Chebyshev 
truncation  window  used  to  realize  the  symmetricalized  shaped 
beam  pattern,  the  lower  the  sidelobe  level  and  pass  band 
ripple,  but,  the  broader  the  shaped  region  in  the  realized 
cosecant-squared  power  pattern.  The  pattern  plotted  in  Figure 
5.5  illustrates  the  improvement  of  the  approximation  by 
increasing  the  size  of  the  array  elements  from  17  to  51,  and 
the  window  function  in  realization  of  the  symmetricalized 
shaped  beam  pattern  is  the  Dolph-Chebyshev  with  R=30dB. 
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u (radian) 

(a) 


u (radian) 

□ ideal  pattern  + rea l i zed  pattern 


(b) 

Figure  5.4  The  desired  and  realized  cosecant-squared  shaped 
power  patterns. 

(a)  The  realized  symmetrical i zed  shaped  beam  pattern. 

(b)  (c) (d)  are  the  realized  cosecant-squared  power 
patterns  corresponding  to  R=15dB,  20dB  and  35dB 
of  Dolph-Chebyshev  truncation  windows. 
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u (radian) 

□ ideal  pattern  + real i zed  pattern 


(c) 


u (radian) 

□ ideal  pattern  + rea I i zed  pattern 


(d) 


Figure  5 . 4 


continued 
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Figure  5.5  The  realized  cosecant-squared  shaped  power 

pattern  with  increasing  number  of  array  elements. 
(N=51) . 

The  excitations  for  the  shaped  beam  cosecant-squared 
power  pattern  in  Figure  5.4(b),  5.4(c)  and  5.4(d)  are 

tabulated  in  Table  5.1,  where  the  desirable  complex  conjugate 
symmetric  property  of  an  array  excitation  is  maintained. 
Actually  the  phase  in  each  array  element  should  be  the 
summation  of  the  tabulated  value  and  its  corresponding  linear 
phase  gradient  which  is  n times  uc,  where  uc  is  the  value  used 
to  symmetricalize  the  shaped  pattern  in  each  realization  case. 
Thus  by  adding  the  linear  phase  nuc  to  each  corresponding 
number  of  element,  we  are  able  to  shift  the  realized  shaped 
patterns  back  to  the  original  required  angular  position  shown 
in  Figure  5.4(b),  5.4(c)  and  5.4(d).  The  characteristics  of 


149 


these  three  realized  cosecant-squared  shaped  beam  power 
patterns  are  also  listed  in  Table  5.2,  where  the  parameter  of 
Jmax/Jmin  t^ie  ratio  of  the  maximum  to  the  minimum  amplitude 
in  the  element  excitations,  which  is  related  to  the  excitation 
(radiation)  efficiency  of  the  array  antenna  system. 


Table  5.1  The  excitations  of  the  realized  cosecant-squared 
power  patterns.  N=17,  J(-n)=J*(n) 


array 

number 

Figure 

5.4(b) 

Figure 

5.4 (c) 

Figure 

5.4(d) 

Amplitude  phase 

amplitude  phase 

amplitude  phase 

0 

.222 

. 000 

.222 

. 000 

.222 

.000 

1 

.202 

68 .486 

.201 

57.078 

.200 

59 . 399 

2 

. 153 

-46.509 

. 151  - 

■69.535 

. 146 

-64.910 

3 

.098 

9.520 

. 094  - 

■26.097 

. 087 

-19.561 

4 

.067 

44.569 

. 063 

-2.579 

. 055 

7.286 

5 

. 068 

83.714 

.057 

27.748 

.045 

37.797 

6 

.054 

-45.359 

.044 

64.942 

.032 

79.734 

7 

.042 

-25.856 

.039  - 

76.331 

.021 

-59.851 

8 

.055 

18.775 

.031 

-8.132 

.012 

-3.971 

Table  5.2  The  characteristics  of  the  realized 


cosecant-squared  power  patterns. 


type  of 
pattern 

(None) 

SLLpeak 

(dB) 

SLRpeak 

(dB) 

Figure  5.4(b) 

5.29 

-18.12 

-21.48 

Figure  5.4(c) 

7.16 

-36.57 

-26.38 

Figure  5.4(d) 

18.50 

-43.87 

-35.16 

The  results  in  Table  5.2  show  that,  as  expected,  the 
lower  the  sidelobe,  the  poorer  the  excitation  efficiency  of 
the  realized  analytic  excitation.  Furthermore,  as  shown  in 
Figure  5.4(c)  and  5.4(d),  the  sidelobe  levels  below  the 
horizontal  (SLLpeak)  of  the  resulting  shaped  beam  patterns 


150 


using  analytic  excitations  can  be  heavily  suppressed.  In  the 
practical  environmental  situation  where  the  clutters  which 
come  from  ground  or  sea-surface  are  more  severe  than  those 
from  the  space,  the  level  of  the  lower  side  sidelobe  (below 
the  horizon)  is  desired  to  be  lower  than  that  of  upper  side, 
which  implies  that  the  asymmetric  sidelobe  levels  shown  in 
Figure  5.4(c)  and  5.4(d)  may  satisfy  the  sidelobe  requirement 
of  an  antenna  system.  Then,  by  using  the  asymmetric  sidelobe 
shaped  beam  pattern,  we  will  be  able  to  gain  the  advantages 
such  as  excitation  efficiency,  narrower  roll-off  region 
compared  to  the  symmetric  sidelobe  shaped  beam  pattern. 
Actually,  this  is  one  of  the  common  and  attractive  properties 
of  an  analytic  excitation. 

5 . 5 Discussion  and  Summary 

A new  design  technology  for  synthesizing  a shaped  beam 
array  radiation  pattern  has  been  introduced  in  this  chapter. 
From  the  example,  we  can  reach  a reasonable  conclusion  that 
the  new  design  method  provides  a great  reduction  in  the 
computational  effort  compared  to  the  Fourier  expansion  and  the 
iterative  numerical  method,  and  provides  a simple  skill  to 
control  the  sidelobe  level  compared  to  the  Woodward's  voltage 
pattern  sampling  method.  The  most  important  idea  and 
contribution  conveyed  by  this  method  is  that  the  well- 
developed  equiripple  design  criterion  can  be  equally  applied 
to  the  design  of  any  desired  shaped  array  antenna  radiation 
pattern  provided  that  Equation  (5.5)  is  satisfied. 
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Furthermore,  it  is  worthwhile  to  be  emphasized  that  a 
desirable  asymmetric  sidelobe  shaped  beam  pattern  [27]  can  be 
easily  generated  by  using  this  design  methodology. 


Chapter  VI 


ANALYSIS  AND  DESIGN  OF  RADIATION  PATTERNS  FOR 
A MONOPULSE  ARRAY  ANTENNA  SYSTEM  USING 
ANALYTIC  SIGNAL  THEORY 

6 . 1 Introduction 

A monopulse  tracking  radar  system  usually  requires  two 
kinds  of  patterns,  a sum  pattern  and  a difference  pattern. 
The  sum  pattern  is  responsible  to  detect  and  measure  the  range 
of  the  target,  and  the  difference  pattern  is  used  to  acquire 
the  angular  information  of  a moving  target.  Traditionally,  the 
sum  and  difference  patterns  of  an  amplitude  monopulse  array 
antenna  radar  system  are  formed  as  the  summation  and  the 
difference  of  two  squint  patterns  [5],  as  shown  in  Figure  6.1, 
where  us  in  radians  is  usually  called  the  squant  angle,  and 
E(u)  and  A(u)  denote  the  sum  and  difference  patterns 
respectively.  However,  as  the  modern  technology  continues  to 
advance,  it  becomes  increasing  attractive  to  produce  the  sum 
and  difference  patterns  using  dual  excitations  in  an  array 
antenna  system  [67-68].  Actually,  this  separation  concept  has 
already  been  widely  adopted  in  dealing  with  the  problems 
related  to  the  difference  patterns  of  a monopulse  tracking 
radar  system  [5,69-75]. 
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S(u)  A(u) 


Figure  6.1  Sketches  of  the  typical  radiation  patterns  of 
an  amplitude  monopulse  array  antenna  radar 
system. 

(a)  The  overlapping  of  two  squint  patterns. 

(b) The  sum  pattern  E(u). 

(c) The  difference  pattern  A (u) . 
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Theoretically,  any  odd  symmetric  excitation  could  be  used 

to  generate  a desired  difference  pattern.1  It  has  been  shown 

in  Chapter  III  that  the  Hilbert  transform  of  an  even  symmetric 

function  is  always  an  odd  symmetric  function.  As  the  matter 

of  fact,  we  will  show  that  the  real  and  the  imaginary  part  of 

a theoretical  analytic  signal  represent  an  ideal  excitation 

pair  for  the  monopulse  array  antenna  radar  system.  It  is 

therefore  interesting  and  worthwhile  to  exploit  the  concept 

using  the  analytic  signal  to  generate  the  sum  and  difference 

patterns  for  the  monopulse  array  antenna  radar  system. 

6 . 2 The  Generations  of  Radiation  Patterns  for  a Monopulse 
Array  Antenna  Radar  System  using  Analytic  Signal  Concept 

Before  studying  the  utilization  of  the  analytic  signal 
concept  to  generate  the  radiation  patterns  for  a monopulse 
array  antenna  radar  system,  let's  figure  out  some 
characteristic  indices  used  to  measure  the  performance  of  the 
sum  and  difference  patterns.  The  characteristic  indices  used 
to  describe  the  performance  of  the  sum  pattern  for  a monopulse 
radar  system  are  exactly  the  same  as  those  of  the  narrow 
beamwidth  radaition  pattern  described  in  Chapter  III,  that  is, 
the  beamwidth  resolutions,  directivity,  excitation  efficiency, 
and  sidelobe  level.  For  a difference  pattern,  they  are 
directivity,  difference  slope  at  the  boresight  and  sidelobe 
level.  The  definition  of  the  directivity  of  the  difference 

Generally,  the  difference  pattern  is  an  odd  symmetric 
pattern,  such  that  the  angular  information  of  a moving  target 
can  be  obtained. 
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pattern,  denoted  as  Dd,  is  the  same  as  that  defined  in 
Equation  (3.15),  with  F(0O,0O)  and  F(0,0)  replaced  by  A(0m,0m) 
and  A (0,0),  where  A(0m,0m)  is  the  mainbeam  maximum  of  the 
difference  pattern.  The  difference  slope  at  the  boresight  (u=0 
or  9 = n/2)  of  a difference  pattern,  denoted  as  A'(0),  is 
defined  as  [5] 

dA  (u) 

A'  (0)  = (6.1) 

du  u=0 

where  A(u)  is  the  difference  pattern,  and  u=kdcos0 , 6 is 

measured  from  the  endfire  direction.  A' (0)  is  a very  important 
performance  index  of  a difference  pattern,  since  it  is  related 
to  the  tracking  accuracy  of  the  monopulse  radar  system[64]. 

Since  the  sum  pattern  can  be  any  broadside  narrow 
beamwidth  radiation  pattern  using  this  analytic  excitation 
concept,  it  is  not  discussed  here.  In  the  following,  we  will 
show  that  by  applying  the  concept  of  analytic  signal  (analytic 
excitation) , we  can  generate  a difference  pattern  which 
possesses  many  desirable  properties. 

It  has  been  shown  earlier  that  a theoretical  analytic 
excitation  J (n)  =Jr  (n) +j  Ji  (n)  will  produce  a positive-sided 
radiation  pattern,  and  its  dual  Jd (n)  =Jr  (n) -j  Ji  (n)  has  a 
negative-sided  pattern,  where  both  Jr(n)  and  J^n)  are  real, 
and  J^n)  is  the  Hilbert  transform  of  Jr(n).  Let  Fp(u)  and 
Fn(u)  represent  the  positive-sided  pattern  corresponding  to 
J(n)  and  the  negative-sided  pattern  corresponding  to  Jd(n) , 
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respectively.  Accordingly,  a desired  sum  pattern  Z(u)  can  be 
the  summation  of  Fp(u)  and  FN(u)  ; that  is, 

2(u)  = Fp  (u)  + Fn(u) 

= Z [ (Jr(n)+jJi  (n) ) + (Jr(n)-jJ1(n))]e-^un 
n 

= 2 Z Jr(n)  e“^un  = 2 Fr(u)  (6.2) 

n 

where  Fr(u)  is  the  typical  broadside  narrow  beamwidth 
radiation  pattern  corresponding  to  Jr(n) . The  difference 
pattern  can  be  formed  as  the  difference  of  Fp(u)  and  FN(u); 
that  is, 

A(u)  = Fp  (u)  - Fn(u) 

= 2j  Z Ji  (n)  e~lun 
n 

= 2 j Ff  (u)  (6.3) 

where  F^u)  is  the  radiation  pattern  corresponding  to  J^n). 
Furthermore, 

A(u)  = 2 j [-jsgn(u)  Fr(u)  ] 

= 2Fr  (u)  sgn (u) 

= Z (u)  sgn (u)  (6.4) 

where  Fi  (u)  = [-jFr(u)  sgn(u)  ] is  applied.  The  results  in  Equation 
(6.2)  and  (6.3)  indicate  that  if  an  analytic  excitation  is 
obtainable,  then  the  sum  pattern  is  just  the  pattern  excited 
by  its  real  part,  and  the  difference  pattern  is  the  one 
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excited  by  its  imaginary  part.  From  Equation  (6.4),  the 
difference  pattern  is  just  that  of  sum  pattern  except  inverse 
in  the  region  of  u<0,  and  both  of  them  have  exactly  the  same 
power  pattern.  The  result  shown  in  Equation  (6.4)  represents 
one  of  the  most  desirable  relations  between  the  sum  and 
difference  patterns  for  a monopulse  radar  system.  Furthermore, 
according  to  Equation  (6.1),  the  difference  slope  of  the 
difference  pattern  is 


A'  (0) 


dA  (u) 


du 


u=o 


d[E (u) sgn(u) ] 

[ 

du 


u=0 


Due  to  the  discontinuity  of  the  pattern  [E(u)sgn(u]  at 
u=0,  the  derivative  of  [E(u)sgn(u)]  can  mathematically  be 
expressed  as 


A'  (0) 


[R (u)  + S( u)] 


u=0 


6(0)  =« 


(6.5) 


where  R(u)  represents  the  derivative  function  of  the  pattern 
[E (u) sgn (u) ] over  all  u,  except  u=0.  Thus,  the  difference 
slope  of  the  difference  pattern  excited  by  the  imaginary  part 
of  a theoretical  analytic  excitation  is  infinite. 

As  mentioned  eariler,  the  real  part  of  an  analytic 
excitation  used  to  generate  the  sum  pattern  for  a monopulse 
array  antenna  radar  system  can  be  a well-designed  excitation 
which  produces  a broadside  narrow  beamwidth  radiation  pattern 
that  may  have  various  desired  characteristics.  Then,  based  on 
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the  relation  between  the  sum  and  difference  pattern  as 
indicated  in  Equation  (6.4),  it  is  easy  to  see  that  the 
characteristic  in  the  sum  pattern  will  always  be  preserved  in 
the  difference  pattern  which  is  excited  by  the  imaginary  part 
of  the  corresponding  analytic  excitation. 

Based  on  the  above  deduction,  let's  summarize  the 
properties  of  the  radiation  patterns  of  the  monopulse  array 
antenna  radar  system  using  a theoretical  analytic  excitation, 
that  is,  the  real  part  of  the  analytic  excitation  is  the 
excitation  for  sum  pattern  and  the  imaginary  part  for  the 
difference  pattern.  First,  both  of  the  sum  and  difference 
patterns  have  the  same  power  patterns,  that  is  both  sum  and 
difference  patterns  cover  the  same  illuminated  volume  in  the 
space.  Secondly,  the  difference  pattern  has  an  infinite 
difference  slope  at  the  boresight.  Third,  the  characteristics 
of  both  sum  and  difference  patterns  are  simultaneously 
controllable  to  meet  the  system  requirments,  such  as 
directivity,  beamwidth  resolution,  and  sidelobe  level. 
Actually,  those  properties  form  the  most  desirable  patterns 
for  a monopulse  array  antenna  radar  system.  Thus,  we  have 
reached  the  following  important  assertion: 

Assertion  6.1  The  most  desirable  radiation  patterns  for 
a monopulse  array  antenna  system  can  be  generated  through  a 
ideal  analytic  excitation. 

The  sum  pattern  and  the  corresponding  most  desirable 
difference  pattern  of  a monopulse  array  antenna  radar  system 
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are  plotted  together  in  Figure  6.2,  where  the  excitation  for 
the  sum  pattern  is  uniform. 


u (rad j an) 

□ sura  pattern  + difference  pattern 


Figure  6.2  The  desirable  radiation  patterns  of  an 

amplitude  monopulse  array  antenna  radar  system, 
(sum  and  difference  patterns) . 


6 . 3 Analysis  of  Radiation  Patterns  of  a Monopulse  Array 
Antenna  Radar  System  using  Realizable  Hilbert 
Transformer 

As  stated  in  Assertion  6.1,  the  most  desirable  odd 
excitation  for  the  difference  pattern  of  a monopulse  array 
antenna  radar  system  can  be  obtained  from  the  even  symmetric 
excitation  for  the  sum  pattern  by  using  an  ideal  Hilbert 
transformer.  It  is  this  elegent  consequence  that  stimulates 
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us  to  investigate  the  properties  of  the  difference  pattern 
generated  from  a sum  pattern  using  a realizable  Hilbert 
transformer.  The  generation  of  an  odd  symmetric  excitation 
from  an  even  symmetric  excitation  by  a realizable  Hilbert 
transformer  has  been  discussed  in  Chapter  III,  where  it  was 
considered  to  generate  an  analytic  excitation  for  a broadside 
narrow  beamwidth  radiation  pattern.  However,  the  purpose  and 
results  here  will  be  totally  different. 

Let's  denote  the  even  symmetric  excitation  for  a desired 
sum  pattern  of  an  N=2M+1  elements  monopulse  array  antenna 
radar  system  as  Js(n).  Then,  according  to  Eguation  (6.3),  the 
normalized  ideal  excitation  for  the  difference  pattern  is 
(jj^n)}  where  J^n)  is  generated  from  Js(n)  through  an  ideal 
Hilbert  transformer.  However,  by  using  a window  function,  the 
realizable  odd  symmetric  excitation  for  the  difference 
pattern,  denoted  as  Jd(n),  can  be  obtained  and  expressed  as 

Jd(n)  = j{[Js(n)*h(n) ]w(n) } (6.6) 

where  h(n)  is  impulse  response  of  the  discrete-time  Hilbert 
transform,  and  w(n)  is  a typical  truncation  window.  As  usual, 
let's  denote  the  sum  pattern  corresponding  to  Js(n)  as  2(u) 
and  the  difference  pattern  corresponding  to  Jd(n)  as  A(u). 
Then,  by  applying  the  convolution  and  multiplicative  law,  the 
realized  difference  pattern  A(u)  can  be  expressed  as 


A(u)  = [E (u) sgn(u) ] *W(u) 


(6.7) 
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where  [-jsgn(u) ] is  the  spectrum  of  the  discrete-time  Hilbert 
transformer,  and  W(u)  is  the  spectrum  of  the  selected 
truncation  window.  Thus,  the  difference  slope  A * ( 0 ) of  the 
realizable  difference  pattern  can  be  evaluated  by 


A'  (0) 


dA  (u) 


du 


u=0 


d{  [2  (u)  sgn  (u)  ] *W  (u)  } 
du 


u=0 


(6.8) 

However,  according  to  the  Fourier  relations  between  the 
radiation  pattern  and  its  excitation,  the  evaluation  of  the 
difference  slope  in  Equation  (6.8)  is  equivalent  to  the 
following 


A ' (0)  = 


dA  (u) 


du 


d{S(j  ( [ J (n)  *h(n)  ]w(n)  })  e~iun} 


u=0 


du 


u=0 


= S(n{ [ Js (n)  *h (n) ]w(n) )) 
n 


(6.9) 


The  result  in  Equation  (6.9)  provides  a neat  formula  to 
calculate  the  difference  slope  of  a difference  pattern  derived 
from  a realizable  Hilbert  transformer.  It  is  said  that  the 
difference  slope  is  the  weight  sum  of  the  element's  number  and 
the  weight  is  determined  by  the  truncation  window  and  the 
excitation  which  is  the  Hilbert  transform  of  the  even 
symmetric  excitation  for  the  sum  pattern.  It  is  interesting 
to  note  that  the  best  truncation  window  for  maximizing  the 
difference  slope  is  the  uniform  rectangular  window. 
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Based  on  the  results  in  Equation  (6.7)  and  (6.9),  for  a 
given  sum  pattern  E(u) , the  characteristics  of  the  difference 
pattern  can  be  controlled  by  the  number  of  array  elements  N 
and  the  selected  trancation  window.  As  just  mentioned,  the 
best  window  to  maximize  the  difference  slope  is  uniform, 
however,  the  best  chosen  window  to  minimize  the  beamwidth  of 
the  difference  pattern  with  a controllable  sidelobe  level  is 
exactly  the  same  as  that  discussed  in  Section  3.3.2,  that  is 
the  truncation  window  should  have  narrowest  mainlobe  width 
with  a specified  sidelobe  level,  which  leads  to  the  Dolph- 
Chebyshev  window.  By  using  the  windowing  method,  we  will  show 
the  flexibility  and  the  properties  of  the  difference  pattern 
which  are  obtained  from  a practical  sum  pattern  by  using  the 
realizable  discrete-time  Hilbert  transformer.  Then,  from  those 
analysis  results,  a specific  difference  pattern  for  the 
monopulse  array  antenna  radar  system  can  be  created. 

6 . 4 System  Analysis  and  Design  Example 

In  this  analysis  example,  let's  assume  the  specific  sum 
pattern  to  be  the  Taylor  pattern  with  peak  sidelobe  30dB 

and  n=4  (refer  to  Example  2,  Section  4.6,  Chapter  IV),  due  to 
its  popularity.  First,  the  directivity  and  the  difference 
slope  of  the  difference  pattern  are  analyzed  as  function  of 
the  number  of  the  array  elements  N.  That  is,  the  number  of  the 
array  element  N is  treated  as  an  independent  variable.  The 
truncation  windows  used  to  generate  the  realizable  odd 
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excitation  as  given  by  Equation  (6.6)  are  assumed  to  be  the 
uniform  and  the  Dolph-Chebyshev  with  R=20dB  for  comparison. 
The  analysis  results  with  N=ll,  13,  15,  17,  19,  21,  23,  25 
are  listed  in  Table  6.1,  where  Dd  and  SLpeak  indicate  the 
directivity  and  the  peak  sidelobe  of  the  difference  pattern. 


Table  6.1  The  characteristics  of  the  difference  pattern 
as  functions  of  the  number  of  array  elements  N. 


number 

of 

elements 

truncation  window 

uniform 

Chebyshev 

(R=20dB) 

N 

A'  (0) 

Dd 

SJjpeak 

A'  (0) 

Dd 

SIjpeak 

(None) 

(None) 

(dfi) 

(None) 

(None) 

(dB) 

11 

17.53 

6.70 

-12.08 

12.01 

6.53 

-14.84 

13 

24.49 

7.90 

-12.28 

17.21 

7.63 

-14.92 

15 

32.61 

9.11 

-12.40 

23.37 

8.68 

-14.67 

17 

41.90 

10.31 

-12.45 

30.51 

9.69 

-14.59 

19 

52.34 

11.52 

-12.45 

38.62 

10.65 

-14.61 

21 

63.94 

12.73 

-12.48 

47.71 

11.58 

-14.56 

23 

76.70 

13.93 

-12.55 

57.79 

12.46 

-14 . 56 

25 

90.62 

15.14 

-12.53 

68.85 

13.33 

-14.43 

The  analysis  results  of  directivity  and  difference  slope  in 
Table  6.1  are  plotted  in  Figure  6.3(a)  and  (b)  respectively, 
where  a linear  interpolation  scheme  is  adopted.  Figure  6.3(a) 
shows  that  the  directivity  of  a difference  pattern  generated 
by  using  the  Hilbert  transform  method  is  linearly 
increasing  with  the  number  of  array  elements  N,  and  in  Figure 
6.3(b)  the  difference  slope  increases  more  prominently  than 
linear  relation  with  increasing  N.  The  result  in  the 
difference  slope  just  reflects  the  property  of  the  nonlinear 
distortion  decreasing  with  the  increasing  number  of  the 
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Nurrtoer  of  Array  Elements  (N) 

□ uniform  + Dolph-Chebyshev 


(a) 


Nurrber  of  Array  Elements  (N) 

□ uniform  + Dolpb-Chebyshev 


(b) 

Figure  6.3  The  directivity  and  difference  slope  of 
a difference  pattern  as  functions  of  the 
number  of  array  enements  N. 

(a) Directivity . (b) Difference  slope. 
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truncation  points  of  the  Hilbert  transformer.2  Therefore,  the 
larger  the  number  of  array  elements,  the  more  advantage  we  can 
obtain  in  the  difference  slope  of  the  difference  pattern  by 
using  the  Hilbert  transform  method.  Moreover,  the  selected 
truncation  windows  control  the  slopes  in  both  directivity  and 
difference  slope,  and  the  uniform  truncation  window  has 
advantages  in  both  pattern  characteristics,  because  the  ratio 
of  the  mainlobe  width  of  the  uniform  window  to  the  mainbeam 
width  of  the  assumed  Taylor  sum  pattern  is  smaller  than  that 
of  the  Dolpf-Chebyshev  with  R=20dB  to  the  Taylor  sum  pattern, 
and  also  the  uniform  window  is  the  best  choice  in  terms  of 
difference  slope.  However,  the  peak  sidelobes  of  the 
difference  patterns  shown  in  Table  6.1  do  not  vary  with  the 
number  of  array  elements  N.  They  are  leaked  and  therefore 
determined  by  the  peak  sidelobes  of  the  selected  truncation 
windows . 

Next,  we  can  further  analyze  the  manner  in  which  the 
truncation  windows  control  the  performances  of  the  difference 
pattern.  Let  the  specific  sum  pattern  remain  the  same  as 
before,  and  the  number  of  array  element  N is  set  to  21.  Then, 
the  uniform,  Dolph-Chebyshev  with  R=20dB,  25dB,  30dB,  and  35dB 
are  used  as  the  truncation  windows  for  generating  the  odd 
excitation  in  Equation  (6.6).  The  characteristics  of  the 

2Recall  that  the  impulse  response  of  the  Hilbert 
transform  is  h (n)  = (1-cos  (n7r) ) /nn  which  is  a decreasing 
function  of  (1/n) . 
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difference  patterns  corresponding  to  those  truncation  windows 
are  tabulated  in  Table  6.2,  and,  again.  Figure  6.4(a)  and  (b) 
plot  the  analysis  results  of  the  difference  slopes  and  the 
peak  sidelobes  of  the  realized  difference  patterns  using  a 
linear  interpolation  method,  where  the  characteristics  values 
corresponding  to  -13dB  sidelobe  on  the  abscissas  of  Figure 
6.4(a)  and  (b)  are  the  results  of  the  uniform  truncation 
window. 


Table  6.2  The  characteristics  of  the  difference  patterns 
as  functions  of  the  truncation  windows. (N=21) 


truncation  windows 

Dd 

A'  (u) 

^■kpeak 

uniform 

12.73 

63.94 

-12.48 

Chebyshev (R=20dB) 

11.58 

47.71 

-14.56 

Chebyshev (R=25dB) 

11.96 

37.27 

-18.98 

Chebyshev (R=30dB) 

11.59 

30.45 

-23.45 

Chebyshev (R=35dB) 

11.04 

25.66 

-27.98 

It  has  been  shown  in  Table  6.2  and  Figure  6.4  that  using 
the  different  peak  sidelobes  of  the  window  functions,  we  can 
control  the  peak  sidelobe  of  the  difference  pattern  generated 
by  the  Hilbert  transform  method,  provided  that  the  sidelobe 
of  the  original  sum  pattern  is  low  enough  (it  is  always  the 
case) . However,  the  advantage  of  the  lower  sidelobe  is 
obtained  by  sacrificing  the  diffence  slope.  Furthermore,  the 
effect  of  the  truncation  windows  on  the  directivity  of  the 
difference  patterns  seem  minute.  Therefore  the  directivity 
ofthe  difference  pattern  largely  depends  on  the  number  of  the 
array  elements  N. 
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Figure  6.4  The  difference  slopes  and  the  peak  sidelobes  of 
a difference  pattern  as  functions  of  the 
truncation  windows. (N=21) 

(a)The  difference  slopes.  (b)The  peak  sidelobes. 
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From  the  analysis  result,  it  is  obvious  that  the 
generation  of  the  difference  pattern  from  a sum  pattern  by 
using  the  Hilbert  transform  exactly  provides  a flexible  and 
meaningful  method  to  control  the  performance  indices  of  the 
difference  pattern.  Normally,  the  number  of  the  array  elements 
is  determined  by  the  system  requirements  such  as  beamwidth 
resolutions,  and  array  size.  A sum  pattern  can  then  be 
specified.  Thus,  given  the  number  of  array  elements  N and  a 
specific  sum  pattern  and  by  applying  the  analysis 
proceduresdemonstrated  above,  it  is  easy  to  establish  a set 
of  data  shown  in  Figure  6.4,  which  is  relevant  to  the 
difference  slopes  and  sidelobe  levels  of  the  difference 
patterns  by  using  the  Hilbert  transform  method.  Thus  a 
difference  pattern  which  satisfies  the  system  requirement  that 
involves  a compromise  between  the  difference  slope  and  the 
sidelobe  level  can  be  obtained.3  As  an  immediately  practical 
example,  let  the  array  number  N and  the  sum  pattern  be  the 
same  as  those  just  used  in  the  analysis  example  above,  and  a 
23  dB  sidelobe  of  the  difference  pattern  assumed  to  be  another 
system  requirement  (ignore  the  difference  slope) . Then, 
checking  from  Figure  6.4(b),  a Dolph-Chebyshev  window  with 
R=30dB  can  be  used  to  generate  the  odd  excitation  as  given  by 


3It  is  possible  to  simultaneously  satisfy  the  system 
requirements  on  both  difference  slope  and  sidelobe  level  of 
the  difference  pattern  by  increasing  the  number  of  the  array 
elements . 
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Equation  (6.6)  for  the  required  difference  pattern.  After 
carrying  out  these  computation  routines,  the  realized  sum  and 
difference  patterns  are  plotted  together  in  Figure  6.5,  and 
the  corresponding  excitations  for  the  realized  sum  and 
difference  patterns  are  tabulated  in  Table  6.3.  The 
performance  of  the  realized  difference  pattern  is  shown  in 
the  fourth  row  of  Table  6.2. 


u (radian) 

□ Taylor  sum  pattern  + difference  pattern 


Figure  6.5  The  realized  sum  and  difference  patterns 
using  the  Hilbert  transform  method. 
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Table  6.3  The  excitations  for  the  realized  sum  and 


difference  patterns  of  a monopulse  array 
antenna  radar  system. 

(where  Js (-n)  =JS (n)  , and  Jd (-n)  =-Jd (n) ) . 


array 

number 

sum  pattern 
Js(n) 

difference  pattern 
Jd(n) 

0 

1.00000 

.00000 

1 

.98578 

. 13523 

2 

.94416 

.25555 

3 

.87801 

. 33899 

4 

.79138 

.39067 

5 

.68929 

.39124 

6 

.57821 

.36785 

7 

.46691 

.30392 

8 

.36666 

.24389 

9 

.29010 

. 16604 

10 

.24846 

.22614 

It  should  be  clear  that  the  design  methodology  presented 
in  this  chapter  for  the  linear  array  monopulse  radar  system 
can  be  easily  extended  to  the  planar  array  system.  The 
difference  pattern  in  either  u plane  or  v plane  can  be 
obtained  by  using  the  Type-II  2-D  analytic  signal  and  the 
associated  Hilbert  transform  (Type  II)  which  are  defined  and 
derived  in  Chapter  II. 

6 . 5 Discussion  and  Summary 

The  conceptual  association  of  the  analytic  signal  with 
the  radiation  patterns  of  a monopulse  array  antenna  radar 
system  has  stimulated  a new  approach  to  design  radiation 
patterns  for  monopulse  array  antenna  radar  systems.  The 
concept  of  ideal  analytic  signal  represents  the  most  desirable 
excitations  (its  real  and  imaginary  parts)  for  the  sum  and 
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difference  patterns  for  the  monopulse  array  antenna  radar 
system.  Using  the  realizable  Hilbert  transformation,  the 
performances  have  been  analyzed,  and  the  design  procedures 
have  been  formalized  for  the  difference  pattern  of  the 
monopulse  array  antenna  radar  systems.  Furthermore,  we  claim 
that  this  concept  can  be  equally  applied  to  the  design  of  both 
linear  and  planar  monopulse  array  antenna  systems. 

Another  attractive  property  of  this  design  method  is  that 
once  an  excitation  for  the  sum  pattern  is  specified,  the 
excitation  for  the  difference  pattern  can  be  analytically 
obtained. 


CHAPTER  VII 


CONCLUSION 

There  are  three  major  results  and  a new  discovery  in  this 
dissertation.  First,  the  analytic  signal  concept  has  been 
extended  to  two-dimensional  signal  domain.1  Secondly,  the 
application  of  analytic  signal  concept  to  the  design  of  both 
broadside  narrow  beamwidth  array  antenna  systems  and  monopulse 
array  antenna  radar  systems  has  been  introduced  and 
extensively  investigated.  Third,  by  associating  the  analytic 
signal  concept  with  the  modulation  theory,  a new  generalized 
design  methodology  has  been  developed  to  generate  practical 
desirable  asymmetric  sidelobe  array  antenna  radiation 
patterns.  The  new  discovery  is  that  a design  concept  has  been 
shown  that  enables  the  equiripple  (minmax)  design  criterion 
to  be  applied  to  the  synthesis  of  asymmetric  shaped  beam 
patterns . 

In  Chapter  II,  the  concept  of  the  analytic  signal  in 
analog  signal  domain  has  been  extended  to  one-dimensional  and 
two-dimensional  discrete-time  signal  domains.1  Based  on  the 

^he  extended  analytic  signal  is  denoted  by  analytic 
signal  (where  the  tilde  is  used) . The  analytic  signal  is 
defined  by  Definition  2.1.  The  analytic  signal  is  not  analytic 
since  its  derivatives  do  not  exist  and,  consequently,  it  does 
not  satisfy  the  Cauchy-Riemann  conditions  given  by  (2.1). 
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definitions,  two  types  of  two-dimensional  analytic  signals  and 
associated  discrete-time  Hilbert  transfroms  have  been 
developed.  They  are  differentiated  as  Type  I and  Type  II. 
Furthermore,  there  are  three  interesting  results  related  to 
the  generation  of  two-dimensional  analytic  signals  either  by 
using  two  stages  of  transformations  or  by  using  one  and  two 
dimensional  discrete-time  Hilbert  transforms,  and  the  results 
are  formulated  by  three  Lemmas  in  Chapter  II. 

It  needs  to  be  emphasized  that  the  2-D  analytic  signals 
are  defined  from  the  application  points  of  view  and  the 
conclusion  in  Lemma  2.1  is  obtained  from  the  case  of  the 
separability  of  the  kernel  of  the  2-D  discrete-time  Hilbert 
transform.  As  for  further  study,  it  would  be  interesting  to 
deduce  a generalized  2-D  analytic  signal  and  to  derive  the 
properties  of  the  2-D  Hilbert  transform  whose  impulse  response 
may  not  be  separable  by  using  the  method  of  z-transfrom, 
although  this  could  be  a challenging  and  difficult  problem. 

The  analytic  excitations  (corresponding  to  the  analytic 
signals)  for  broadside  narrow  beamwidth  array  antenna  systems 
are  discussed  and  stated  by  three  Assertions  presented  in 
Chapter  III.  These  results  are  based  on  the  hypothesis  that 
an  ideal  Hilbert  transform  is  available.  By  using  windowing 
functions,  practical  design  procedures  are  developed. 
Radiation  patterns  resulting  from  the  use  of  a realizable 
Hilbert  transformer  are  examined  and  the  performance  is 
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evaluated,  where  the  measures  of  performance  include  half- 
power beamwidth,  first-null  beamwidth,  antenna  directivity, 
excitation  (radiation)  efficiency  and  sidelobe  level.  The 
results  show  that  the  array  radiation  patterns  excited  by  the 
realizable  analytic  excitations  are  inherently  asymmetric  and 
have  the  advantages  in  all  measures  of  performance  compared 
to  the  typical  broadside  symmetric  array  radiation  patterns, 
except  the  sidelobe  level.  The  cause  of  the  high  sidelobe 
level  in  the  realized  asymmetric  radiation  patterns  is 
attributed  to  the  realizable  Hilbert  transformer  which  can  not 
faithfully  transform  the  radiation  patterns  if  the  original 
patterns  possess  large  spectral  components  around  the  u=0 
direction. 

A new  generalized  design  methodology  for  desirable 
asymmetric  sidelobe  radiation  patterns  has  been  developed  for 
narrow  beamwidth  linear  and  planar  array  pattens  (Chapter  IV) 
as  well  as  shaped  beam  patterns  (Chapter  V)  . The  kernel  of 
this  design  method  is  the  philosophy  of  associating  the 
analytic  signal  concept  with  the  modulation  theory.  The 
excitation  generated  by  this  design  method  is  called  modulated 
analytic  excitation.  However,  there  is  a minor  diference  in 
the  design  procedure  used  to  generate  an  asymmetric  sidelobe 
broadside  narrow  beamwidth  and  the  shaped  beam  patterns. 

In  Chapter  IV,  the  purpose  of  the  modulation  process  for 
designing  asymmetric  sidelobe  narrow  beamwidth  array  radiation 
patterns  is  to  shift  the  mainbeam  of  the  broadside  array 
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pattern  off  u=0  direction,  so  as  to  minimize  the  distortion 
caused  by  a realizable  Hilbert  transformer.  By  adjusting  the 
modulation  angle,  such  that  the  first  null  of  a symmetric 
broadside  pattern  is  exactly  positioned  in  the  u=0,  we  can 
obtain  various  broadside  narrow  beamwidth  asymmetric  sidelode 
radiation  patterns. 

In  Chapter  V,  the  procedures  for  designing  an  asymmetric 
shaped  pattern  are  to  symmetrical ize  the  desired  shaped  beam 
pattern  and  to  shift  the  shaped  beam  region  off  u=0  direction, 
so  as  to  minimize  the  distortion  caused  by  the  realizable 
Hilbert  transformer  and  to  suppress  one  of  the  sidelobes.  It 
is  the  symmetrical izat ion  process  that  enables  us  to  apply  the 
preferred  equiripple  (minmax)  design  criterion  to  the 
realization  of  an  asymmetric  shaped  beam  pattern,  although  the 
realized  asymmetric  shaped  beam  pattern  may  not  obtain 
equiripple  result  because  of  the  process  of  Hilbert 
transformation . 

Finally,  in  Chapter  VI,  by  association  of  the  analytic 
signal  concept  with  the  required  patterns  of  an  amplitude 
monopulse  array  antenna  radar  system,  a new  design  method  has 
been  created  for  the  difference  pattern  of  the  monopulse  radar 
system.  That  is,  the  difference  pattern  can  be  generated  from 
the  sum  pattern  by  using  the  Hilbert  transform.  Based  on  the 
system  requirements  on  difference  slope  at  boresight  direction 
and  the  sidelobe  level,  the  design  procedures  have  been 
formalized. 
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In  summary,  new  design  methods  have  been  developed  using 
the  analytic  signal  theory  for  the  design  of  array  antenna 
radiation  patterns.  The  features  of  these  design  methods  are 
conceptual  simplicity,  and  computational  effectiveness  while 
yielding  satisfactory  results.  Furthermore,  these  design 
methods  allow  the  use  of  existing  theory  of  symmetric 
broadside  array  antenna  design  that  is  well  developed.  This 
makes  these  design  methods  even  more  straightforward  and 
attractive. 


APPEXDIX 


THE  DOLPH-CHEBYSHEV  LINEAR  ARRAY  SYSTEM 

The  Dolph-Chebyshev  linear  array  system  was  originally 
developed  for  optimal  synthesization  of  an  in-phase  excitation 
for  an  equal  interelement  distance  linear  broadside  array 
system.  It  is  optimal  in  the  sense  that,  once  the  sidelobe 
ratio  is  specified,  the  mainbeam  width  is  minimun,  or  if  the 
mainbeam  width  is  specified,  the  sidelobe  ratio  is  minimum. 
This  is  due  to  the  properties  of  the  Chebyshev  polynomial 
stated  as  follows.  A Chebyshev  polynomial  is  expressed  as 

Tn(x)  = cos(n  cos-1x)  for  |x|  < 1 

= cosh(n  cosh-1x)  for  |x|  >1  (A.l) 

Let  0= cos  -Lx,  it  can  be  recursively  expressed  as 

Tn+1(x)  = 2 xTn ( x ) - Tn.j(x) 

where  T0(x)=l  and  T1(x)=x.  Normally,  TN(x)  is  a Nth  order 
polynomial  of  x (or  cos <p)  , let  it  be  expressed  as 


N 

Tn(x)  = £ c (n)  xn  (A. 2) 

n=0 
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where  c(n)  is  the  coefficient  of  the  xn  term.  The  Chebyshev 
polynomial  of  TN(x)  is  characterized  by 

1.  The  polynomial  oscillates  between  ± 1,  for  |x|  <1 

2.  It  has  N roots  in  the  interval  -l<x<l,  those  roots 
are  located  at  xk=cos  0k,  0k={  [ (2k-l)  7r]/2n} , k is  from 
1 to  n. 

3.  It  has  N— 1 extreme  values (±1)  in  the  interval  -l<x<l, 
those  extreme  values  are  located  at  xm=cos0m, 
0m=(rn7r)/n,  m is  from  1 to  n-1. 

4.  Tn(x)  monotonically  increases  or  decreases 

as  xN  for|x|>l.  (depending  on  whether  N is  even  or  odd) 
The  pattern  of  an  2N+1  array  elements  (for  2N  array 
elements  follow  the  same  procedures  [19])  real  even  symmetric 
excitation  J(n),  denoted  as  F(u),  is  proportional  to 

N 

F(u)  = E J(n)cos(un)  (A. 3) 

n=0 

Note  that,  it  has  been  implicitly  assumed  that  the  weight  of 
the  excitation  in  the  phase  reference  element(n=0)  is  [2J(0)] 
in  representing  the  pattern  of  an  odd  number  of  array  elements 
as  given  by  Equation  (A. 3).  Let  the  pattern  in  Equation  (A. 3) 
be  alternatively  expressed  as 

N 

F(u)  = E J (n) cos [2n(u/2) ] (A. 4) 

n=0 

Note  that  the  following  two  equations  always  hold 
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cos(un)  = ( £ ) (cos(u)  )n_lc(sin(u)  )k,  for  k even  and  K < n 

(A. 5) 

and 

sin2u  = 1 - cos2u  (A. 6) 

Substituting  Equation  (A. 5)  and  (A. 6)  into  Equation  (A. 4), 
the  array  pattern  can  be  expressed  as 

N 

F(u)  = E d(2n) (cos2n(u/2) ) (A. 7) 

n=0 

where  d(2n)  is  function  of  array  excitation  J(n)  and 
represents  the  coefficient  of  the  cos2n(u/2)  term.  After  those 
mathematical  manipulation,  the  pattern  representation  in 
Equation  (A. 7)  has  exactly  the  same  form  as  the  Chebyshev 
polynomial  in  Equation  (A. 2)  with  the  order  of  the  polynomial 
increased  from  N to  2N,  except  the  independent  variables  x and 
cos(u/2).  Two  problems  remain:  we  need  to  transform  the 

appropriate  region  of  x into  cos(u/2),  and  to  equate  the 
coefficients  between  the  corresponding  two  polynomials.  For 
transformation  the  sidelobe  region  of  the  pattern  of  an  array 
system  is  identified  with  the  oscillated  region  of  the 
Chebyshev  polynomial,  and  the  mainbeam  region  is  associated 
with  the  monotonically  increasing  region.  The  Dolph's 
transformation  [19]  is  given  by  cos(u/2)  = (x/x0)  , and  x0  is 
chosen  such  that 


T2n(xo)  ~ R 


(A. 8) 
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where  R is  the  desired  ratio  of  mainbeam  maximum  to  the  peak 
sidelobe.  The  Dolph's  transformation  between  the  variables  of 
x and  u can  be  seen  as  follows:  as  u varies  from  -n  to  ir , x 
varies  from  0 to  1 in  the  sidelobe  region,  up  to  x0  at  the 
mainbeam  maximum  at  u=0,  and  then  back  through  1 and  other 
sidelobes  to  0 at  u=7r.  Once  R is  determined,  x0  can  be 
calculated  from  Equation  (A.l)  as 

x0  = cosh (cosh-1R/n)  (A. 9) 

Then,  subsituting  x=x0cos(u/2)  into  Equation  (A. 2)  with  N 
replaced  by  2N,  gives  the  pattern  of  the  Dolph-Chebyshev 
linear  array  system 

F (u)  = T2N(x0  cos  (u/2 ) ) (A. 10) 

By  equating  Equation  (A. 7)  and  (A. 10),  the  coefficient  d(2n) 
as  well  as  the  array  excitation  J(n)  can  be  calculated,  and 
the  corresponding  pattern  will  have  the  equal  peak  sidelobe 
with  the  ratio  R.  Also,  it  has  been  shown  that  the  mainbeam 
width  of  the  pattern  is  minimized.  For  a large  number  of  array 
elements,  the  coefficients  equating  procedure  is  tedious. 
However,  the  procedure  can  be  easily  carried  out  by  use  of 
computer  [76].  The  comparison  of  the  mainbeam  width  between 
the  Dolph-Chebyshev  and  the  uniform  linear  array  antenna 
systems  is  plotted  in  Figure  A.l,  where  the  number  of  array 
elements  is  21.  From  Figure  A.l,  given  a specific  sidelobe 
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ratio  R (in  unit  of  dB)  on  the  abscissa,  the  mainbeam  width 
of  the  corresponding  Dolph-Chebyshev  array  system  is  equal  to 
the  mainbeam  width  ratio  on  the  ordinate  times  (47r/N)  , where 
N is  the  number  of  array  elements  and  (47T/N)  is  the  mainbeam 
width  of  the  uniform  array  system.  For  instance,  with  the  same 
number  of  array  elements  and  a (-16.3dB)  sidelobe  ratio,  the 
corresponding  Dolph-Chebyshev  array  system  has  the  same 
mainbeam  width  as  the  uniform  array  system  (corresponding  to 
the  unity  of  the  mainbeam  width  ratio  in  Figure  A.l). 


Sidelobe  Ratio  R COB) 


Figure  A.l  The  mainbeam  width  of  the  Dolph-Chebyshev  linear 
array  system  with  various  peak  sidelobe  ratios  R 
compared  to  that  of  the  uniform  linear  array 
system. 
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